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Summary 



Solar Coronal mass ejections (CMEs) are large-scale ejections of plasma and 
magnetic field from the corona, which propagate through interplanetary space. 
CMEs arc the most significant drivers of adverse space weather on Earth, but 
the physics governing their propagation through the Heliosphere is not well 
understood. This is mainly due to the limited fields-of-view and plane-of-sky 
projected nature of previous observations. The Solar Terrestrial Relations Ob- 
servatory (STEREO) mission launched in October 2006, was designed to over- 
come these limitations. 

In this thesis, a method for the full three dimensional (3D) reconstruction 
of the trajectories of CMEs using STEREO was developed. Observations of 
CMEs close to the Sun (< 15 Rq) were used to derive the CMEs trajectories 
in 3D. These reconstructions supported a pseudo-radial propagation model. 
Assuming pseudo-radial propagation, the CME trajectories were extrapolated 
to large distances from the Sun (15-24Oi?0). It was found that CMEs slower 
than the solar wind were accelerated, while CMEs faster than the solar wind 
were decelerated, with both tending to the solar wind velocity. 

Using the 3D trajectories, the true kinematics were derived, which were free 
from projection effects. Evidence for solar wind (SW) drag forces acting in 
interplanetary space were found, with a fast CME decelerated and a slow CME 
accelerated toward typical SW velocities. It was also found that the fast CME 
showed a linear dependence on the velocity difference between the CME and the 
SW, while the slow CME showed a quadratic dependence. The differing forms 
of drag for the two CMEs indicated the forces responsible for their acceleration 
may have been different. Also, using a new elliptical tie-pointing technique the 
entire front of a CME was reconstructed in 3D. This enabled the quantification 
of its deflected trajectory, increasing angular width, and propagation from 2 
to 46i?0 (0.2 AU). Beyond 7 Rq, its motion was shown to be determined by 
aerodynamic drag. Using the reconstruction as an input for a 3D magnetohy- 
drodynamic simulation, an accurate arrival time at the LI Lagrangian point 



near Earth was determined. 

CMEs are known to generate bow shocks as they propagate through the corona 
and SW. Although CME-driven shocks have previously been detected indirectly 
via their emission at radio frequencies, direct imaging has remained elusive due 
to their low contrast at optical wavelengths. Using STEREO observations, 
the first images of a CME-driven shock as it propagates through interplanetary 
space from 8 Rq to 120 Rq (0.5 AU) were captured. The CME was measured to 
have a velocity of ~ 1000 kms~^ and a Mach number of 4.1ibl.2, while the shock 
front standoff distance (A) was found to increase linearly to ~2Oi?0 at 0.5 
AU. The normalised standoff distance (A/Do) showed reasonable agreement 
with semi-empirical relations, where Dp is the CME radius. However, when 
normalised using the radius of curvature (A/Rq), the standoff distance did not 
agree well with theory, implying that Rq was underestimated by a factor of 
~3-8. This is most likely due to the difficulty in estimating the larger radius 
of curvature along the CME axis from the observations, which provide only a 
cross-sectional view of the CME. The radius of curvature of the CME at 1 AU 
was estimated to be ~ 0.95 AU 
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Chapter 1 

Introduction 



In this chapter I will introduce the fundamental physics and concepts that are discussed in 
this thesis, beginning with a short introduction to the Sun, its various layers, atmosphere, 
and activity. This is followed by an introduction to coronal mass ejections (CMEs) which 
consists of a historical account of CME observations and their interpretation and a brief 
review of the modern perspective on CMEs. Finally, I discuss some of the open questions 
surrounding CMEs. 
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1. INTRODUCTION 



"Physicists say we are made of Stardust. Intergalactic debris and far-flung 
atoms, shards of carbon nanomatter rounded up by gravity to circle the sun. 
As atoms pass through an eternal revolving door of possible form, energy and 
mass dance in fluid relationship. We are Stardust, we are man, we are thought. 
We are story." 

Glenda Burgess (The Geography of Love: A Memoir) 

Since the dawn of mankind the Sun has been a source of great wonder and fascination. 
Some early civilisation worshipped the Sun as deity, indeed some of our oldest antiqui- 
ties such as Newgrange and Stonehenge are believed to be solar observatories of a form. 
Building these monuments required detailed knowledge of the Sun's motion relative to the 
Earth, and ever since then we have been increasing our knowledge of the Sun. As science 
and technology advanced we began to probe the Sun's structure, composition and energy 
source. This led to a picture of the Sun, in terms of stellar equations, as a hot ball of gas in 
equilibrium, with gravitational contraction being balanced by the energy released during 
nuclear fusion. The modern era of space borne observations has revealed the dynamic and 
active nature of the Sun in the form of sunspots, filaments, prominences, flares and coronal 
mass ejections (CMEs). The Sun is the power source for all life on Earth and will also 
ultimately be the cause of its death. 



1.1 The Sun 



The Sun, our nearest star, is a main sequence star of spectral type G2V. The Sun has a 
total luminosity Lq = (3.84 ±0.04) x lO^^W, mass Mq = (1.9889 ± 0.0003) x 10^° kg and 
a radius Rq = (6.959 ± 0.007) x 10*^ m (|Foukal| [2004| . As all stars, the Sun was born from 



a giant molecular cloud of approximate mass 10^ — 10^ Mq which began to gravitationally 
collapse and fragment. The process of collapse and fragmentation continued until one 
of these fragments attained a central temperature large enough to start hydrogen fusion. 



about 4.6 X 10^ years ago rtPrialnikl 2009^. At this point fusion taking place in the core 



produced enough energy to counterbalance the gravitational collapse. Currently the Sun 
is in a stable configuration, on the Main Sequence, where it is in hydrostatic equilibrium 
(VP = —pg)- The Sun will continue to maintain this stable state for about another 5 x 10^ 
years before entering the red giant phase. At this point the Sun will expand to about 100 
times its current size and begin shedding its outer layers, due to successive nuclear burning 
in ever more distant shells. This ultimately leads to the total loss of the outer envelope 
exposing a degenerate core, in which all nuclear burning has ceased, called a white dwarf 



(Phillips, 1995) 
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1.1 The Sun 



1.1.1 Solar Interior 

The modern picture of the Sun's interior structure has been buih up over time, the three 



most important contributions to this have been the 'standard solar model' (SSM; Bahcall 



et al. 1982), helioseismology, and solar neutrino observations. As we cannot directly ob- 
serve the interior of the Sun we have to model its structure and then compare the model 
to observed properties, iteratively changing the model parameters until they match the 
observations. The SSM is essentially several differential equations derived from fundamen- 
tal physics principles, which are constrained by boundary conditions (mass, radius, and 
luminosity). The SSM treats the Sun as a spherically symmetric, quasi-static system which 
is powered by nuclear reactions at its hottest part, the core. The system is assumed to 
start out as a cloud of primordial gas which collapses under gravity. The composition of 
this cloud and hence the Sun can only be altered by the process of 'nuclear burning' in the 
core. All the energy generated by the nuclear burning is transported by radiation except 
where convection is a more efhcient process. 

The process of 'nuclear burning' or nuclear fusion occurring in the core at about 1.5 x 
10^ K, is the power source of the Sun. The proton-proton (p-p) chain describes the basic 
process which occurs: 

\R + \H^lR + e+ + ue (1.1) 
?H + lH^^He + 7 (1.2) 
iHe + iHe ^ ^He + 2}H (1.3) 

where \IL is hydrogen, is an isotope of hydrogen (deuterium), |He is helium, gHe is an 
isotope of helium with one neutron, a positron, fg an electron neutrino and 7 a photon. 



The p-p chain splits off into three branches at (1.2) which operate simultaneously, their 
reaction rates are determined by the density, temperature, and elemental abundances in 
the core. The net result of the p-p chain is the fusion of 4 hydrogen nuclei to from a Helium 
nuclei or a-particle: 

4iH^|He (1.4) 

with the products weighing less than sum of the initial masses by 0.02866 amu or 4.8 x 
lO"'^^ kg. Using Einstein's mass energy relation E = mc^ this results in a energy output of 
4.28 X 10"^2 J (26.73 MeV) of which varyingly small proportions are carried away by the 



neutrinos (Foukal, 2004; Phillips, 1995). Comparing the total energy output of the Sun 
(Lq) to the energy released during one cycle of the p-p chain we can estimate the reaction 
rate to be 8. 97 X 10^7 g-i. The region where the temperatures and densities are high enough 
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1. INTRODUCTION 




Figure 1.1: At the centre of the Sun is the core (<O.25i?0) where temperatures reach 
^l.SxlO^K, high enough for fusion to take place. The energy generated at the core from 
the fusion process is transported towards surface via thermal radiation in the radiative zone 
(0.25-0.70 i?©). At this point the solar plasma is cool enough to from highly ionised atoms 
and becomes optically thick. As a result it is convectively unstable and energy is transported 
through mass motions in the convection zone (0.7-1.0 Rq). The visible surface of the Sun, the 
photosphere, is a thin layer in the atmosphere where the bulk of the Sun's energy is radiated, 
its spectra is well matched to a blackbody with peak temperature of 5,600 K. Above the Sun's 
visible surface lies the chromosphere and finally the corona where the temperature soars back 
up to 1-2 X lO^K. Image courtesy of Steele Hin.NASA/GSFC, 
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1.1 The Sun 



for fusion to occur is known as the core, and extends out to about 0.25 Rq (Figure [lT| ) . 
Outside of the core is the radiative zone (0.25-0.70 where thermal radiation is the 
most efficient means of transporting the intense energy generated in the core, in the form of 
high energy photons, outward. The temperature drops from about 7 x 10^ K at the bottom 
of the radiative zone to 2 x 10^ K just below the convection zone. As the radiation is 
thermal it is described by the Planck (blackbody) equation, the specific radiative intensity 
is given by: 



BxiT) 



2hc^ 



1 



1 



W 



m^ sr m 



(1.5) 



where h is the Planck constant, c is the speed of light, k is the Boltzmann constant. The 
peak of this function is given by Weill's (lispl3>ceiiient Xdnw \rnax — 

2.8977 X IQ-^T-imK 

thus the radiative zone is dominated by x-ray and gamma-ray photons. Due to the still 
high densities (2 x 10^-2 x lO^kgm"-^) in the radiative zone the mean free path of the 
photons is very small (~ 9.0 x 10~^ cm) hence it can take tens to hundreds of thousands 



of years to escape (Mitalas &: Sills, 1992). As the temperature continues to fall, highly 



ionised atoms begin to form once an appreciable number of atoms have formed the plasma 
becomes optically thick, and as a result, becomes unstable as indicated by the upper limit 
of the Eddington luminosity: 

KF<47rcGm (1.6) 
where k is the opacity, F the radiation flux, G the gravitational constant, c the speed of 



light and m the mass coordinate (Prialnik, 2009). The sudden increase in opacity above 



the radiative zone causes the plasma to become convectively unstable. In this region called 
the convection zone (0.7-1.0 ii©), convection is the dominant form of energy transport. If 
the mass motions are rapid enough to assume they are adiabatic, and radiation pressure 
negligible, the Schwarzchild criterion for stability against convection may be used: 



dT 



dr 



star 



7-1 

7 



dP 



dr 



(1.7) 



star 



where 7 = Cp/Cy is the ration of speciflc heats (Prialnik, 2009). This gives a lower limit 
on the conditions necessary for convection to occur. In other words for convection to occur 
the temperature gradient must be larger than the adiabatic gradient. 

Helioseismology allows us to probe the solar interior by studying the propagation of 
sound waves in the Sun. Solar pressure waves (p-modes) are believed to be generated by the 
turbulence in the convection zone near the surface of the sun. Only certain allowed modes 
(spherical harmonics) can persist, as a result the Sun 'rings' like a bell. As these acoustic 
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1. INTRODUCTION 




Figure 1.2: This plot shows the relative differences between the squared sound speed in the 
Sun inferred from two months of Michelson Doppler Imager (MDI) data and from the SSM. 
The feature at about 0.7 r/R corresponds to the tachocline. The horizontal bars show the 
spatial resolution, and the vertical bars are the error estimates. Image courtesy of [Stanford] 
Solar Oscillations Investigation (SOI) , 



waves travel, they are refracted due to the waves' speed dependance on temperature, hence 
depth. This means that for all allowed modes the waves will be propagating normal to 
the solar surface when they reach it. This motion of the wave can be detected as Doppler 
shifts at the surface. Different modes penetrate to different depths and by combining a large 
number, the entire solar interior can be studied. The comparison of the properties derived 
from helioseismology, and that of the SSM can be used to change the model parameters to 



better fit the data, and the modern comparisons are extremely good as shown in Figure 1.2 



During the fusion process a large number of neutrinos are produced and escape the 
Sun. These neutrinos can be detected on Earth and the flux compared to predictions 
from the SSM. When this was first done for the rare ®B neutrino flux (p-p III chain; 



Prialnik, 2009), there was a large discrepancy between the predicted flux and the observed 



6 



1.1 The Sun 
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Figure 1.3: Left: Shows the Sun's bipolar field. Middle: The magnetic field is being twisted 
by differential rotation. Right: Loops of magnetic field begin to break the surface forming 
sunspots. (From The Essential Cosmic Perspective, by Bennett et al) 



flux, which was only ~ 0.4 of the predicted value and led to the so-called 'solar neutrino 
problem' (Bahcall, 2003). Results from helioseismology strongly suggested that the SSM 
was correct, but the experimental neutrino results were rigorously tested and found to be 
correct, leaving a serious issue with either the standard model for particle physics, or with 
the SSM. The problem was resolved when it was theorised that neutrinos could oscillate, 
that is an electron neutrino {ug) could become a muon neutrino (t'^) as it propagated 
from the Sun to the Earth. This was a large step as it required the neutrino to have a 
small but finite mass. It became clear that the early experiment to measure the neutrino 
flux was only sensitive to electron neutrino (ve) and so could have 'missed' some of the 
flux. Today the measured neutrino flux (from the 3 flavours, electron, muon, and 
tau) is (5.44 it 0.99) x 10® neutrinos/cm^ s in good agreement with the SSM prediction of 



(5.05: 



-1.0\ 
-0.8^ 



lO'^neutrinos/cm^s (Bahcafl 2003). 



The core and the radiative zones of the Sun rotate rigidly (as a solid body) but the 
convection zone rotates differentially, there is a thin interface between the two regions 
known as the tachocline. Due to the meeting of the two bodies rotating at different rates 
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this region is subjected to large shear flows. These flows are believed to be the mechanism 
that generates the Sun's large-scale magnetic field and powers the solar dynamo. The Sun's 
magnetic field is mainly dipolar and aligned to the rotation axis, thus each hemisphere has 
an opposite dominant polarity (Figure 1.3 left). The differential rotation of the convection 
zone winds-up this field. This large scale twisting which transforms poloidal field to toroidal 



field is know as the 0-effect (Figure 1.3 middle). As the field is twisted up the magnetic 
pressure increases and bundles of magnetic field lines (flux ropes) can become unstable and 
rise up in the from of loops. Due to solar rotation, the Coriolis effect twists these loop back 
towards north-south orientation reinforcing the original poloidal field, this is known as the 



a-effect (Figure 1.3 right) and completes the afi-dynamo. When magnetic loops become 
buoyant and rise up through the surface they are visible as sunspots on-disk and mark the 
footprints of large loops which extend into the solar atmosphere. In a given hemisphere the 
leading sunspot and trailing sunspot will have opposite polarities, this order is reversed in 
the other hemisphere (Hale's Law). Also the tilt angle of the sunspots pairs have a mean 
value of 5.6° relative to the solar equator (Joy's Law). Sunspots are known to migrate from 
high latitudes towards the equator over an 11 year cycle (Sporer's Law; see Figure 1.4). 
The net affect is an increase in opposite polarity field at the poles, ultimately the majority 
of the field will be oppositely oriented and the dipole will flip. This occurs every 11 years, 
thus a complete cycle takes 22 years (N to S to N). The activity of the Sun, in the form 
of active regions, flares, transient events, and other associated phenomenon, is modulated 



by this cycle (see Figure 1.4 lower). 



1.1.2 Solar Atmosphere 

The Sun's atmosphere is composed of all the regions above the photosphere. Until now 
we have referred to the photosphere as the visible surface of the Sun, it is in fact a very 
thin layer of the solar atmosphere. The solar atmosphere is usually separated into three 
regions, the photosphere, chromosphere and corona based on their density, temperature. 



and composition as shown in Figure 1.5 However, this separation is a simplification as 
the atmosphere is an in- homogenous mix of different plasma properties due to up-flows, 
down-flows, heating, cooling and other dynamic processes. The density of the plasma gen- 
erally decreases through these regions with increasing height. The temperature decreases, 
reaching a minimum in the chromosphere, then slowly rises until there is a rapid increase at 
the transition region which continues into the corona. This rapid increase in temperature 
leads to the so-called 'coronal heating problem'. 

An important parameter in describing the solar atmosphere is the plasma-/? term, the 
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Figure 1.4: The position of the sunspots in equal area latitude strips, averaged over a solar 
rotation with respect to time (top) . The butterfly pattern is clear as the decrease in the upper 
limit of sunspot latitudes with time, (bottom) The average sunspot area as a function of time. 
The 11 year modulation is clear in both of these plots. Image courtesy of NASA MSFC, 



ratio of the thermal to magnetic pressures: 



^=^ = ^. (1-8) 

Pmg -D /2^0 

where n is the number density and ^uq the permeability of free space. In the photosphere 



the plasma-/? is large and the plasma motions carry the field with them (Figure 1.6). 
Ascending into the chromosphere and corona the plasma-/? becomes small and the plasma 
is constrained to follow the magnetic fields. Continuing upwards the plasma-/? drops again 
and the magnetic field is advected out with the solar wind plasma flow and ultimately 



forms the Parker spiral (Figure 1.6) 
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Figure 1.5: A ID model of the electron density Nf. [cm~^] and temperature Tf. [K] profile 
thought the solar atmosphere from Gallagher (20001 after Gabriel & Mason (1982). Neutral 
atoms are present in the photosphere and chromopshere but the plasma is fully ionised in the 
corona due to the higher temperature. 



1.1.2.1 Photosphere 

The photosphere is the visible surface of the Sun and is defined as the height where the 
optical depth, at visible wavelengths, equals 2/3 (tsqoo ~ 2/3, I = loe""^). This is the 
mean optical depth at which the photospheric radiation is emitted or where the effective 
temperature {T^jf = 5776) and blackbody temperature of the photosphere match. This 
can be seen substituting B{T) = a/nT'^ and F = (^T^fj, where F is the total radiative 
flux, in the general solution of the radiative transfer equation: 

i?(r) = I (r + 2/3)- (1.9) 

4 TT 



which gives, 



aT' = l{T + 2/3)aT,'fj (1.10) 



implying that r = 2/3 (Foukal, 2004). The temperature drops from 6,400 K at the base of 
the photosphere to 4,400K at the top. The spectrum of photospheric radiation is that of 
a blackbody with a large number of absorption features, Fraunhofer lines, due the upper 
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Figure 1.6: Plasma-/3 in tiie solar atmospliere as a function of lieight for two magnetic 
field strengths of 100 G and 2500 G. The layers of the atmosphere are segregated by the dotted 
lines. The corona is the only region in which the magnetic pressure dominates over the thermal 
pressure, a low f3 plasma (Aschwanden 2006). 



layers of the atmosphere superimposed on it. The photospheric number density ranges 
from ~10^^-10^^ m~^ over the depth of the photosphere (500 km). 

One of the main observable features in the photosphere is granulation due to the con- 
vective motions. Granules are small-scale features made up of brighter regions isolated by 
darker lanes, this interpreted as the upflow of hot, bright material to the surface which 
then flows horizontally and cools, flowing back down in the dark lanes. The plasma-/3 is 
much larger than one throughout the photosphere which means the magnetic fields are tied 
to the flows. Typical granules are of the order of 1,000 km in diameter and have lifetimes of 
5-10 minutes with vertical flow velocities of hundreds to thousands kms^^. There are also 
larger scale flow patterns know as mesogranulation and supergranulation. Mesogranules 
are typically 7000 km in diameter, have lifetime of hours with vertical flows of the order of 
tens of ms~^. Supergranules are larger still at diameters of 3 x 10^ km, and consequently 
have longer lifetimes of days, they have large horizontal flows and smaller vertical flows 
of the order of 0.5kms~^. Sunspots are also found in the the photosphere they appear as 
darker regions due to their lower temperature (4,000 K) as convection is suppressed by the 
strong magnetic fields (kG). Sunspots play an important role in the activity of the Sun as 
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they are the source of solar flares and many CMEs. 
1.1.2.2 Chromosphere 

The chromosphere lies above the photosphere, the temperature initially decreases to a 
minimum of ~ 4,500 K before increasing to ~ 20,000 K with increasing height. It occupies 
a region approximately 2,000 km thick and with a density of about 10^^ m~^ but the mass 
density decreases by a factor of 10^. The structure of the chromosphere is split between the 



hot bright magnetic network and the cooler darker internetwork (Gallagher et al, 1999). 
Jet-like structures, called spicules, with diameters of hundreds of kilometers and attaining 
heights of tens of thousand of kilometers, with flows of the order of 30 km s~^ lasting 5-10 
minutes are ubiquitous. 

The source of the chromosphere's increasing temperature is not fully understood, but 
looking at the details we can infer some physical properties. The initial decrease in tem- 
perature is due to the decrease in the density of H~ ions, decreasing the plasmas' ability to 
absorb radiation from below, thus the temperature falls. Further out, some non-radiative 
form of energy is deposited, energy which ionises the hydrogen. The free electrons produced 
excite atoms which de-excite by line emission such as H-a, Ca ii and Mg ii. There is a 
balance between the energy input and radiative losses, forming a broad plateau in temper- 
ature at about 6,000 K. This balance is limited by the supply of neutral hydrogen, as this 
decreases the number of neutral or partially ionised atoms decreases, and the temperature 
rapidly rises. At about 20,000 K there is thought to be another plateau due to Lyman-a, 
but as height increases the ionisation of hydrogen increases and Lyman-a emission can no 
longer balance the energy input and the temperature rises rapidly. This transition marks 
the edge of the chromosphere and the start of the transition region. 

While the nature of the heating mechanism is unclear, from the observations it is 
clear there must be some form of energy deposition occurring. Neither radiation nor 
conduction can not be the source as the temperature is lower at the base of the lower 
chromosphere and photosphere than in chromosphere proper (and would thus violate the 
laws of thermodynamics). Mass motions are neither observed nor applicable since the 
chromosphere is in hydrostatic equilibrium. The most likely source of the energy (heat flux) 



is the dissipation of compressional or sound waves as proposed by Biermann ( 1946 ) and 



Schwarzschild (1948). In this paradigm, the convective plasma motions of the photosphere, 
launches sound waves into the chromosphere which travel upwards with little dissipation. 
As the density drops, the waves steepen and form shocks which rapidly dissipate energy, 
heating the chromosphere. This type of acoustic heating is not appropriate in the network 
regions where the strong magnetic field suppress the convective motions which drive the 
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waves. This led to the idea of Alfven wave heating, first introduced by (Osterbrock, 1961 ). 
Aflven waves are magneto- hydrodynamic waves which propagate along magnetic fields, the 
restoring force is provided by magnetic tension and the ion mass provides the inertia. Here 
the magnetic field itself is responsible for transporting and depositing the energy from the 
photospheric motions. This type of heating matches well with observations of plage and 
emerging flux regions, which both show strong heating, implying the heat flux is related 
to the magnetic fleld strength. 

Filaments are seen as dark channels in on-disk Ha observations often over active regions 
or as prominences when observed on the limb as bright features. Spicules which are jets 
of plasma are also observed on the limb, typically reaching heights of ~3,000- 10,000 km 
above the solar surface and lasting only ~5-15 minutes. The transition region lies be- 
tween the chromosphere and corona, here the temperature rapidly jumps (over 100 km) to 
above 1 MK. Above the transition region the magnetic fleld dominates and determines the 
structures. The high temperatures result in prominent emission from carbon, oxygen and 
silicon ions in the UV and EUV. 

1.1.2.3 Corona 

The tenuous, hot, outer layer of the atmosphere is known as the corona. The electron 
density of the corona ranges from ~10^^m~^ at its base, 2,500 km above the photosphere. 



to ^10^^ m ^ for heights >1 Rq (Aschwanden, 2006). The density varies depending on the 



feature, the open magnetic structures of coronal holes can have densities in the region of 
(0.5-1.0)xl0^^m~^, streamers (3-5)xl0^^m~^ while active regions have densities in the 
region of 2 x 10^^ - 10^^ m~^. The temperature in the corona is generally above 1 x 10^ K but 
again varies across different coronal features. Coronal holes have the lowest temperature 
(less than 1 x 10^ K) followed by quiet Sun regions at l-2xl0^K, and active regions 
are the hottest at 2-6xl0^K with flaring loops reaching even higher temperatures. The 
high temperatures reached in the corona give rise to EUV and X-ray emission with highly 
ionised iron lines being a prominent feature. The visible corona during eclipses is due to 
Thomson scattering of photospheric light from free electrons in the coronal plasma. The 
corona has a number of components: 

• K-corona (kontinuierliches spektrum) is composed of Thomson-scattered photospheric 
radiation and dominates below ^2Rq. The scattered light is strongly polarised par- 
allel to the solar limb as a result of the Thomson scattering mechanism. The high 
temperatures mean the electrons have high thermal velocities which wash out (due to 
thermal broadening) the Fraunhofer lines, producing a white- light continuum. The 
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intensity of the K-corona is proportional to the density summed along the line-of- 
sight. 

• F-corona (Fraunhofer corona) is composed of photospheric radiation Rayleigh-scattered 
off dust particles, and dominates above ~2i?0. It forms a continuous spectrum with 
the Fraunhofer absorption lines superimposed. The radiation has a very low degree 
of polarisation. The F-corona is also know as Zodiacal light and can be seen with 
the naked eye at dawn or dusk under favourable conditions. 

• E-corona (Emission) is composed of line emission from visible to EUV due to various 
atoms and ions in the corona, containing many forbidden line transitions, thus it 
contains many polarisation states. Some of the strongest lines are Fe xiv 530.3 nm 
(green-line; visible), H-a at 656.3 nm (visible), and Lyman-a 121.6 nm (UV). 

• T-corona (Thermal) is composed of thermal radiation from heated dust particles. 
It is a continuous spectrum according to the temperature and colour of the dust 
particles. 



The earliest descriptions of the corona were based on a static model with a heat input 



at some level rg which is described by Tq, in a spherically symmetrical corona (Chapman 



& Zirin 1957). The goal then was to derive the density, pressure, and temperature with 



respect to this reference level. The equation of hydrostatic equilibrium is: 



dp 
dr 



-P- 



[l.U) 



where the density of the plasma is p = nnip and the pressure due to electrons and protons 
is p = InksT. Rewriting the pressure in terms of density, and substituting into ( |1.11[ ) 
gives: 



dp GMqitIp dr 
p IksT 



(1.12) 



and integrating: 



p{r) = pq exp 



GMQinip 

2kB 



dr 



(1.13) 



Due to the high coronal temperatures, conduction should play an important role. The 
temperature distribution is determined by the conservation of conductive flux q = kVT 
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where k is the thermal conductivity which in the absence of sources or sinks reduces to: 



V • « = 0. 



In the case of spherical symmetry, this can be written: 



(1.14) 



11 (,2^dZ 
dr \ dr 



0. 



:i.i5) 



this implies that: 



2, dT 
r K—— = constant. 
dr 



For a fully ionised hydrogen plasma to a good approximation we have k{T) 



(Spitzer, 1962) thus we can write (1.16): 



dr 



integrating yields: 



T{r) = To 



ro\2/7 



The integral from (1.13) may be evaluated: 

dr 



ro r 



^ 7_J_ r _ /ro\5/7 

27-g (^rsiy/'^ 5Toro V 



and from 1.13 pressure is given by: 



, GMQm„7 1 
p{r)=poew' 



and similarly for density: 



p{r) = Po (^) 



2kB 5 Toro V V r 



r / \ 2kB 5 Toro V ^ 1^ 



ro\2/7 



(1.16) 

(1.17) 
(1.18) 



:i.l9) 



(1.20) 



(1.21) 



Inspecting (1.21) and (1.20), as r — )• 00, the pressure tends to a constant value while 



the density goes to infinity which is clearly unphysical. Even ignoring the density issue, 
the pressure value the static model tends to is some seven orders of magnitude greater 
than pressure in the interstellar medium (ISM) which must provide a boundary condition 
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at large radii. 



1.1.3 Solar Wind and Heliosphere 

Eugene Parker was one of the next scientists to tackle the problem and ultimately solve 
it in the form of the solar wind. To explain the observation that comet tails always point 
away from the Sun, both when approaching it and receding from it, Biermann, in 1953 



suggested that there must be a continuous outflow from the Sun (Biermann, 1957 and 



references therein) . Parker (1958) was the first person to make the connection between 
Biermann's and Chapman's work: that the heat flow in Chapman's static corona could 
drive the stream of particles (solar wind) Biermann speculated must exist. 

The Parker model is a spherically symmetric, static, isothermal model of the corona 
with only radial flows. The conservation of mass equation (V • (pv) = 0) for this system is: 

|-(rV) = 0, (1.22) 

therefore 

r'^pv = constant. (1-23) 

The radial component of the momentum equation (p(v • Vv) = —VP + pg) can then be 
written: 

dv _ dp GMqp 

P ^j. fj^j. J.2 

Jv^_ldp_GA^_ (1.24) 
dr p dr 



As the model is isothermal we can write the equation of state p = 2pkBT/mp so ( 1.24 ) can 
be written by eliminating p using the equation of state: 

2kBTl\dv AkeT GMq 

(1.25) 



rup V J dr rm„ 



A critical point occurs when dv/dr — )• hence we define: 

r.-°^ (1.20) 



IlksT GMq 



m„ 2v1 
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Figure 1.7: The solar wind velocity v as a function of radius r for various values of the 
constant C. The five different classes of solution are indicated by the labels TV. cs is the 
sound speed Vc and rc is the sonic radius 7'c 



and rewrite (1.25): 



1 dv 
V dr 



(r - Tc) 



(1.27) 



Integrating (1.27) equation yields: 

2 



In 



4 In 



+ 4^ + C, 
r 



:i.28) 



Parker's equation. 



The various possible solutions to (1.28) depend on the constant C and are shown 



in Figure 1.7, While they are all mathematically valid most have unphysical properties. 



Solution I is double- valued which is unphysical as it implies the SW leaves the solar surface 
with sub-sonic velocities, reaches a maximum radius and returns to the surface at super- 
sonic velocities. Solution II is double-valued and never reaches the solar surface so is 
clearly unphysical. Solution III starts at super-sonic velocities at the solar surface which is 
not observed so this solution is neglected. Solution IV know as the 'solar breeze' because 
inspecting Figure 1.7 we see that as r — )• oo that u — )■ thus at large distances we may 



17 



1. INTRODUCTION 



approximate (1.28) with 



In 



4 In 



V 

Vr. 



r 



r^Vv = constant. 



(1.29) 



Inserting this into the mass continuity equation p = const/r^?; = const/r^Vc the density 
tends to a constant and as the plasma is isothermal the pressure must also be a finite 
constant. This is not physical as the extremely small ISM pressure cannot balance the 



SW pressure. Solution V crosses through a critical point at r 



Vn, V 



Vc called the 



sonic point. This point can be used to constrain the integration constant C = —3. Again 



inspecting Figure 1.7 we see as r — )• oo that v t>c so (1.28) is approximated by 



4 In 




(1.30) 



From the mass continuity equation the density is given by 

const const' 

P 



j.2y 



(1.31) 



so the density will tend to zero, and as the plasma is isothermal so will the pressure thus this 
solution can match the ISM pressure. So the Parker solar wind solution is given by solution 
V which starts off sub-sonically at the solar surface, the velocity increases monotonically 
with height, reaching the sound speed at the critical point (or sonic point), and propagates 
super-sonically thereafter. We can estimate the critical radius and velocity by assuming 
a typical temperature for the solar wind of T ~ 10^ K, deriving Vc ~ 120kms~^ and 
Tc = 12 Rq which corresponds to a solar wind velocity at 1 AU of 430kms~^ very close to 
the measured velocity of ~400kms~^. In 1959 the first indications of the presence of the 
solar wind came from the Russian Lunik III and Venus I spacecraft and was confirmed in 
1962, just four years after its prediction, by Mariner II spacecraft measurements analysed 



by Neugebauer & Snyder (1962). 



As the plasma-/? of the solar wind is much larger than one, (see Figure 1.6 ) the magnetic 
field will be advected outward with the solar wind as the field lines are 'frozen in' (see 



Section 2.1.3). The combination of this radial advection with the rotation of the Sun forms 



what is known as the Parker spiral. If we assume that solar gravitation and solar wind 
acceleration can be neglected beyond some distance rg, then the radial outfiow velocity 
[vr) can be approximated by a constant v. In a spherical coordinate system which rotates 
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with the Sun we can write the velocity components as: 

Vr = v, vg = 0, Vcf, = uj {r — ro) sinO (1.32) 

where w is the angular velocity of the Sun {uj = 2.7 x 10~^ rads~^). A differential equation 
for the velocity stream lines can be obtained from v x dS = 0: 



dr rdO r sin 



Vr V0 

Integrating this equation from tq to r and from to (j) gives: 



(1.33) 



il_l_ln(il) =^(<^-0o), (1.34) 



when r > tq this equation can be approximated by: 



(r-ro) «-(</>- 0o) (1.35) 

UJ 



which is in the form of an Archimedean spiral shown in Figure 1.8 



As the magnetic field 'frozen in' to the flow and considering that V • B = we may 
write: 

ro\2 



Br{r,e,(t))=B{e,(^ 
B^{r, d, 0) = B{e, 00) (^) - ^o) (^) ' sin ( 



(1.36) 



where B{9, (po) is the magnetic field at r = tq. The angle of the magnetic field with respect 
to the radial direction can be obtained from: 

tan^/;= ^ = (-) (r-ro)sin6' (1.37) 



Br \V 

and when r is large this is approximated by: 



tanV = ^ = f— ')sin0. (1.38) 



Br \ V 

Inserting typical values of the solar wind at Earth of v = 400kms~^, co = 2.7 x 10^^ 
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Sun O Start Time # End Time 

# Jupiter Saturn # Uranus # Neptune# Pluto 

Figure 1.8: The interplanetary magnetic field showing the Parker spiral geometry for a solar 
wind speed of 600kms~^ and a field line starting at a Carrington longitude of 0°. 



rads ^ and using r = 1 AU = 1.5 x 10^ km and sin(90°) = 1 we find tp ~ 45° close to the 
measured value (M. Goossens, 2003). 

It is now known that the solar wind consists of two components: the slow solar wind 
with typical 1 AU values for velocity, density, and temperature of ~400kms~^, ~10cm~^, 
~1.4xl0^K respectively; and the fast solar wind with typical 1 AU values for velocity, 
density, and temperature of 800kms~^, ~3cm~^, and ~1.0xlO^ K respectively originating 
from open magnetic field regions (see Figure |1.9| ) . The two different speed streams can 
interact to form co-rotating interaction region (CIRs) where the fast wind ploughs into the 
slow solar wind and can form shocks. Another feature of the solar wind is the heliospheric 
current sheet (HCS) which separates the two opposite polarities of the Sun's magnetic, 
field forming what is know as the 'ballerina skirt'. 

The solar wind cannot expand forever and eventually runs into the ISM at what is 
known as the termination shock, where the solar wind transitions back to sub-sonic veloc- 
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a Ulysses First Orbit 



b Ulysses Second Orbif 



Ulysses Third Orbit 




Figure 1.9: (a-c) Polar plots of solar wind speed with the magnetic field polarity indicated by 
colour for three Ulysses orbits. The background image are blended composites from SOHO/EIT 
195 A, the Mauna Load K coronameter and SOHO/LASCO C2. (d) Show smooth sunspot 



number and heliospheric current sheet tilt angle ( McComas et al. 2008 ) 



ities (Figure 1.10). The Voyager II spacecraft recently (and previously Voyager I) passed 



though this region at some 70-90AU, details of these shock crossings can be found in 



Richardson et al. (2008), Burlaga et al. (2005) and Decker et al. (2008). Outside of the 



termination shock lies the heliosheath where the ISM and SW are in pressure balance, its 
outer boundary is the heliopause which also marks the edge of the Heliosphere. In this 
region the interaction between the ISM and SW causes turbulence and heating of plasma. 
As the Sun travels around our galaxy and encounters the ISM a bow shock is believed to 



from ahead of the heliopause (Figure 1.10) 



1.2 Coronal Mass Ejections 



"We define a coronal mass ejection to be an observable change in coronal 
structure that (1) occurs on a time scale between a few minutes and several 
hours and (2) involves the appearance (and outward motion) of a new, discrete, 
bright, white-light feature in the coronagraph field of view. " 



Hundhausen et al. (1984) 



Coronal Mass Ejections (CMEs) are large scale eruptions of plasma and magnetic field 
which propagate from the Sun into the Heliosphere. A typical CME has a magnetic field 



strength of tens of nT, a mass in the range of 10^^-10^^ g (Vourlidas et al., 2002) and 
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Figure 1.11: The 3D structure of the hehospheric magnetic field. Figure courtesy 
of Steve Suess, NASA/MSFC. 
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Figure 1.10: Depiction of the Heliosphere and its prominent features. Shown are the spiraUing 
magnetic field hnes, termination shock, hehopause, and bow shock. Image courtesy of Steve 
Suess, NASA/MSFC. 
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is("ng/or STEREO 



1. What Is the question this proposal addresses? 

Solar coronal mass ejections (CMEs) are spectacular eruptions of plasma and 
magnetic fields that drive space weather in the near-Earth environment. Despite nearly thirty 
years of study, the basic physics that expels these plasma clouds into the solar system is still 
not well miderstood. The Smi-Earth Connection Coronal & Heliospheric Investigation 
(sECCHI ) aboard NAsA s recently launched Solar Terrestna] Relations Observatory 
(oTEREO ) IS designed to explore vanous manifestations of the CME process in tliree 
dimensions. Tlie aECCHI instrument coraplmient combines comparable solar disk, 
coronagraph. and heliosphenc observations from two distinct views and is well suited to 
explormg the physics ot CMEs. both at the source and dunuG their propagation to Earth 
(Howard el at. 2002). This therefore offers us a umque opportunity to investigate the 
detailed physics involved in imhahon and accelerahng CMEs. We propose to develop 
advanced image processing methods to extract the evolving morphology and kinematics oi 
CMEs and to compare these results with as yet unconfirmed predictions oj theory. 

I.I Theoretical Models 

From a theoretical perspective, several models have been proposed to describe the 
observed properties of CMEs shown in Figure I. The two-dimensional flux-rope mode] was 
first proposed by Priest & Forbes (1990) and subsequently developed by Isenberg, Forbes & 
Demoulin (1993) and Forbes & Priest (1995). In this model, the CME is assumed to be 
initially located at the centre of a bi-polar lield configuration as shown in Figure 2a. The 
field foot-point separation (cf\ is Liraduallv decreased, and an eruption is triggered by a loss 
of equilibrium or instability in the field. The flux-rope then accelerates away from the 
surface as overlying fields are sequentially disconnected from the surface by magnetic 
reconnection (Figure 2b-d). 

A more recent model tliat builds on the basic | 
features of these two-dimensional ideas, is the three- 
dimensional magnetic flux-rope model of krall et al. 
(2001) and Chen & Kjall (2003). They assumed that the 
kinematics of an erupting flux-rope could be described 
using the force-balance equation. 



— VP- 



0) 



where h is the height, and the terms on the right describe 
gas pressure, gravity and the Lorentz force. Dropping 
the small gravitational force and initially small drag, the 
acceleration can be expressed as 



d'h 



dl^ [Rln(8R/a^]]= 
where 0^ is the (poloidal) magnetic flux inside the tube, 
J'r the force in the radial direction, af the flux-rope 
radius, and R its radius of curvature (see Chen & Krall 
2003 for details). The flux-rope acceleration is therefore 
dependent on its geomefrical properties, including 




Figure 1.11: The typical components of a CME observed by the LASCO coronagraph on 
SOHO spacecraft. 



velocity between ~10- 2,000 km s ^ some times even reaching 3,500 km s ^ close to the 
Sun dYashiro et aI||2004D. At 1 AU, CME velocities (300 -1,000 km s"^) tend to be closer 



to the solar wind speed (Gopalswamy 2006 Lindsay et al, 1999: Wang et al, 2005). The 



energies associated with CMEs are of the order of 10 -10 J making CMEs the most 



energetic events on the Sun ( Vourlidas et al. . 2002 ) . Although CMEs often exhibit a three^ 



part structure which consists of a bright front followed by a dark cavity and bright core 
(see Figure [l.ll[ ), they may also exhibit more complex structures (Pick et al. 2006). In fact 



less than ~30% of CME events possess all the three parts (Webb & Hundhausen, 1987). 



CMEs are know to be the most important driver of adverse space weather on Earth and 



in the near-Earth environment as well as on other planets ( Prange et al. , 2004 Schwenn 



et al. , 2005 ) . The most famous phenomena associated with space weather is the Aurora 



Borealis or Northern Lights (see Figure 1.12). The Aurora is caused by energetic particles 



traveling along the Earth's magnetic field lines interacting with atoms (mainly nitrogen 
and oxygen) in the upper atmosphere producing emission (resonance or recombination). 
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1.2 Coronal Mass Ejections 



One of the most extreme space weather events in recent history occurred on 2 September 



1859 (Carrington, 1859). This event was associated with a white light flare observed by 
Carrington. About 18 hours after the flare a severe geomagnetic storm occurred causing 
widespread sightings of the Aurora down to latitudes as low as 18° North and the loss 



of a significant portion of the telegraph service for many hours ( Green et al. , 2006 ) . As 
society has become more reliant on technology it has also become more susceptible to 
the adverse affects of space weather. On the 13 March 1989 a geomagnetic storm caused 
large geomagnetically induced currents (GIC) which caused the failure of a transformer 
that ultimately led to collapse of the Hydro-Quebec network. This left some six million 
people without power for nine hours causing substantial economic losses upwards of $13.2 



million (Bolduc, 2002). The storm was caused by a CME ejected from the Sun on 10 
March 1989 impacting the Earth some ~ 50 hours later. Other CMEs have knocked-out or 
caused damage to satellites, most recently on 5 April 2010, when Galaxy 15 or 'Zombiesat' 
(Intelsat) stopped responding to ground commands. The satellite continued to broadcast 
while it drifted forcing other satellites to take evasive action to avoid interference from 
'Zombiesat'. Galaxy 15 was subsequently recovered and placed into a safe mode. More 
generally, the increased radiation due to space weather poses hazards to astronauts and 
passengers on long distance flights, especially those over the poles. Also, polar flights may 
be restricted due to communications blackouts caused by space weather. In the modern era, 
society's dependence on GPS, communication satellites and inter-connected power grids 
mean space weather is an increasing concern. The cost of a severe space weather event 
was estimated to be up to $2 trillion in a recent report by the National Research Council, 



( Committee On The Societal Sz Economic Impacts Of Severe Space Weather Events , 2008 ) . 



As such ,the monitoring and prediction of space weather is of the utmost importance to 
society. Understanding the formation, acceleration and propagation of CMEs is vital to 
this goal as they are the main drivers of adverse space weather. 



1.2.1 Historical Observations 



The earliest observation of a CME probably dates back to the eclipse of 1860 in a drawing 
recorded by G. Temple shown in Figure 1.13[ It took 113 years for the CME to be for- 



mally discovered, the flrst deflnitive observation being made by Tousey (1973) using the 



coronagraph on-board the seventh Orbiting Solar Observatory (OSO-7). Following this 
a number of space-borne coronagraphs recored numerous CMEs. While the flrst imaging 
observation of a CME was in 1973, it is now apparent that CMEs and their effects had 
been observed much earlier, in a number of different types of observations. For example. 
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Figure 1.13: Drawing of tiie 1860 eclipse recorded by G. Teraple from Torreblanca , Spain 



with what is probably a CME (Eddy 1974). 



geomagnetic disturbances caused by CMEs had been recorded as early as 1724. CMEs 
were also observed in radio observations via interplanetary scintillation (IPS) from the 
1960s, however it was not until 1980s that the IPS could be directly related to the CMEs. 
Also fast CMEs produce shock waves and these shock waves then produce radio emissions 
at the local plasma frequency. These are known as Type II bursts and lead indirectly to 
the first height-time plot of a CME or at least the shock it drove as shown in Figure |1.14[ 
These shock waves can also be detected in measurements of the solar wind plasma in situ. 
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Figure 1.14: A height-time plot derived from Type II and III radio burst from [Wild et al 



It became clear, soon after the sunspot cycle was first discovered, that there was a link 
between the sunspot cycle and geomagnetic activity. One of the first people to suggest that 
transient ejections of plasma from the Sun could be the source of geomagnetic storms was 



Lindemann (1919), and later Chapman and Ferro also suggested that coronal transients 



could account for the geomagnetic activity (Chapman & Ferraro, 1931a[ b 1932). At this 
time there were two main interpretations of these transients: flare induced transient ejec- 
tions of plasma from formerly closed regions which drag the magnetic field outward with 



them (Gold 1962 Piddington, 1958); and shocks formed in open regions by the rapidly 
increasing expansion speed of the solar wind caused by the flare which would not drag any 
magnetic field with them (e.g., Parker, 1963). 

There were many signatures of these interplanetary transients such as solar energetic 
particle (SEP) events detected at Earth (Forbush 1946). Wild et al. (1963) suggested a 
two stage mechanism was required to accelerate the particles up to the observed energies: 
an initial flare acceleration followed by further acceleration by an outwardly-moving fast 
magnetohydrodynamic shock. The radio signatures of such moving shock fronts were 
commonly detected as Type II radio bursts. 

IPS is caused by the diffraction of radio- waves, from distant pulsars, by inhomogeneities 
in the solar wind plasma: as the solar wind is moving it produces a time- varying intensity 
signal which can be detected at Earth. This technique led to the detection and study of 



many transient events from the 1960s onward (e.g. Houminer, 1973 Houminer & Hewish 
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1972). Correlations were found between these events and in situ measurement at Earth 



(Houminer & Hewish, 1974 Rickett 1975) and in geomagnetic activity indexes (Vlasov 



1982). The exact nature of the transient density variations and their relationship, or 



lack thereof, to coronal transients was unclear. Attempts were made to relate the IPS 



observations with solar source features (Hewish 1988: Tappin et al, 1983) but with limited 
success. 

Entry into the space age allowed direct measurements of the properties of these coronal 
transients. The in situ signatures of the outwardly propagating shocks were fist discovered 



in Mariner 2 observations (Sonett et al. 1964), and later in Vela 3 observations (Gosling 



et al. 1968). It was found that the plasma driving these shocks had different properties 



compared to that of the normal solar wind. The material often showed enhanced helium 



content (Hirshberg et al. 1972), and low proton and electron temperatures (Gosling et al. 



1973). These early in situ observations indicated the Parker model was incorrect, and 



that the disturbances were associated with previously closed field regions. In some cases, 
streaming super-thermal electrons were found, suggesting the field lines threading the 
plasma were either connected to the Sun at both ends, or disconnected from it entirely. 



Hundhausen et al. (1970) used observations of shock disturbances to estimate the mass 
and energy associated with a large shock to be ~ 10^'^ kg and ~ 10^^ 



The magnetic signature of these transient events was first observed by Burlaga et al. 



(1981) in observations from five different spacecraft. The signature was a smooth rotation 
of the magnetic field vector following the shock, which they called a 'magnetic cloud' 
(MC) citing earlier theoretical work. However, the link between this signature and the 
coronal transients was not made until a following paper (Klein & Burlaga, 1982), which 
also identified the basic combination of properties that defines a MC which are still used 
today namley; low temperatures, high magnetic field strength, and a smoothly-rotating 
magnetic field vector. 

Around this time, the first space-borne coronagraphs were appearing. This led to the 



first remote imaging observation of these coronal transients (Tousey, 1973). Soon after, a 



direct link between these transients and radio Type II bursts was made (Stewart et al. 



1974a|b ). The near continuous monitoring of the solar corona by various coronagraphs in 



the 1970s provided some unexpected results. Far more ejections were found than would 
have been expected on the frequency of occurrence of shock wave disturbances, and the 
ejections were far more commonly associated with eruptive prominences than impulsive 



flares (Gosling et al., 1974 Munro et al., 1979). However many were followed by long- 



duration soft X-ray events lasting hours (Sheeley et al. 1975; Webb et al. 1976). Also, 



'These values are remarkably close to the modern measurements of CME mass and energy. 
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surprisingly, most of the material in the ejection was of coronal origin rather than flare 



or prominence ejecta (Hildner et al. 1975). The ejections also showed a wide range of 



speeds from hundreds of kilometres a second to greater than 1,200 km s ^ (Gosling et al. 



1976), morphologies, and energetics (Howard et al. 1985). The unexpected results led to 



the so-called 'Solar Flare Myth'. The first use of the term 'coronal mass ejection' appears 



to be from Gosling et al. (1976), and likely arose from the more tentative 'mass ejection 



coronal transient' which first appeared in (Hildner et al.. 1975). 



1.2.2 The 'Solar Flare Myth' 

Up until the early 1990s, despite all the unexpected results, all the activity discussed above 
was primarily attributed to solar flares, and CMEs were believed to be the result of flare- 
driven shock waves. However the evidence against such an interpretation had been growing 
throughout the 1970s and 1980s. For example, CMEs and geomagnetic storms were often 
not associated with flares, and the energy required to launch the CME was greater than the 



flare (MacQueen, 1980). Also, coronagraphic observations demonstrated that the flanks 



of the CMEa did not move laterally as the loop top moved outward through the corona 



which was not consistent with the shock interpretation (Sime et al., 1984). Further, in a 



number of papers which back-projected CMEs to find their onset times, none were found 
to be coincident with flares and typically, the flare occurred some time after the CME onset 



(Harrison & Sime, 1989 Harrison et aZ. , 1990). The debate continued and intensified with 



publication of Gosling's 1993 paper entitled 'The solar flare myth' in which the author 
demonstrated that the source of interplanetary shocks and of most geomagnetic storms 
were CMEs not flares, and that the CME flare relationship was secondary at best. In this 
paper Gosling proposed a "modern paradigm" describing the relationship between flares. 



CMEs, and geomagnetic activity shown in Figure 1.15 



1.2.3 Current Understanding 

CMEs originate wherever flares and prominence eruptions occur. Both flares and promi- 
nences are associated with active regions (ARs), which are regions of high magnetic field 
with or without sunspots. The most energetic CMEs seem to come from ARs which contain 
sunspots of opposite polarity. Cool prominences suspended above ARs also often produce 
CMEs, prominences can also form above neutral lines between sunspots of opposite polar- 
ities in ARs. Studies of in situ measurements of CMEs at 1 AU indicate CMEs probably 
remain attached at the Sun. The counter-streaming (bi-directional) particles in CMEs 



indicate that the 'legs' remain anchored on either side of the neutral line (Farrugia et al 
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CAUSE AND EFFECT IN SOLAR-TERRESTRIAL PHYSICS 
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Figure 1.15: Gosling's interpretation of the relationship between flares, CME and geomag- 
netic activity. Note the flare is at best from a common process "evolving solar magnetic fields" 



(reconnection) or secondary to the CME launch process Gosling (19931. 



1993a|bt [Kahler k Reamesj |1991| >. So closed magnetic structures seem to be the basic 
characteristic of CME-producing regions on the Sun. Thus the energy required to power 
the CME must ultimately come from the magnetic field itself. 



The general appearance of a CME was shown previously in Figure 1.11, but not all 
CMEs share this appearance. CMEs are also associated with on-disk 'dimming' regions, 



typically on either side of the neutral line (see review by Sterling in Wilson , 2003 ) . Dim- 
ming is the reduction in the intensity of radiation due to physical changes in the plasma 



(mass motion, density, temperature) typically observed in X-rays (Bastian et al. . 1999) or 



EUV (Gopalswamy & Thompson, 2000). The long term observations of LASCO have also 



discovered a number of unclassified morphologies, flux ropes ( Chen et al. , 2000 Plunkett 



et al., 2000), prominence- less CMEs (Gopalswamy et al, 2001; Yashiro et al. 2003), and 
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Figure 3 | Kinematics and morphology of the 3D reconstruction of the 
CME front of 12 December 2008. (a) Shows the velocity of the middle 
of the CIVIE front with the corresponding drag model and, inset, the early 
acceleration peak. Measurement uncertainties are indicated by one 
standard deviation error bars. Ch) Shows the declinations from the ecliptic 
(0°) of an angular spread across the front between the CME flanks, with 
a power-law fit indicative of non-radial propagation. It should be noted 
that the positions of the flanks are subject to large scatter; as the CME 
enters each field of view, the location of a tangent to its flanks is prone 
to moving back along the reconstruction in cases in which the epipolar 
slices completely constrain the flanks. Hence the 'Midtop/Midbottom of 
Front' measurements better convey the southward-dominated expansion, 
(c) Shows the angular width of the CME with a power-law expansion. 
For each instrument, the first three points of angular width measurement 
were removed as the CME was still predominantly obscured by each 
instrument's occulter. 



The early acceleration phase results from the rapid release 
of energy when the CME dynamics are dominated by outward 
magnetic and gas pressure forces. Different models can reproduce 
the early acceleration profiles of CME observations, although it is 
difficult to distinguish between them with absolute certainty'^'^^. For 
this event, the acceleration phase coincides with a strong angular 
expansion of the CME in the low corona, which tends towards a con- 
stant in the later observed propagation in the solar wind. Although, 
statistically, expansion of CMEs is a common occurrence^', it is 
difficult to accurately determine the magnitude and rate of expan- 
sion across the two-dimensional plane-of-sky images for individual 
events. Some studies of these single -viewpoint images of CMEs 
use characterizations such as the cone model-"-' but assume the 
angular width to be constant (rigid cone), which is not always true 
early in the events'^'^*. Our 3D front reconstruction overcomes 



Figure 4 | In situ solar wind plasma and magnetic field measurements 
from the WIND spacecraft. From top to bottom, the panels show proton 
density, bulk flow speed, proton temperature and magnetic field strength 
and components. The red dashed lines indicate the predicted window of 
CME arrival time from our ENLIL with Cone Model run (08:09-13:20 UT 
on 16 December 2008). We observed a magnetic cloud flux rope signature 
behind the front, highlighted by the blue dash-dotted lines. 



the difficulties in distinguishing expansion from image projection 
effects, and we show that early in this event there is a non-constant, 
power-law, angular expansion of the CME. Theoretical models of 
CME expansion generally reproduce constant radial expansion, 
based on the suspected magnetic and gas pressure gradients between 
the erupting flux rope and the ambient corona and solar wind'^'^ ''. 
To accoimt for the angular expansion of the CME, a combination of 
internal overpressure relative to external gas and magnetic pressure 
dropoft's, along with convective evolution of the CME in the diverg- 
ing solar wind geometry, must be considered'^. 

During this early-phase evolution, the CME is deflected from a 
high-latitude source region into a non-radial trajectory, as indicated 
by the changing inclination angle (Fig. 3b). Although projection 
effects again hinder interpretations of CME position angles in single 
images, statistical studies show that, relative to their source region 
locations, CMEs have a tendency to deflect towards lower latitudes 
during solar minimum'''^". It has been suggested that this results 
from the guiding of CMEs towards the equator by either the mag- 
netic fields emanating from polar coronal holes''' or the flow pattern 
of the background coronal magnetic field and solar wind/streamer 
influences'''^'. Other models show that the internal configuration 
of the erupting flux rope can have an important effect on its propa- 
gation through the corona. The orientation of the flux rope, either 
normal or inverse polarity, wiU determine where magnetic recon- 
nection is more likely to occur, and therefore change the magnetic 
configuration of the system to guide the CME either equator- or 
poleward'". Alternatively, modelling the filament as a toroidal flux 
rope located above a midlatitude polarity inversion line results in 
non-radial motion and acceleration of the filament, because of the 
guiding action of the coronal magnetic field on the current motion". 
Both these models have a dependence on the chosen background 
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Figure 1.16: Top to bottom, the panels show proton density, bulk flow speed, proton temper- 
ature and magnetic field strength and components. The red dashed lines indicate a predicted 
window of CME arrival time. The magnetic cloud flux rope signature is behind the front. 



highlighted by the blue dash-dotted lines (Byrne et at. 2010). 



jet-like CMEs (Yashiro et al. 2003) which contain no aspect of the three-part structure. 

Much of the material contained in CMEs is already present in the corona so it is 
expected to have coronal temperatures and densities. This is not true of the cool (400- 



8,000 K), dense (10^^ -10^^ cm prominence material of the core. The cavity is believed 
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1 shock only 

2 shock + sheath 

3 shock* sheath + MC UMConly 
Figure 1.17: Possible spacecraft tracks through a fast shock driving CME, slow CME and 



corresponding features which would be sampled Gopalswamy (2006). 



to be of coronal temperature but lower density than the front and core. This is often 



interpreted as a flux rope geometry often seen in 1 AU in situ measurements, Figure 1.16 
shows an example. The decreased density, temperature, increased magnetic field strength 
and smooth rotation of the magnetic field are the classic signatures of a MC. Not all 
CMEs measured at 1 AU have MC structures, one possible explanation of this is that the 
spacecraft trajectory though the CME cloud only samples one part and misses the MC as 



shown in Figure 1.17 



CMEs have a wide range of speeds, accelerations and widths (see Gopalswamy, 2004 



for an overview). The kinematic evolution of CMEs is complex, consisting of multiple 



phases. Sheeley et al. (1999) first classified CMEs into two classes: (1) slow gradual CMEs 



associated with prominences; and (2) impulsive fast CMEs often associated with flares. 
Studies based on larger samples indicate that slow CMEs arise from prominence lift-offs 



and streamer blowouts, and fast CMEs from ARs and flares ( Gonzalez-Esparza et al. , 2003 



Gopalswamy & Thompson 2000 ) . Zhang & Dere ( 2006 ) ; Zhang et al. ( 2001 ) suggest a three 



phase model of initation, acceleration and propagation and relate this to other observable 
properties such as soft X-rays. Numerous studies have attempted to relate the observed 
projected kinematics to the true kinematics: dal Lago et al. (2003) found v^ad = 0.88?;exp; 
while Schwenn et al. (2005) found the 1 AU travel time followed Tar = 203.0 — 20.771nT;p_ 



These type of statistical studies only provide a rough estimate of the corrections (Vrsnak 



et al. 2007), which are not accurate enough for detailed comparisons to models. Detailed 



studies of single events are necessary to make such comparisons. 



One such study was conducted by Zhang et al. (2001) who found a strong correlation 



between the soft X-ray flux and CME accleration. Figure 1.18 shows the velocity evolution 
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Figure 1.18: Evolution of a CME velocity (LASCO) in terms of a three phases profile (dashed 
line) and the flare (GOES) soft X-ray flux temporal profile (solid line) for the event on i998 
June 11 (iZhang et al.\\2001l. 



for a CME, and the x-ray light curve: notice that the two vary nearly in step. Gallagher 



et al. (2003) conducted a similar study, but extended the data set back towards the Sun 



using UCVS and TRACE. They used a double exponential profile: 



a{t) 



1 



1 



arexp{t/Tr) adex.p{t/Td) 



(1.39) 



to fit the acceleration where a^, are the initial accelerations and r^, are the e- 
folding times for the rise and decay phase. Using this they found an acceleration peak 



of ~ 1,500 km s ^ which corresponded to the soft X-ray rise phase (Figure 1.19) and sup' 



ported a relationship between the X-ray flux and CME acceleration. This relationship. 



between X-ray flux and CME acceleration, was fully realised in a paper by Temmer et al. 
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Figure 1.19: The (a) height, (b) velocity, (c) acceleration profiles and (d) the GOES soft 
X-ray flux (Gallagher et al. 2003). The velocity was obtained by taking the first and second 



numerical derivatives: the first-difference values are also shown as filled circles. The solid line 
gives the best fit to the data using a double exponential acceleration profile (1.39). 



(2008). They compared the CME kinematics to the hard X-ray flux and found the peak 



acceleration corresponded to the peak X-ray flux Figure 1.20 This is consistent with the 
previous observation as the derivative of the soft X-ray flux is often used as a proxy for 
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Figure 1.20: Event kinematics from top to bottom, distance-time profile d{t), velocity v{t), 
and acceleration a{t) of the CME as observed by GOES SXI and LASCO and (bottom) the 
RHESSI 50 -lOOkeV hard X-ray flux of the associated flare ( [Temmer et~al] [2008]). The left 
panels show the full CME height range covered by the LASCO FOV the right panels zoom 
into the early acceleration phase of the CME. 



hard X-ray flux. The relationship between the peak acceleration and hard X-ray flux indi- 
cates a strong relationship between magnetic reconnection and acceleration. This could be 
interpreted as reconnection occurring in a current sheet behind the flux rope as suggested 



by a number of models (Section 2.2.1). Lin et al. (2010) tried specifically to test the CME 
kinematics against a number of models. They found peak accelerations of ~l,500kms~^ 
and ~600 kms~^ close to the Sun (< 3 Rq). The authors show it is difficult to give merit to 
one model above another, as within the inherent scatter they all reproduce the kinematics 
to a good degree. 

The above studies are subject to large uncertainties due to projection effects. A number 
of efforts have been made to overcome this limitation, such as forward modelling, tomo- 
graphic, and polarimetric techniques. Forward modelling uses a pre-assumed geometry 



such as the cylindrical model (Cremades & Bothmer, 2004), or cone model (Xie et al. 
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Figure 1.21: (A) Total brightness a CME occurring on 31 October 1998 at 4:56UT (front 
view), (B) a reconstructed side view of the CME in the (z, y) plane, (C) a reconstructed top 
view of the CME in the (x, z) plane, and (D) a topographical map of the CME displaying 
distance from the (x, y) plane ( Moran & Davila 2004 1 . The colour bar indicates distance from 
the sky plane in Rq. The solar disk is outlined. 



the 2D observations. The kinematics can then be derived from these best-fit models, how- 
ever they are subject to the large and unknown errors about the pre-supposed geometry. 



The polarimetric technique of Moran &; Davila (2004) uses the ratio of unpolarised 



to polarised brightness of the Thomson-scattered K-corona to estimate the average line- 
of-sight distance from the instrument plane of sky. An example of this reconstruction 



technique is shown in Figure 1.21 This technique is only applicable up to 5 Rq as beyond 
this the F-corona can no longer be considered unpolarised. Yet other studies attempt to 



tie on-disk or in situ signatures of CMEs together: Demoulin et al. (2008); Howard et al. 



(2008b). These studies have made progress but the fundamental problem of uncertainties 



in the derived kinematics still remain a problem. 



From comparing CME speeds close to the Sun to those at 1 AU (see Figure 1.22), and 



by studying the derived speed from radio observations (see Figure 1.23), it is clear that 
CMEs are accelerated as they propagate. This has been studied in terms of a 'drag' force 



(Cargill et a/., 1996 Tappin, 2006 Vrsnak Sz Gopalswamy, 2002) and an extended Lorentz 



force (Chen, 1996). The limited spatial nature of coronagraphic and in situ measurements. 



and both temporal and spatial limitation of radio and heliospheric imagers, mean that 
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Figure 1.22: Speed-distance profiles plotted on a log-log scale (Manoharan 20061. Average 



initial speed of the CME increases from bottom (~300kms~^) to top (~2500kms~^). The 
date and the start time of each CME are shown. These profiles converge toward the ambient 
solar wind speed at the Earth's orbit. The data points at i? < 30 i?© and R > 80 Rq are from 
LASCO and IPS measurements, respectively. 



little can really be said about the exact nature of the forces governing CME propagation. 
However, it is clear there is some force at play which tends to equalise the CME velocity 
to that of the background solar wind speed. 

The STEREO mission was proposed to remove or diminish many of the issues outlined 
above. The STEREO mission aims to achieve this by providing dual perspective views 
along the entire Sun-Earth line. It will facilitate the tracking of CME signatures close 
to the Sun all the way to Earth, and their association to the in situ measurements. The 
high resolution and cadence images from dual perspectives will allow the 3D trajectory 
of CMEs to be reconstructed over an extended range, allowing accurate kinematics to 
be extracted and compared to theory in an attempt to answer some of the fundamental 
questions outlined above. 
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Figure 1.23: The speed distribution of CMEs (left) and ICMEs or MCs (right) for a set of 59 



CME-ICME pairs (Gopalswamy 20061. The average of the distribution is indicated on each 
panel. 



1.2.4 Open Questions 

While our knowledge of CMEs has greatly expanded with improving observations and 
theoretical interpretations, new questions have been raised. Some of the most fundamental 
questions about CMEs still remain unanswered. How is the energy required to launch a 
CME built up and stored? While a number of models have been developed to answer this 



(see Section 2.2.1 ), there is still no consensus on the matter. Another key question is what 
leads to the release of this energy and the eruption of a CME? A related question is: are 
CMEs pre-existing structures, or are they formed during the eruption? Again a number 
of models attempt to answer these questions but, as of yet, still fall short. Some authors 



claim (Zhang et al. , 2001 ) that there are two (or more) types of CME: slow, gradual CMEs 
which are accelerated slowly over large distances; and impulsive fast CMEs accelerated low 
in the corona, often associated with flares. It is not clear if these two distinct types of CME 
exist, and are due to different processes, or if they are the extreme ends of a continuous 
spectrum of CME properties. 

Once CMEs leave the Sun they are accelerated by an interaction with the solar wind; 
very little is known about this interaction other than that it tends to accelerate most CME 
towards the solar wind speed while others seem to be unaffected. Attempts have been made 
to model this interaction as a form of drag, or by the action of an extended Lorentz force. 
A combination of limited observations, and uncertainties in what observations there are, 
have left this particular question wide open. Also at 1 AU why do some CMEs contain MC 
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structures, and others not, and if not, are they still connected to the Sun? This could be an 
observational effect, but it could also be a manifestation of different underlying structures 
and mechanisms at play in some CMEs. 

1.3 Thesis Outline 

The work undertaken for this thesis enhances the understanding of the kinematics and 
morphologies of CMEs as they propagate through the inner Heliosphere. Until now, stud- 
ies of these properties have been hampered for a number of reasons. The most significant 
of these are projection effects which affect not only the morphology but also the kine- 
matics. Also, the lack of contiguous high resolution and high cadence CME observation 
through the inner heliosphere has lead to significant uncertainties in the interplanetary 
CME kinematics. 

Chapter 2 discusses CME theory and related phenomena. Chapter 3 details CME 
observations, instrumentation, as well as the data reduction, and conversion to physical 
units used in this work. Chapter 4 details the methods used to analyse the observations. 
Crucial to this work are the 3D reconstruction techniques which are outlined as well as the 
drag modelling and fitting. The 3D reconstructions are the basis for the work undertaken 
in Chapter 5, where the trajectories of a number of CMEs were reconstructed and are 
presented. Chapter 6 presents results from drag modelling efforts on a number of CMEs. 
Also presented arc the results from a new 3D reconstruction method called 'elliptical tie- 
pointing' and the drag modelling results from this. In Chapter 7, a CME-drivcn shock is 
studied and compared with semi-empirical shock relations. Finally, Chapter 8 presents the 
main conclusions of the thesis and details possible future work that could follow on from 
these new developments. 
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Chapter 2 

Coronal Mass Ejections and 
Related Theory 



In this chapter the details of the theoretical framework used to describe CMEs and CME 
related phenomena are presented. We have seen in the previous chapter that magnetised 
plasmas are ubiquitous on the Sun and throughout the Heliosphere. The complex interac- 
tion of these plasmas m,ay be understood in terms of magnetohydrodynamics (MHD). The 
application of this theory to the various aspects of CME evolution is presented. Some of the 
current CME models associated with initiation, acceleration, and propagation are reviewed. 
The relevant theoretical description of shocks and its relation to CME-driven shocks is also 
presented. 
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2.1 Magnetohydrodynamic (MHD) Theory 

Magnetohydrodynamics combines Maxwell's equations with the theory of fluid mechanics 
in an attempt to describe the interplay between the plasma flow's effect on the magnetic 
field, and the altered magnetic fields effect on the plasma fiow. Simply put, moving charges 
(plasma flow) will generate a magnetic field, this field will affect any other changes, altering 
their motion and thus a feedback between the fiow and magnetic field is set up. In MHD 
this feedback is described by the induction equation. 

2.1.1 Maxwell's Equations 

Maxwell's equations form a closed set of equations and describe the interaction of magnetic 
(B) and electric (E) fields. They can be written in a derivative form in vacuum as follows: 

1 5E 

VxB = ^oj + ^^ (2.1) 

V-B = (2.2) 
(9B 

VxE = -- (2.3) 

V • E = -/) (2.4) 

where j is the current density, p is the charge density, (j,q is the magnetic permeability of 
vacuum, eo is the permittivity of free space and c is the speed of light. If the typical plasma 



velocities are much less than the speed of light, the second term in Ampere's law (2.1 ) may 
be neglected giving: 



V X B = /ioj. (2.5) 



2.1.2 Fluid Equations 



The fiuid equations, or Navier-Stokes equations, arise from applying Newton's Second Law 
[F = ma) to a fiuid (continuum), and can be written in the most general form (Cauchy 
momentum equation; momentum conservation): 

p{^+^- = + V • T + f (2.6) 

where v is the velocity vector, p is the pressure, T is the stress tensor and f represents 
other body forces. Using the convective derivative {D/Dt = d/dt + v • V), makes it more 
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clear that it is just a statement of Newton's Second Law: 



Dv 

p— = -Vp + V-T + i. (2.7) 

Generally, the mass continuity equation (mass conservation): 

^+V-(pv) = 0, (2.8) 

is used in conjunction with the fluid equations. This is an incomplete description, the stress 
tensor T is still unknown, however using knowledge of the viscous behaviour of the fluid, 
we may make some assumptions on the form and properties of T. For an incompressible 
Newtonian fluid the equations may be written: 

p(^^+^r■V^r^ =-Vp + nV\ + i (2.9) 

where // is the (constant) dynamic viscosity, and the mass continuity equation reduces 
to V • V = 0. The Navier-Stokes equations are strictly an expression of conservation of 
momentum and do not fully describe the system: additional information such as an energy 
equation, or an equation of state is necessary. An energy equation which is often used is: 

where L is the total loss function. Assuming the fluid is described by an ideal gas p = pRT, 
then for adiabatic processes £ = 0, thus (2.10) becomes: 

^ + V • Vp = -7pV • V (2.11) 

The Navier-Stokes equations can be non-dimensionalised by scaling each of the quantities 
by a characteristic measure: 

' ^ / V y , 1 

"=L' "=y' '=l'' ^=L^ 
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where L is characteristic length, V characteristic velocity, L/V characteristic time scale. 
Then, depending on the regime under study, the equation can be normalised by: 

p = —rjTjP, if viscous effects are small, (2.12) 
or by p = —rjp, if viscous effects are large. (2-13) 

The Navier-Stokes equations thus become: 

1^ + V . VV = -Vp' + iv'^v' + f (2.14) 



and: 



3? (^1^ + v' • W j = -Vp' + V'^v' + f (2.15) 

where 3f? is the Reynolds number (JR = pVL/p). In the low Reynolds number regime 
where viscous effects dominate, the pressure is provided by the viscosity of the fluid {p ~ 
pL/V); in the high Reynolds number regime where inertial effects dominate, the pressure 



is provided by dynamic or ram form [p ~ pV^). Equations (2.14) and (2.15) can be solved 
to give the force acting on a body in a hight 3? flow: 

F = ^CoApV^ (2.16) 

where Co is the drag coefficient, A is the cross-sectional area of the object perpendicular 
to the flow direction and a low flow or Stokes flow: 

F = GirpaV (2.17) 

where a is the radius of the object. The drag coefficient (Cd) is defined as the drag force 
normalised by pU'^ and the area A. 



2.1.3 MHD equations 



The ideal and resistive MHD equations are obtained by taking the adiabatic, inviscid 
(viscosity neglected) fluid equations, assuming that the only body forces are due to gravity 
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and the Lorentz force, namely: 



p,- + v.Vv 



dp 
di 



-Vp + j xB + pg 



+ V • Vp = — 7pV • V 



dp 
dt 



+ V • (pv) = 



these may be combined with Maxweh's equations through Ohm's law: 



j = c7(E + v X B) 



(2.18) 

(2.19) 
(2.20) 

(2.21) 



where a is the conductivity of the plasma. Ohm's law can be rewritten using Ampere's 



law (2.4) to obtain: 



E = -V X B H V X B, 

poa 



(2.22) 



and so Faraday's Law (2.3) becomes: 



aB 

~dt 



V X (v X B) - V X (?7V X B) 



(2.23) 



where rj = 1/ poa is the magnetic diffusivity. Assuming that r/ is constant (2.23) becomes: 

aB 



dt 



V X (v X B) - r?V X (V X B) 



using a vector identity this can be written: 

SB 



dt 



V X (v X B) - [V^B - V(V • B)] . 



Using the solenoid constraint V • B = we obtain the induction equation: 

aB 



dt 



V X (v X B) - r?V^B. 



(2.24) 



(2.25) 



(2.26) 



This equation is vital for any model that considers a magnetised plasma, as it describes 
how a magnetic configuration will respond to fluid motions and vice versa - which is more 
dominant depends on the ratio of the terms. In a similar fashion to the fluid equations. 
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we may define a dimensionless quantity Rm, called the magnetic Reynolds number to be: 

V X (v X B) , , 

Rm = \ (2.27) 

VL 

= , 2.28 

This can thought of as providing two different timescales (1) changes due to fluid motion 
dominating when Rm is large: 

B V ^ L 

~B Tmotion ~ TT (2.29) 



'^motion L V 

where t = V/L is the characteristic time scale, and (2) changes due to diffusion of the field 
when Rm is small: 

B Bo 

~ r/yg- =^ Tdif fusion ~ (2.30) 

'^dif fusion ^ V 



Let us consider a typical sunspot with L = 10 m, rj = 1 m^ s"^ and V = 10^ ms ^, so 
the magnetic Reynolds number will be: 

lO'^lO'^ 1, , , 

Rm = — — = 10^^ > 1, (2.31) 

so the sunspot is dominated by plasma motions and we can estimate the diffusion time 
scale to be: 

Tdif fusion = — = ^-[^ ~ 10^^ - 31, 000 years (2.32) 



In the case of a perfectly conducting fluid [Rm — ^ co) the magnetic field lines must 
move with the plasma, or are 'frozen in' to the plasma (Alfven, 1943) and the induction 
equation reduces to: 



5B 

~dt 



V X (v X B) 



(2.33) 



which forms the basis of ideal MHD. In this case, the magnetic field is tied to the plasma 
motions. 
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Sweet-Parker model 




Petschek model 




Figure 2.1: Geometry of the Sweet-Parker (top) and Petshek (bottom) reconnection models 
(Aschwanden 20061. The diffusion region has a length A and width S. 



2.1.4 Magnetic Reconnection 



It is clear from a comparison of the lifetimes of sunspots of days to weeks, and our estima- 



tion of the lifetime of a sunspot, thousands of years (2.32) that some other process must 



play a role in the dynamic evolution of magnetised plasmas. That process is magnetic 
reconnection. This is also indicated from the rapid energy release associated with CMEs 
and flares. Magnetic reconnection is generally defined as a change in the connectivity of 
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field lines with time, where the energy stored in the magnetic field is converted into par- 
ticle, thermal, and kinetic energies. In an ideal plasma, the field lines are 'forzen-in' and 
thus coupled to the plasma motion. However when regions of opposite polarity fiux are in 
close proximity a boundary layer with large currents must form to separate the two regions 



(this can be seen from (2.5)). The small but non-zero resistivity (or large but not infinite 
conductivity) in the boundary layer opposes these currents. At the very centre of the 
boundary layer the neutral layer the magnetic field must go to zero, so the field smoothly 
and continuously changes across the region. The total pressure in this region must balance 
on both sides, Bi + pi = pni = B2 + P2 where p is thermal pressure and B the magnetic 
pressure. Thus, in the neutral layer ,the plasma-/3 becomes large and non-ideal effects, 
such as diffusion can take place. This region is known as the diffusion region, grey area 
in Figure 2.1 This formalism was first developed by Sweet (1958), who showed diffusion 



could occur in the boundary layer, and then Parker (1957) who derived the scaling laws 
and together is known as Sweet-Parker reconnection. 



2.1.4.1 Sweet-Parker Reconnection 



The same MHD equations cannot be applied simultaneously to describe inside and outside 
the diffusion region. However, using some simplifications, separate solutions for the two 
regions can be found and matched. The Sweet-Parker model is a steady-state 2D model 
where the diffusion region is assumed to be much longer than it is width (A ^ 5; Fig- 
ure [2T] (top)). For a steady state incompressible fiow (V • v = 0), we can see that the 
inflow and outflow must balance = Vout^- Outside of the diffusion region, ideal MHD 



applies so there are no currents and the z-component of Ohm's law (2.21 ) is E^+VmBx = 0. 
Inside the diffusion region there are large currents and the flow has stagnated (w = 0), thus 
the z-component of Ohm's law is Ez = ijJz- For a steady state flow, this implies that the 
electric fleld inside and outside the diffusion zone must be the same by Faraday's Law 



(2.3). Thus we can write VinBx = rjJz or Vin = rjJz/Bx and integrate Ampere's Law (2.5) 



around the region to obtain B^ = fJ-Jz^, and thus we arrive at: 

V^n = 'l. (2.34) 


Now, if we assume the x-component of the magnetic field is completely destroyed, then by 
conservation of energy Bx/2^ = l/2pv'^^ and thus Vout = B^/ JJip = va and we can write 
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in inflow velocity in terms of the outflow velocity (or Alfven velocity): 

vfn = (vin^) (?) or vf^ = v'^^t ( -\ ) • (2.35) 



'A \ X m ~OUl 1 A 

The rate of reconnection is given by the ratio of the inflow to outflow velocities: 

( v^V'' ^ ( jlV'' 1 



M = ^ =[^] = ^ (2.36) 

where S = Lva/i] is the Lundquist number (or magnetic Reynolds number). The Lundquist 
number can be estimated from typical values of L = lO'^m and va = lO^ms"^ to obtain 
S 10^3 and hence M ^ lO"'^. 

For a typical solar flare/CME some ~ 10'^'^ J are released over a time scale of about 
10^ s, which means an energy release rate of ~ 10^^ Js~^. Assuming all this energy comes 
from magnetic reconnection we can write: 

dEM B^dV B\^ B^L^VA 

L Vout ~ 7, 7^- (2.37) 



dt 2^ dt 2^ 2^1 ./S 

where Em the magnetic energy. Taking a fleld strength of 10"'^ T we get an energy release 
rate of ~ 10^^ Js~^, which corresponds to a time of tens of days for a flare rather then 
minutes much too long compared to observations. Sweet-Parker reconnection is thus too 
slow, with the limiting factors being the ratio of the width to length of the diffusion region 
and the Alfven velocity. 



2.1.4.2 Petschek Reconnection 



Petschek (1964) put forward an altered version of the Sweet-Parker model, this altered 



conflguration is shown in Figure [23] (bottom). In this configuration the size of the diffusion 
region is very compact {5 ~ A) compared to Sweet-Parker (A >> 5). Due to the smaller 
size, the propagation time through the diffusion region is reduced and reconnection can 
proceed at a faster rate. However, as the size is reduced so is the amount of plasma that 
can flow through the diffusion region. The result of this is that much of the inflowing 
plasma turns around outside of the diffusion region and two slow mode shocks form due to 
the abrupt changes from to Vouf These shock waves are the main site where inflowing 
magnetic energy is converted to heat and kinetic energy. By assuming a potential field 



in the infiow region Petschek (1964) found the external field scales logarithmically with 



distance L. Using this he estimated the reconnection rate, at a distance L, where the 
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2 4 6 8 10 



Figure 2.2: The emerging parameter space for the various reconnection regimes. S is the 
Lunquist number, and Lgp is the length of the reconnection region normalised by density 
(Daughton & Roytershteyn 2011). 



magnetic field dropped to half its original value, to be: 

The reconnection rate only logarithmically depends on the Lundquist number (or Rm) 
using our previous estimate of S this gives M w 0.01 which is four orders of magnitude 
faster than Sweet-Parker. 

Both of the models discussed have no external driver, such as an inflow. They also 
require a potential field in the inflow region. Relaxation of these simplifications has led 
to a family of solutions. These generalisation of 2D reconnection can be summarised in 
terms of internal Alfven Mach number at the entrance to the diffusion region Mai, and 
the Alfven Mach number at the exterior inflow M^e jas: 



(^) = 1 - ^(1 - 6) [0.834 - In tan [\s-^ M^'' )\ (2.39) 



This contains the Sweet-Parker, Petschek (6 = 0) and other hybrid solutions (Aschwanden 
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20061 iPriest k Forbesl 120001 ISchrijver & Siscoel [20091) . Figure O shows how as the length 



scales (Lgp) and Lundquist number (S) vary, different regimes become important. 3D 
reconnection is significantly different from 2D reconnection in many regards for example 
in the 3d case x-point less reconnection can occur, as well as reconnection where the 
magnetic field lines have the same topology. Reconnection plays a vital role in many CME 
models and in the interpretation of observations, although not all models require magnetic 
reconnection. 



2.2 Coronal Mass Ejection Theory 



CMEs are known to be associated with flare and filament eruptions, but the exact mecha- 
nism of their driver is unknown. Additionally, there is a debate about whether the flux rope 
(which becomes the CME) is a pre-existing structure, or is formed during the eruption. A 
number of theoretical models have been put forth to describe the forces and mechanisms 
responsible for the initiation and acceleration of CMEs. In nearly all models some instabil- 
ity is the trigger for the eruption, these instabilities may be understood using the cartoon 



mechanical analogies shown in Figure 2.3 



Thermal Blast: Early models proposed that the rapid increase in thermal pressure 
caused by a flare cannot be constrained by the magnetic fleld and thus drives the 
CME outward. Observations suggest this is not correct as many CMEs have no 



associated flare, or the flare occurs after the CME (see Section 1.2.1 and 1.2.2) 



• Dynamo: These models require the rapid generation of magnetic flux by stressing 



of the magnetic fleld (Klimchuk, 1990). This driver, known as flux injection (Krall 



et oZ. , 2000), corresponds to one of the following cases; (1) pre-existing fleld lines 



become twisted; (2) new ring shaped fleld lines rise up and become detached from 
the photosphere; or (3) arch shaped fleld lines rise up while remaining rooted in the 
photosphere. The first case requires footpoint motion some two orders of magnitude 



faster then observed ( Krall et al. , 2000 ) . The second case has not been observed and 



there is no obvious force to lift the mass. The third case is the most plausible for 
this model. 



Mass Loading: The slow build up of mass which is then removed causes an eruption. 
In terms of observations it is consistent with growing quiescent and eruptive filaments. 
Theoretically, the magnetic energy pre- and post-eruption are compared to show the 



plausibly of the transition from a higher to lower energy state (e.g. Low, 1999). Two 
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1J 





Tether Release 




,.H<i'..S.„....^1i.f. 




Tether Straining 

Figure 2.3: Cartoons showing the mechanical analogues for different CME eruption mecha- 



nisms (Klimchuk 2001). 
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possible forms of mass loading are: (1) by prominences with dense cool material 
contained in a compact volume; (2) by relatively less dense material spread over a 
large volume, which is susceptible to the Rayleigh- Taylor instability. The first concept 
is supported by numerous observations of coincident filament and CMEs, where the 



derived mass of the filament is a critical constraint ( Zhang & Low , 2004 ) . The second 



concept is supported by low density cavities observed in some helmet streamers, but 



there are also many which do not contain low density cavities (Hundhausen, 1999) 



Tether Release: Magnetically dominated loops, such as those in the corona, remain 
stable due to a balance between the upward force of magnetic pressure (VB^/2/io) 
and the inward force of magnetic tension ((B- V)B//xo)- The field lines which provide 
the tension are sometimes called tethers. In the mechanical analogy, once a tether is 
released then the tension on the remaining tethers increases, this can continue until 
the stain becomes so large that the remaining tethers fail and the spring (CME) 
is launched upwards. A 2D model recreates this type of behaviour by driving the 
footpoints towards each other, until a catastrophic loss of equilibrium occurs and 



the x-point jumps discontinuously to a new height (e.g. Isenberg et al, 1993). In 
non- ideal MHD, reconnection at the x-point after the tethers are cut would launch 



the CME (e.g. Amari et al 2000 Lin k Forbes 2000 Mikic fc Linker 1999). 



Tether Straining: This is very similar to the tether release model except the strain 
is increased by some external force rather than a constant strain being redistributed 
among fewer and fewer tethers. A physical model of this is the 'magnetic breakout' 



model (Antiochos et al., 1999 Aulanier et al., 2000). This type of behaviour can 



also found using mass loading (Zhang & Low, 2002). The equilibrium loss model of 



Forbes & Priest (1995a) is also of the tether straining type, as are the sheared arcade 



model of Mikic & Linker (1994); Linker Sz Mikic (1995) and fiux rope models of Wu 



et al. (1995, 2000) 



The last three models are known as storage and release models where magnetic stress is 
slowly built up over time before the eruption, some of which is then rapidly released during 
the eruption. The tether straining and release models seem the most likely scenarios, as 
they are able to reproduce a number of observable features through 2D and 3D modelling. 
The full dynamical evolution of CMEs includes three phases: initiation, acceleration and 



propagation (Zhang et al. 2001). Below I will outline a representative number of models 



which describe the initiation and acceleration (Section 2.2.1) of CMEs and, separately, 
their propagation (Sectioning. 
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2.2.1 Initiation and Acceleration 



One of the initial problems facing storage and release CME models was the apparent 
paradox that the pre-eruption closed field system was of lower magnetic energy than post- 



eruption open field system (Sturrock et al, 1984). It was argued that the relaxation 



of the pre-existing field releases more energy than it takes to open the field or the free 
magnetic energy stored in the field should be greater than the energy required to open 



the field (Barnes Sz Sturrock 1972). Following this Kopp & Pneuman (1976) proposed a 



three stage model for eruptive fiares: (1) energy is stored in a force free arcade or flux 
rope; (2) the eruption occurs, fully opening the field; and (3) the open field is closed via 
reconnection to a nearly force-free state. The transition from 1) to 2) would be ideal and 2) 



to 3) non-ideal. Aly (1984) conjectured that this scenario was energetically unfeasible. He 
argued that the open field configurations must always have higher magnetic energy than 
corresponding force-free configurations, as long as the field is simply connected. In 1991 



both Aly (1991) and Sturrock (1991) developed proofs of Aly's conjecture. These results 



seemed to imply that CMEs were energetically impossible, however there are a number 
of ways to avoid this problem: the field may not be simply connected but contain x- and 
o-points, the field lines are not opened up to infinity or only a portion of the closed field 
lines are opened. 



2.2.1.1 Catastrophe Model 

In these models a 2D fiux rope catastrophically loses equilibrium and erupts due to foot- 



point motions in the photosphere ( 


Forbes & Isenberg 


1991 


Forbes & Priest 


1995b 


Isen- 


berg et al. 


1993 


Priest & Forbes 


1990 


) . Figure 


2.4 illustrates one such model where a 



flux rope is suspended above the photosphere due to an equilibrium between the upward 
magnetic pressure gradient force and the downward magnetic tension force. As the foot- 
point separation distance (A) is reduced, the flux rope moves downward due to increasing 
magnetic tension, and as the magnetic pressure grows the magnetic energy of the system 
also increases. When the critical footpoint separation is reached and equilibrium is lost, the 
pressure gradient force becomes larger than the tension force and the flux rope is propelled 
upwards. If all the kinetic energy is dissipated then the flux rope will stabilise at the upper 
equilibrium, as shown in Figure [2^ If not, the flux rope will oscillate between the critical 
height and some point above the upper equilibrium height. If reconnection occurs in the 
current sheet which forms below the flux rope, it will increase in height without limit. An 
analytical expression for the kinematics of the flux rope, provided it is thin (i? <C h), and 
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Figure 2.4: The theoretical evolution of a 2-D flux rope of radius Rq/Xq = 0.1 (Forbes & 



Priest 1995b ) . (b-c) The footpoints are slowly moved towards each other and flux rope height 
decreases, (d) At some critical distance Xq the flux rope loses equilibrium and abruptly jumps 
in height, (a) Shows the resulting height evolution as a function of the footpoint separation A. 
All coordinates are in units of Xq. 



before the formation of the current sheet h/Xo <^2 (Priest &: Forbes, 2000), is: 

h 



h'. 



VAO 



In 



XqJ 2 Xq 



A/2 



(2.40) 
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where h is the flux rope height, h is the velocity of the flux rope, ho is an initial perturbation 
velocity, Aq is the critical footpoint separation, and vao is the Alfven speed at /i = Aq. 
Simplified expressions may be obtained by separating the evolution into an 'early' (i <C 
^o/vAo, h/Xo <^ 1) phase: 



/ ■ \ 3/2 
v37r \ -^0 / 



(2.41) 



and a 'late" (/i/Aq ^ 1, In/i <C Ini?) phase: 



h 



VAO 













~ 2 



1/2 



(2.42) 



Once the current sheet forms, analytical solutions of the CME evolution can no longer be 
derived. 



2.2.1.2 Toroidal Instability Model 



This model consists of an extension of the 2D fiux rope to a 3D magnetic flux rope sus- 
pended above the photosphere, which is anchored in the photosphere at both ends Fig- 
ure [231 (p^ein p89l [19961 IChen et al] [20001 IKrall et al.\ [2000l [200T|. It is initially in 



equilibrium, and erupts as a result of poloidal magnetic flux being injected into the flux 
rope, triggering the toroidal instability (TI). There are two dominant forces acting on the 
flux rope: the outward Lorentz self- force (hoop force), and the inward Lorentz force due 
to the background field. The fiux rope is susceptible to an eruptive instability if the back- 
ground magnetic field (Bs) decreases sufficiently quickly and vanishes at infinity. A sudden 
increase in the poloidal flux would then lead the flux rope to rapidly and continually ex- 
pand, triggering a CME eruption. The condition on the background magnetic field for this 
to occur is: 



dlnBe: 
-fl— ^ > 3/2 
dr 



(2.43) 



where R is the major radius of the flux rope see Figure 2.5 (top). An equation of motion 
for this system may be written: 



M 



dZ 

~dt 



Fr + F„ + Fd 



(2.44) 
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Figure 2.5: (top) A 3D flux rope modelled as a current loop with subscripts "t" and "p" 
referring to toroidal and poloidal directions, (bottom) Geometry of the flux rope as a model 
for a prominence as indicated by vertical hash lines. The ambient field Bs marks the boundary 



of the flux rope (Chen 1996) 
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Figure 2.6: Toroidal instability simulation, (a) The height of the apex (thick), centre of 
mass (thin), and traihng edge (d ashed), (b) Corresponding speed, (c) The minor radius a{t) 



in 10^ km and w{t) in lOkms"! (Chen 19961. 



where Fji, Fg and F^ are radial, gravitational and drag forces respectively. The radial 
component can be written: 



Fr 



cm 



In 



8R 



1 



2 Bpa^ 



+ 2 



Br, 



1 + 



+ F, + Fd 



(2.45) 



where is the internal inductance, /3p is the ratio of the average gas to magnetic pressure 
at the flux rope boundary. This equation is too complex to find analytical solutions, but 



has been numerically solved an example is shown in Figure 2.6 



Kliem & Torok (2006) were able to derive an analytical solution based on a simplified 



model. They assumed only two counteracting Lorentz forces acted on the flux rope, ignored 
changes in the external field its effects on the flux rope, and assumed a simple profile for 
the external field {Bs{H) = BqH~^). They found that at the beginning the instability 
can be approximated by: 



h{T) = ^ sinh(Pir), /i = ^ - 1 « 1 



(2.46) 



58 



2.2 Coronal Mass Ejection Theory 



592 




19:50 20:00 20:10 20:20 UT 

Fiu. 8. — Nearly conslant-jcrk rise profili: for an unperturbed lorus-unstablc 
flux rope equilibiium with steeper field decrease above the flux rope than in Fig. 7. 
Thefielddecay index in this ease is n > 2.85, i.e., nearthe value for the farfield in 
the dipolar ease (see text for other parameter differences for aspect ratio and initial 
torus depth). Solid lines show the rise profile of the apex point of Ihe magnetic 
axis as in Fig. 7; dashed hnes show the rise profile of a fluid element below the 
apex, initially at A = 0.65Ao. The simulation data for this lower fluid element are 
scaled to the rise profile of the 2005 June 16 filament eruption, and the resulting 
Alfven lime, Alfviin speed, and footpoint distance are given. 



ideal MHD instability (the helical kink in this case) and recon- 
nection similar lo the run shown in Figure 8, but with an initial 
velocity of ^0.01 Ta, the rise was clearly exponential. 

While all data in Figure 7 and the solid line in Figure 8 monitor 
the apex of the magnetic axis of the flux rope, the dashed lines in 
Figure 8 show the rise of a fluid element near the bottom of the 
flux rope, which is a likely location for the formation of filaments. 
Lying initially at Q.SSfiQ, it belongs to an outer flux surface of the 
rope. Although the flux rope in the simulation expands during 
the rise, both the axis and the bottom part show an approximately 
constant jerk and no significant timing differences between the 
acceleration profiles. 

3.1. Scaling Simulation to Observation 

Figure 8 presents a scaling of the simulation data to the rise 
profile of the 2005 June 16 eruption, determined in three steps. 
First, the time of the velocity minimum near IOta in the simula- 
tion is associated with the onset time, Iq, of the rapid-acceleration 
phase, 19:54:58 UT, as obtained in § 2. Second, the time t[ of 
maximum simulated velocity is associated with a time halfway 
between the final MLSO data points, which yields a substan- 
tially better match between the acceleration profiles than is 
found assuming that the acceleration ceased at or after the final 
MLSO data point. These two choices yield ta = 32.5 s. Third, 
the simulated and observed heights are matched at t[ , resulting 
in a length unit for the simulation of /?o = 44.4 Mm, an Alfven 
speed Fa = ho/TA = 1370 km s~', and a nonnalization value 
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Fit;. 9. — Side view of a TI simulation (see Fig. 8). The field lines of the torus 
are shown lying in a flux surface at half the minor torus radius. Sample field lines 
for die overlying field are also shown. The starting points in the bottom plane for 
the traced field lines are the same for all panels. The times (expressed in Alfven 
crossing times, as in Figs. 7-11) are 0, 20, 30, and 40, respectively. The motion of 
the loop apex marked by an asterisk is shown in Fig. 1 1 . 

for the acceleration of oo = VJ'^k- Figure 8 shows the observed 
heights on a linear scale, with derived velocity and acceleration 
data (based on central differences, with seven-point boxcar aver- 
aging to smooth the heights and velocities and five-point boxcar 
averaging for the accelerations). 

Both the rise of the magnetic axis of the flux rope ( Fig. 8, solid 
line) and the rise of a fluid element originally below the magnetic 



Figure 2,7: Side view of a 3D toroidal instability (Schrijver et al. 2008). The field lines of 



the torus are shown lying in a flux surface at half the minor torus radius. Sample field lines 
for the overlying field are also shown. The starting points in the bottom plane for the traced 
field lines are the same for all panels. The times (expressed in Alfven crossing times) are 0, 
20, 30, and 40, respectively. 
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Figure 2.8: Model for Kink instability the magnetic field under study is modeled by a force- 
free circular flux tube with the total current /, a pair of magnetic charges — q, q and a line 
current Jq (Titov & Demoulin 1999). 



where H is the height of the flux rope, Hq is the height of the flux rope at on-set of the 
instability, r is time normalised by Alfven time, Pq comprises initial parameters of the flux 
rope and Pi is associated with the external magnetic field profile. A distinctive feature 
of their expression and simulation was that the acceleration shows a fast rise and a more 



gradual decay. However, Schrijver et al. (2008) was able to demonstrate that the evolution 
of the instability can be changed by tuning the initial conditions. The hyperbolic height 
dependance can change to a polynomial form, which changes the acceleration profile to a 
more gradual one. One issue with these models is that they contain no twist as is often 



observed in CME eruptions (see Figure 2.7). Still, a number of studies have shown good 
agreement between simulations and observations (Krall et al. , 2000, 2001 ). Recent attempts 
have tried to establish the amount of energy available from the photosphere to create the 



instability. Using an ideal flux rope event Schuck (2010) found there was insufficient energy 
to launch the CME. 



2.2.1.3 Kink Instability Model 



The kink instability (KI) aries from a similar set up as the TI, except it involves the twisting 
of the photospheric footpoints (increasing the toroidal flux) of an initially potential flux 



bundle see Figure 2.8 (Fan & Gibson, 2003; Rachmeler et al. 2009; Torok & Kliem 2003 



2005; Torok et al., 2004). The twisting motion serves to both form the flux rope, and also 



to move it into a region of possible instability. The continual twisting of a flux rope of 
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finite length must lead to a KI. The threshold for such an instability depends on the ratio 
of the azimuthal to axial field components {B^/Bz), the length-to- width ratio of the flux 
rope (L/r), and the radial profile of the flux rope. It can be expressed in terms of the flux 
rope twist: 



LB^rp) 
roBz{ro) 



(2.47) 



where tq is a characteristic radius of the flux rope. For a straight cylindrically symmetric 
flux rope with fixed ends of uniform twist (Gold-Hoyle equilibrium), the threshold was 
numerically found to be ^gh = 2.497r (Hood & Priest, 1981). There are no analytical 



solutions for the kinematics of the KI model, but it has been studied extensively using 
numerical simulations. Once the flux rope achieves "supercritical twist" two current sheets 
form, the first helical wrapping the rising kinked flux rope, and the other a vertical sheet 



below the flux rope (Fan k Gibson, 2003, 2004 Titov k Demoulin, 1999). The formation 



of the vertical current sheet enables the eruption to proceed with, or without, reconnec- 
tion. The KI model recreates some of the phenomena observed in CME-related eruptions 



including soft x-ray sigmoids ( Kliem et al. , 2004 ) and the twisted structures which often 



appear in CMEs and prominences (see Figure 2.9). Significantly it also enbles the build 



up of substantial mass and can release up to 25% of the stored energy (Torok k Kliem 



2005 ) . However not all simulations of the KI form current sheets and so remain confined 



i.e., with no CME eruption (Fan, 2005), and the amount of twist required to produce an 



explosive eruption may be non-physical. 
2.2.1.4 Magnetic Breakout Model 

In the magnetic breakout (BO) model the CME eruption is triggered by shearing which 
causes reconnection between the overlying field and the multipolar field below it, as shown 



in Figure 2.10 (Antiochos et al, 1999 DeVore k Antiochos, 2005; Lynch et al. 2004 2008 



MacNeice et al. 2004). This configuration has four distinct flux systems: a central low- 



lying arcade straddling the equator (Figure 2.10 a-blue); two low-lying side arcades (one 



on each side of the central arcade Figure |2.10 a-green) ; and a large scale (polar) arcade 



|2.10 


a-green) ; 


gure 


2.10 


a-re 



the central arcade. Shearing concentrated at the equatorial neutral line causes the central 
arcade to rise, distorting the null point into an x-point (Figure 2.10| b). Continued shearing 



causes the central arcade to rise even more, stretching the x-point to form a current sheet. 
Reconnection occurring in the current sheet transfers flux from the overlying field and 



the un-sheared field to the side arcades (Figure 2.10 c), thus creating a passage for the 
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Figure 2.9: Simulation of the kink instability (Fan 2005). The black lines corresponds to the 
range of r = l-3.6i?0 and is not the full simulation domain. The marked times are in units 
of Rq/vao- 
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Fig. 1. — Left: Schematic showing the main four stages of topological evolution during the axisymmetric (2.5D) breakout CME scenario, (a) initial topology and the 
shear channel, (fi) Shearing phase required to energize the system and the onset of magnetic hreakout reconnection at the distorted null line, (c) Beginning of flare 
reconnection that starts deep in the shear channel and creates the magnetic flux rope, and (d ) Antibreakout reconnection pbase that describes the system relaxation and 
topological restoratioTL Right: Corresponding field line plots fiom the MacNeice et al. (2004) MHD simulation. See text ibr further details. 



processed, the side arcades come together and reconnect at the 
flare current sheet. These side-lohe fields tend to be much less 
sheared and their reconnection proceeds in a slower, less ener- 
getic manner corresponding to the "gradual" phase often observed 
in flares. When the {green region) side-lobe fields reconnect, they 
restore bolh the inner {light blue region) and overlying (red region) 
flux systems — denoted by the green arrows showing the side-lobe 
separatrix surfaces shrinking. Since this topological evolution is 
the reverse of the original breakout reconnection, one can think of 
the posteniption restoration of the original flux system configura- 
tion as "antibreakout" reconnection. The antibreakout reconnection 



eventually dissipates the radial current sheet and eventually re- 
stores the coronal null point. 

The empting flux rope structure continues to propagate through 
the background dipole field lines, and in our axi symmetric picture, 
this means continual breakout reconnection at the interface of the 
oppositely oriented flux rope configuration and the overlying 
field. However, there is an important topological distinction to be 
made between the overlying field that has yet to encounter the 
ejecla and that which has either been reformed behind the flux 
rope ( from antibreakout reconnection) or has recently encoimtered 
the flux rope and has undergone reconnection such that it now 



Figure 2.10: Schematic showing topological layout (left) and evolution of breakout model 
(Lynch et al. 2008). (a) Initial multipolar topology, (b) shearing phase which energises the 



system causes magnetic breakout reconnection at the distorted null line, (c) flare reconnection 
occurs low down and forms the flux rope, and (d) the system relaxing after the eruption. 
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Figure 2.11: Running difference image from the breakout model density data (Lynch et al. 



2004j). The arrows indicate the leading edges of the bright CME front, the dark cavity, and the 



central core region. The time stamp indicates the date of analysis and the elapsed simulation 
time. 



CME without opening the field. As the central arcade continues to rise a current sheet 
forms behind it (Figure 2.10 d). A disconnected flux rope is created due to reconnection 
in this current sheet which will result in flare activity. The feedback between the outward 
expansion drives faster breakout reconnection which, in turn, causes more expansion and 
leads to an explosive eruption. After the secondary reconnection cutsoff the current sheet, 
the side arcades will move inward and a third reconnection phase begins to restore and 
reform the magnetic fields (flare ribbons). There are no analytical expressions for the 
kinematics of the BO model though a number of simulations have produced kinematic 



profiles. Lynch et al. (2004) used a 2.5D simulation to show that the kinematics could be 
well represented with a constant acceleration profile: 



h{t) = ho + vo{t) + -aot'^ 



(2.48) 



3D simulations by Lynch et al. ( 2008 ) resulted in more complex kinematics which consisted 



64 



2.2 Coronal Mass Ejection Theory 



of three phases of constant acceleration: 1) extended low acceleration during shearing 
and breakout; 2) followed by a short period of high acceleration due to the secondary 



reconnections; 3) and the final phase of zero acceleration or constant velocity. DeVore & 



Antiochos (2008j) found a similar but not identical profile with three phases of constant 
acceleration: 1) extended low acceleration during shearing; 2) short fast acceleration with 
breakout and secondary reconnections; 3) and a fast deceleration during the restoration 
phase. 



The DeVore Sz Antiochos ( 2008 ) results may not be applicable to CMEs as they led to 



confined events. The BO model reproduces some observed CME and CME related phenom- 
ena ( van der Hoist et al. , 2007[ ) and flares occur during the secondary (flare reconnection) 
and also during the third phase (ribbon flares). It also reproduces the classic three-part 



CME structure (Lynch et al., 2004) as shown in Figure 2.11, It avoids the Aly-Sturrock 



problem (Section 2.2.1 ) by opening only a portion of the field and not the entire field. The 
BO model can also proceed without reconnection using an ideal instability as shown by 



Rachmeler et al. (2009). 



2.2.2 Propagation 



During the propagation phase the dynamics of CMEs are not very well understood and the 
forces causing the residual acceleration which is sometimes observed are not known. Some 
models have been developed in an effort to understand it in the context of a "drag" force 



for example, the 'snow plough' model (Tappin 2006), or the aerodynamic drag model 



(Cargill, 2004), and others model it in terms of a Lorentz force "flux rope" model (Chen 



1996 ). Examining the physical properties of the CME and the environment into which it is 



launched can give some insight into the mechanisms that might be at play. The solar wind 
is a high-/3 plasma which means that magnetic effects such as reconnection and magnetic 
pressure are less of an influence than the hydrodynamics of the plasma. This also means 
that Alfven waves transporting energy can more easily be converted to fast mode waves 
which efficiently dissipate their energy in the corona. However, CMEs are typically low-/3 
structures even at 1 AU, with plasma-/? of 0.1 common. This means treating CMEs as solid 
bodies moving through a fiow is not a bad approximation. CMEs often drive shocks which 
build-up a compressed plasma region ahead of them, called the sheath. Many CMEs also 
contain the structured magnetic field of a flux rope or, magnetic clouds. 
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2.2.2.1 Aerodynamic Drag 

The aerodynamic drag model invokes the drag caused by fluid flow over a fixed body. If 
we define the Lorentz, gravitational, pressure gradient and drag forces as Fl, Fq, Fp and 
Fu respectively, we can rewrite the equation of motion {F = ma) of a CME moving in the 
solar wind as: 

= Fl + Fg + Fp + Fd 

at 

= Fl + Fg + Fp- PeACD{Vi-Ve)\Vi-Ve\ (2.49) 

with = J X B, Fg = —GMm/r"^ and Fp = — Vp, where A is the cross-sectional area 
of the CME, C/) is the drag coefficient and subscript i or e refers to internal CME or 
external solar wind quantities respectively, and = M + My where My is the virtual 
mass My « pe/lr, and r is the CME volume. The gravitational force Fg, and pressure 
gradient Fp far from the Sun are assumed to be negligible, and the Lorentz force F^ has 



been shown to be negligible in full MHD simulations of the CME dynamics Cargill et al. 



(1996). Equation 2.49 can be re-written in the following manner: 



-r = TT- = —lCn(vi — Ve)\vi — Ve\, where 7 = —. — — 1-^- (2.50) 

this is a first-order differential equation in terms of w , and can be numerically integrated 
given the initial conditions. 

In order to arrive at the aerodynamic drag model, a number of implicit assumptions 
have been made. The first of these is that the CME can be treated as a solid body in a 
fluid flow, and that there is no feedback between the CME and solar wind; the second that 
the magnetic nature of the CME and solar wind plays no significant role in the interaction. 
Finally, that the assumption that the plasma of the solar wind can be represented as a 
typical fluid. Simulations of flux-ropes in plasma flows have shown that, in most cases, the 
feed-back between the flux- rope and solar wind is small, and they also show that Cd ~ 1- 
Note that the drag coefficient determined from MHD simulations will include some of the 
magnetic and plasma effects that would be left out of the aerodynamic equation. 

It should be noted that in the above equations that pe, Pi, t, A and Ve are all functions 
of r. If we wish to describe the total CME evolution then we would need to use the equation 
of state describing the internal pressure, temperature, etc. in balance with the external 
parameters, and let the CME evolve accordingly. This would yield the theoretical evolution 
of the parameters mentioned above, but as of yet none of the models have encompassed 
this much detail. However, statistical studies have been performed, and values extracted. 
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by fitting a function of the form < 7 >= ajiR to velocity-time data and best-fit 



parameters were obtained (Vrsnak, 2001). In Figure 2.12 we see a number of simulations 



of the deceleration a CME due to aerodynamic drag, which were obtained by numerical 



integration of Equation 2.51 using a combination of the best-fit parameters from Vrsnak 



(2001) and realistic estimates. Various simulations were carried out for different initial 



distances and velocities Figure 2.12 



dr 



(2.51) 



Reiner et al. ( 1998 ) pointed out that, at least for some fast CMEs, a quadratic form of drag 



was inconsistent with the observations. Specifically they found that the velocity profile de- 
rived from CME-driven shocks producing Type II radio bursts could not be reproduced by 
a quadratic drag model. Instead the authors suggested a linear model was more appropri- 



ate in reproducing the observed kinematics and Vrsnak &: Gopalswamy ( 2002 ) and Vrsnak 



(2006) also investigated a linear model. A generalised drag equation can be written which 



will take both possibilities into account: 



dr 



r 



(2.52) 



where the sign of the force must be taken into account positive if Vsw > Vcme, negative 
otherwise. The a and (3 values will be different depending on the form of the drag which 
is specified by setting 6 = 1 for linear and 6 = 2 for quadratic models. 



2.2.2.2 'Snow Plough' Model 

The 'snow plough' model is based on the conservation of momentum. As the CME sweeps 
up the solar wind ahead of it, theismaterial must be accelerated, thus momentum is trans- 
ferred from the CME to the swept up material. This process gives rise to two coupled 
differential equations: 

dv dM {vi — Ve) 



dt dt M 
dM 



(2.53) 



PeA{vi - Ve) (2.54) 
We can combine these two equations and substitute in the equation for mass to get: 

= - — A{vi - Ve){Vi - Ve) (2.55) 
at rpi 
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Figure 2.12: Numerical integration of Equation 2.51 for various parameters. All simulations 
use a model solar wind with asymptotic speed of 400 kms^^, different distances for the driving 
force cut-off were assumed. Bold lines show deceleration for = 2 x 10"'^ and = 1.5 
for initial velocities vq =1000, 600, 400 and 200 kms""'^ (curves la, 2, 3 and 4 respectively). 
The curves labelled lb and Ic correspond to the deceleration with a n = 10~^, /Sr = 1.5 and 
an = 2 X 10~^, Pr — 1 respectively both with vq =1000 kms~^ (Vrsnak & Gopalswamy 
2002). The a;- axis units are Rsun- 



which is a firstorder differential equation describing the motion of a CME due to a 'snow 



plough' interaction. Figure 2.13 shows a comparison of modelled CME propagation to the 
observations for the 'snow plough' and aerodynamic drag models. 



By looking at (2.55) and (2.50) we can see that while the physical mechanisms behind 
the snow plough and aerodynamic drag models are different, they have very similar math- 



ematical forms. From (2.50) we can see that in the limit when pe Pi we recover the 
simpler snow plough model. 



2.2.2.3 Flux Rope Model 



Chen (1996) has extended his model to include the interplanetary propagation of CMEs: 



close to the sun the forces on the flux rope are described by |2.45 and far from the Sun 



the forces due to gravity, coronal pressure and the solar magnetic field are not significant. 
The equation of motion is reduced to: 



Fr 



It 
c^R 



In 



8R\ 1 Bt 6 
a J 2 Bpa^ 2 



(2.56) 
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Figure 2.13: Modelled height-time profiles for a transient using the 'snow plough' (dashed 
curve) and aerodynamic drag (dash-dot curve) models compared with the LASCO, SMEI and 
Ulysses observations (heavy lines and + signs). The Ulysses shock speeds are indicated by the 



dotted line segments through their locations ( Tappin 2006 ) 



far from the Sun. He assumed that the drag force was an aerodynamic type of the form: 



V V. 



(2.57) 



where pa is the solar wind density and mi is the internal mass of the flux rope. Figure 2.14 



shows the evolution of CME in this model out into interplanetary space. The CME is 
accelerated by the Lorentz force over the first ~24 hours before drag takes over and begins 
to equalise the CME's speed to the background solar wind. For the initial speed used 
in this simulation this corresponds to a distance of about 70Rq or an elongation of 20°. 
While the exact source of long duration acceleration is not known it has been observed 



in a number of CMEs at large distances from the Sun ( Howard et al. , 2007 Manoharan 



2006 Manoharan et al. 2001; Tappin, 2006). 



2.3 Shocks 

A shock wave or simply a shock is a disturbance across which properties of the medium 
such as pressure, density, velocity, temperature change in a nearly discontinuous manner. 
Shock waves may be generated through a number of mechanisms: 1) abrupt changes in the 
properties of the medium, for example caused by explosions; 2) propagation at supercritical 
velocities such as supersonic aircraft; and 3) a non-linear wave steepening, for example 
waves in the ocean. Shocks are formed across many scales and in different conditions, 
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Figure 2.14: Toroidal instability simulation of CME propagation (Chen 1996). (a) Apex 
height, (b) apex velocity and (c) minor radius dynamics a{t) in 10^ km and 'w{t) in 50kms~^. 
Dotted curve simulation with more faster flux injection, (botton) The Lorentz and drag forces 
acting on the CME, dashed curve is net force. The units are in lO^dyn. 



from astrophysical shocks such as planetary bow shocks (Slavin &; Holzer 1981), or the 



shock at the edge of the Hehosphere (Decker et al. 2008 van Buren et al, 1995) to 



shocks generated by the re-entry of the Apollo mission capsules (Glass, 1977). Figure 2.15 



illustrates the bow shock caused by a supersonic projectile in a wind tunnel and from a star 
moving through the ISM. Shocks can be categorised in terms of the angle between the flow 
direction and the shock normal as, parallel (or normal), perpendicular, or oblique, and as 
moving or stationary shocks. Shocks caused by supersonic flows over bodies form attached 
to the body or, if the flow is deflected more than some critical angle a detached shock 
will form. The distance at which the detached shock forms called the stand-off distance is 
a complex function of both the shape of the object and properties of the medium. Bow 
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Figure 2.15: (top) A shadowgraph of a bullet traveling through air at about 1.5 times the 
speed of sound (or Af = 1.5) a bow shock can be seen ahead of the bullet and a turbulent wake 
behind it. (bottom) A Hubble Space Telescope observation of the bow shock around the very 
young star, LL Ori in the Orion nebula. Images courtesy of NASA.GQV| 



shocks are detached shocks which occur when a blunt object moves relative to a medium at 



supersonic (supeaflvenic) speeds ( Rathakrishnan , 2010). The difference between a moving 
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and standing shock can be understood as follows: in a moving shock the shock moves 
through a medium supersonically, while in a standing shock the shock is at rest and the 
medium flows supersonically. The difference between them is only the frame of reference 
so the system is invariant under Galilean transforms. 

In gas-dynamic shocks (coUisional shocks), the important physical processes are the 
collisions between the molecules. They allow temperature and density perturbations to 
propagate, and the associated viscous forces lead to dissipation. Collisions also allow the 
temperature equalisation of different species of molecules. Gas-dynamic shocks are often 
understood in terms of the macroscopic fluid equations containing quantities such as den- 
sity, temperature and bulk (rather than thermal) velocity. Space plasmas, on the other 
hand, are rarefied and as a result collisions are very infrequent. The lack of collisions means 
that different particles can have very different temperatures which can be non-Maxwcllian, 
and in the presence of magnetic fields can even lead to anisotropic temperature distribu- 
tions. The dissipation mechanisms are complex, involving the interaction between fields 
and particles. In these collision-less shocks the ensemble interaction of the particles and the 
fields is the important physical processes. MHD describes these interactions between the 
large scale field and particles in terms of macroscopic density, pressure and bulk velocity, 
similar to the gas-dynamic case. 

In either gas dynamic or MHD cases, these continuum descriptions can not describe 
the shock itself, as they can not represent the kinetic microscopic processes that control 
the shock; particle collisions in the gas-dynamic case and wave-particle interactions in 
plasma shocks. However, the continuum descriptions are applicable on either side of the 
shock. Shocks are often considered to be an infinitely thin interface, this this is a good 
approximation when the shock is physically thin compared to the length scales involved, 
but thick with respect to the mean free path or the Debye and ion gyro-radius, in gas- 
dynamic and plasma shocks respectively. By applying the mass, momentum and energy 
conservation laws across the shock, the jump conditions or Rankine-Hugoniot equations 
can be derived. In the following sections, the notation [X] = — Xa is used to give the 
difference between quantity X upstream and downstream of the shock. 

2.3.1 Gas-dynamic Shocks 

The simplest description of a shock is in the shock rest frame, where gas moving superson- 
ically (faster than information can be transmitted) flows into the shock from upstream. At 
the shock irreversible processes alter the speed, density and temperature of the medium 
and ,as a result the out flow downstream is subsonic. A shock is thus an entropy-increasing. 
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Figure 2.16: (left) Normal incidence shock frame and (right) de Hoffmann- Teller frame 



(Kallenrode 2004). 



or irreversible, wave that causes the transition from supersonic to subsonic flow. The Mach 
number (M = v/vc) is defined as the ratio of the shock speed in the upstream medium to 
the sound speed. The Mach number is always greater than unity upstream, and less than 
unity downstream. In the rest frame of the shock, the application of the conservation of 
mass energy and momentum gives the Rankine-Hugoniot equations: 



[pUn] = 
[pul +p] =0 
[pUnUt] = 



— + — P)n„ 







(2.58) 
(2.59) 
(2.60) 

(2.61) 



where u is the flow speed, and Un {ut) is the normal (tangential) component to the shock and 
the other symbols have their usual meaning. Combining (2.58) and (2.59) gives [ut\ = 0, so 
the tangential component of the flow is continuous and thus we can consider a coordinate 
system moving along the shock with speed ut. In this the normal incidence frame (see 



Figure 2.16 (left)) the flow speeds u and n„ are identical. Making a Galilean transform into 
the laboratory frame, the mass continuity equation (2.58) can bewritten [p{vs — Un)] 
This can be rearranged in terms of the shock speed Vg'- 



0. 



PdUn,d PuUn,u 
Pd- Pu 



(2.62) 



It should be noted that the shock speed alone is not a good indicator of the energetics 
involved, as small density change can result in a high shock speed but the total energy in 
terms of the compression ratio and mass motion may be low. 
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2.3.2 Magnetohydrodynamic Shocks 



A crucial difference between the gas-dynamic and MHD shocks is the magnetic field. This 
must be accounted for in the conservation equations and also in terms of both the flow 
direction and magnetic field direction, as these do not have to be parallel. In the normal 



incidence frame Figure 2.16 (left) the upstream flow is normal to the shock and oblique to 
the magnetic field, downstream the flow is oblique to both the flow and shock normal. A 
transformation can be made under which the flow becomes parallel to the magnetic field. 
This is known as the de Hoffmann- Teller transformation or frame. This reference frame 
moves parallel to the shock at the de Hoffmann- Teller speed {vht x B = — E). Thus, the 
induced field in the shock front vanishes. The Rankine-Hugoniot equations for an MHD 
shock are: 



, , B2\ (B-n)B 
pu{u • n) + p + - — I n 

2/io 



[pu • n] = 

= 



u.n.f^ + - 

V 2 7 



7 ^ I (B-n)(B-u) 

- 1 /^o y /^o 







[B • n] = 
[n X (u X B)] = 0. 



(2.63) 
(2.64) 

(2.65) 

(2.66) 
(2.67) 



The normal component of the magnetic field (Bn) must be continuous from (2.66), and 



from (2.67) the tangential component of the electric field must also be continuous. Thus, 



the Rankine-Hugoniot jump conditions are a set of five equations for the five unknowns p, 
u, p, Bn and Bt. The jump conditions allow the calculation of the downstream parameters 
from knowledge of the upstream conditions, or vice versa. 

The solutions to the jump conditions actually describe a number of discontinuities 
which are not necessarily shocks. A contact discontinuity is formed when there is no 
flow across the discontinuity {un = 0) and is associated with a density jump while all 
other parameters remain unchanged. The magnetic field has a component normal to the 
discontinuity {B^ 7^ 0), so the two sides are not completely separate, but tied to move at 
the same tangential speed ut. A tangential discontinuity completely separates two regions 
no flux crosses the boundary {Bn = and n„ = 0), and the tangential components change 
{[Bt] 7^ and [ut] 7^ 0). The plasma and field properties can change arbitrarily across the 
boundary, but static pressure balance is maintained, i.e. [p + B"^ /2pq\ = 0. A rotational 
discontinuity requires pressure equlibrium but flux flow across the boundary ti„ 7^ and 
Bn 7^ 0. The normal flow speed is u„ = Bn/ yJppo and the change in the tangential flow 
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speed is related to the change in the tangential magnetic field [ut] = [B^ / y/pno] ■ These 
rotational discontinuities can be viewed as large amplitude waves and are related to the 
transport of magnetic signals through their dependance on the Alfven speed. 

A shock differs from the discontinuities above as there is a flow across the boundary 
{un 7^ 0) with compression and changes to the flow speed. An important parameter for 
MDH shocks is the angle between the magnetic field and shock normal Obv, and this leads to 
classification as perpendicular {Obti = 90°), parallel {Obu = 0°) or oblique for intermediate 
values. Oblique shocks are often subdivided into quasi-parallel (0° < Obti < 45°) and quasi- 
perpendicular (45° < Obtl < 90°) shocks. In a parallel shock Bt = and the magnetic 
field is unchanged by the shock and behaves like a gas-dynamic shock except the collective 
interactions are mediated by the field, and not through collisions. In a perpendicular shock 
Bn = 0, and both the pressure and magnetic field strength change. For the MHD jump 
conditions we can write the shock speed as: 

Pd^d - Pu^u 

Vs = n 

Pd - Pu 

In a gas-dynamic shock there is only one critical speed called the sound speed, Vc, 
while in the MHD case different wave modes have different critical speeds resulting in 
three critical speeds - the slow, fast, and Alfven speeds. Only the slow and fast mode are 
compressive and form true shocks; the Alfven mode only forms shocks in an anisotropic 
plasma, while in an isotropic plasma it results in a rotational discontinuity. The phase 
speed of the slow and fast modes is given by: 

^vha = i^c + 4) ± ^{v^c+v\?-4.v^cvlcos'^0 (2.69) 

taking the positive value gives the fast mode, and the negative gives the slow mode. Taking 
6 = 90° gives the fast magnctosonic speed as v„is = \J'^c ~^ "^a^ the other hand taking 
9 = 0° gives two solutions va and Vg . As in the gas-dynamic case we can define critical ratios 
in terms of the characteristic speed, the fast magnetosonic Mach number {M^g = v/vms) 
and the Alfven Mach number {Ma = v/va)- 



(2.68) 
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Chapter 3 

CME Observations and 
Instrumentation 



In this chapter the details of CME observations, instrumentation, the data reduction pro- 
cess, and conversion of data to physical coordinates are described. This begins with an 
introduction to the mechanism which allows CMEs to he imaged namely Thomson scatter- 
ing and its effect on the observations. Next, the design principles behind the first corona- 
graph are outlined. Following this detailed descriptions of the relevant instruments from the 
Solar and Heliospheric Observatory (SOHO) and Solar Terrestrial Relation Observatory 
(STEREO) missions are given. Finally, a discussion on the conversion of data coordinates 
to physical coordinates and the various coordinate systems is presented. 
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3.1 Observations of Coronal Mass Ejections 



The majority of the radiation from a CME is due to the scattering of photospheric radia- 
tion by free electrons in the CME plasma. While some CMEs contain bright prominence 
material which primarily emits at Ha (6563 A), most CMEs are dominated by the scattered 
radiation. The scattering mechanism is Thomson scattering which is a special case of the 
general theory of the scattering of electromagnetic waves (Jackson, 1975). Thomson scat- 
tering applies when: (1) the coherence length of the radiation is small compared with the 
separation of the particles or, in the case of an incoherent source such as the photosphere, 
the wavelength must be small compared to the typical separation of the particles; and (2) 
the rest mass energy of the scatterers greatly exceeds the photon energy. This is the case 
throughout the corona for white- light. In the case of radio waves condition (1) is not met 
and coherent theory must be used. 



3.1.1 Thomson Scattering 

An unpolarised monochromatic plane electromagnetic (EM) wave incident on an electron 
will accelerate the electron which will then radiate symmetrically about the direction of the 
incident wave. Since the electric field of an EM wave is always perpendicular the direction 
of propagation, the acceleration of the electron will be confined to the plane perpendicular 
to the propagation direction. An observer at a scattering angle (x) of 0° or 180° would 
see unpolarised light, while one at x = 90° would see only linearly polarised light, though 



both observe the same intensity (see Figure 3.1 ). The differential cross-section for Thomson 
scattering by an electron is: 

where a is the cross-section, du is solid angle element at the scattering angle Xi ^ is the 
electron charge and me is the mass of the electron. Integrating over all solid angles gives 
the total cross-section for scattering as: 



8^ / y _ 85 

3 {AireomeC^ J 3 



^t=^(7-— ^) =^^e [m'] (3.2) 
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Figure 3.1: Schematic demonstrating how the angular variation in Thomson scattering arises 

(a) the conceptual set up, the scattering angle x is included for the oblique observer (O2), 

(b) (d) the scattered electric vectors as seen by observers at Oi at % = 180°, O2 at x = 60° 
andOs at x = 90° respectively (Howard & Tappin 2009). 



where rg is the classical electron radius. The differential cross-section for perpendicular 
scattering, cTg, can be written as: 

Thus far, we have only considered the scattering of light from a point source by a single 
electron. The solar photosphere is neither a point source nor uniform, and so it is necessary 
to integrate the scattering component over the light from the whole visible disk. The 
intensity of radiation from the solar photosphere decreases towards the limb, this is known 
as limb darkening and can be characterised by: 

I = -u + ucoscf)) (3.4) 
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where u is the Hmb-darkening coefficient (a function of wavelength), and (p is the angle 
between the emitted radiation and the radius vector of the emitting point. The tangential 
(It) and radial (//j) vector components of Thomson-scattered radiation from a single 
electron in the solar corona may be calculated in terms of a small number of measurable 



parameters (Bilhngs, 1966 Howard & Tappin 2009; Minnaert, 1930; van de Hulst 1950) 



as: 



and 



where 



It = lo^ [{I - u)C + uD] 



Ip = sin^ X [(1 - u)A + uB] 



A = cos ft sin 
1 



B 



cos^$7,, „ . 9^,, /l + sinri 



l-3sin2jl — (l + 3sin2l^)ln , 

sin U V cos U 



4 cos^^l 

C = cos S I , 

3 3 

cos^ O 



D 



1 



5 + 3 sin^ Q 



sinO 



(5 - sin^ n) In 



1 + sin 
cos Q 



(3.5) 

(3.6) 

(3.7) 
(3.8) 

(3.9) 
(3.10) 



and is the solid angle of the Sun as seen by the scatterer and Ip = It — Ir- The total 
scattered light intensity is therefore: 



hot = {It + Ir) = 2/t - Ip- 



(3.11) 



The locus of all the points that maximise the scattered intensity form what is known 



the Thomson sphere (TS) as shown in Figure [3^ (left ) . Figure 3.2 (right) shows the ratio 



of the correct to assumed plane-of-sky brightness and indicates it is a valid assumption 
out to about 70 Rqir, after which the brightest feature may be far from the plane-of-sky. 
There is an often quoted misconception that scattering efficiency is maximised on the TS 
which it is not in fact the opposite occurs as the cross-section for scattering is minimised 
on the TS. However the TS also indicates the point of closest approach, thus the amount 
of incident radiation is maximised as is the density of scatterers, and hence the scattered 
intensity is maximised for the above constraint. 

While the effects of Thomson scattering on the plane-of-sky assumptions may be neg- 
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Figure 3.2: Thomson sphere geometry (left) and the range of vaUdity of the plane-of-the- 
sky assumption (right). The ratio Bumb/Bo is the ratio of the brightness calculated using the 



plane-of-the-sky assumption over the brightness derived from the full treatment (Vourlidas & 



Howard 



2006). 



ligible close to the Sun they severely impede the ability to unambiguously identify CME 



structures in observations, as Figure |3 . 3| shows . The observations (Figure 3.3 left) the show 
the graduated cylindrical shell (GCS) model (Thernisien et al. , [2006 ) almost edge-on and 



nearly side-on, but in the synthetic data (Figure 3.3 right) or observations it is nearly im- 



possible to disguising between the two orientations and certainly not unambiguously. This 
highlights the problems in interpreting coronagraph observations and trying to compare 
them to theory. 



3.1.2 Projection Effects 

Close to the Sun, the plane-of-sky assumption will only affect the magnitude of the de- 
rived kinematics, but not their profile shape. For example, any acceleration seen is real, 
but its magnitude is subject to possibly large uncertainties. Using on-disk features and 
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assumptions about the propagation direction can allow accurate kinematics to be derived 
for selected events, especially those on the limb. Even then, care must be taken in their 



interpretation (Byrne et al. , 2009). However, far from the Sun the Thomson scattering and 



its effects on the plane-of-sky assumption become much more important. Figure 3.4 shows 
both bubble and shell models of a CME at different times in its propagation. The differ- 
ences between the observed (x), true (Q) and inferred (+) fronts are clear, and become 
extremely large as the CME reaches observer- like distances. The effects of making the 



plane-of-sky assumption when deriving kinematics are shown in Figure 3.5 which shows 
simulations of the true kinematics (black), as well as those that would be observed from 
STEREO A (red) and B (green) like positions for a point-like CME propagating 12° east of 
the Sun-Earth line. The CME is close to the plane-of-sky for STEREO B but far from A as 



seen in the inset (Figure 3.5 d). Not only are the kinematics significantly under estimated 



close to the Sun from and B, as the CME distance increases apparent acceleration, due to 
the plane-of-sky assumption begins to take place. It should be clear that the plane-of-sky 
assumption is unsuitable for studying CMEs at large distances from the Sun. Other ap- 
proaches have been developed such as the "point-P" method or the "fixed-i;^>" method. The 



"point-P" method (Howard et al., 2006) assumes the CME is intrinsically a very broad. 



uniform, spherical front, centred on the Sun, with the relation ship between the observed 
elongation to height given by: 



dsine 



(3.12) 
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Figure 3.4: Basic CME structures, (a) The expanding bubble and (b) a simple shell with 
an angle of 30 at three different times during their propagation. The tangent drawn from the 
observer O across the CME surface shows the location of the relative leading edge. The x 
symbols represent the location of the leading edge seen by the observer the Q symbols show 
the true location of the leading edge at the central location and the + symbols the inferred 



location of the leading edge based on the central location ( Howard & Tappin 2009 1 



where d is the Sun-observer distance, and e is the elongation angle. The "fixed-i?ii" method 



(Kahler & Webb, 2007) assumes that the CME is a relatively narrow, compact structure 



travelling on a fixed, radial trajectory at an angle, (j) relative to the observer's line-of-sight 
to the Sun and the derived height is given by: 



dsin e 
sin e + 



(3.13) 



More recently another method known as the "harmonic mean" has been developed ( Lugaz 



et al. , 2009 ) . This is so called as it corresponds to the harmonic mean of values calculated 



from the other two methods the derived height is given by: 



2d sin e 



1 + sin e + 



(3.14) 



3.2 Coronagraphs 



A coronagraph is a telescope which aims to reproduce the spectacular images of the solar 
corona during eclipse observations. The two main problems with developing a coronagraph 
are: the scattered light in the telescope, and the sky brightness itself, both of which will 
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Figure 3.5: Plane-of-Sky projected kinematics representative of STEREO A (Red), STEREO 
B (green) and true (black) for a CME propagating 12° east of the Sun-Earth line, (a) Height, 
(b) velocity (c) acceleration and (d) geometry of simulation. 



overwhelm the weak corona. Bernard Lyot solved these problems and constructed the first 



coronagraph in 1930 (Lyot 1939). The optical layout of this coronagraph (see Figure 3.6) 
is designed to reduce internal scattering and blockofF the solar disk. The first element was 
an objective lens that was free from inclusions, and highly polished to reduce scattering 
and reflections. Directly behind the objective lens (01), and in its focal plane, is the 
artificial moon or occulter which reflects away the bright solar disk image. The first field 
lens (Fl) forms an image of the objective lens, and a screen placed here would show an 
image of the objective with a bright halo and central bright spot due to refraction, and 
double reflection, respectively. Lyot introduced what are now known as the Lyot spot and 
Lyot stop to remove these aberrations. Finally the second objective lens (02) forms an 
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Figure 3.6: Optical layout of Lyot internally occulted coronagraph adapted from Lyot ( 1939 ) 



image of the occulter (blocking the solar disk) and corona on to the detector plane. 

There are two different configurations for coronagraphs, internally occulted (occulter 



is inside the first optical element, 01 in Figure 3.6), and externally occulted (occulter 



outside of the first optical element, 01 in Figure 3.6). Each configuration has its own 



unique advantages and disadvantages. The internally occulted coronagraph allows the 
very inner corona to be imaged, but as the front lens and aperture are directly illuminated 
by the solar disk, scattered light within the telescope is a problem. The externally occulted 
coronagraph cannot image the inner corona due to the diffraction limit (in space, the sizes 
of the telescopes are limited) , but scattered light is not as much of a problem as the aperture 
and objective lens are not directly illuminated by the solar disk. Modern coronagraphs 
use the same basic design, but have added baffles to reduce internal stray light further and 
polarisers to separate the signature of CMEs from that of the corona. The coronal light is 
a combination of polarised and unpolarised light, while CMEs contain only polarised light. 
Polarisers block some of the coronal light, allowing the CME to be more easily identified. 

Coronagraphs are limited to observe only a small distance (<3O-R0) from the Sun 
and cannot image CMEs as they propagate to the Earth at ~215i?0, though this isone 
of the 'holy grails' of solar physics. Extended coronagraphs or heliospheric imagers aim 
to fulfil this goal and two instruments have already validated the ability to measure the 
electron scattered CME signal against the strong zodiacal light and stellar background. The 



Zodiacal Light Photometer (Pitz et al. 1976) on the Helios spacecraft launched in 1974, 
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and the Solar Mass Ejection Imager (SMEI; Eyles et a/.||2003" ) instrument, on the Coriohs 
spacecraft launched in 2003, have demonstrated the ability to detect and track CMEs far 
from the Sun (Tappin et al., 2004[ ). The Heliospheric Imagers on-board STEREO image 
from ~4-88° elongation (15-215 and thus, for the first time, allow CMEs to be 

tracked from the Sun to the Earth from two perspectives. The results in this thesis are 
largely based on STEREO SECCHI observations in particular simultaneous (A and B) 
observations of CMEs in CORl, C0R2, and HI. A summary of the instruments and their 



important features and properties is listed in Table 3.1 



3.3 Solar and Heliospheric Observatory (SOHO) 



SOHO (Domingo et al., 1995) was a joint ESA NASA mission which was launched on 
December 2, 1995 by an Atlas II rocket. SOHO was launched into an orbit around the 
first Lagrange point (LI), located along the Sun-Earth line, about 1% of the distance to 
the Sun, thus allowing continual monitoring of the Sun. The main scientific goals of the 
mission were to: (i) study the solar interior using helioseismology techniques; (ii) study the 
heating of the solar corona; and, (iii) investigate the solar wind and its acceleration. There 
are twelve complementary instruments onboard: 3 devoted to helioseismology probing 
the inner structure of the Sun by observing solar oscillations; 3 measuring in-situ plasma 
properties (densities, velocities, magnetic field, composition, etc.) of the solar wind; and 
6 telescopes with imagers or spectrometers studying the solar disk and atmosphere. A 
schematic of the SOHO spacecraft and its instruments is shown in Figure [377| 



Both Extreme-Ultraviolet Imaging Telescope (EIT) and Large Angle Spectrometric 
Coronagraph (LASCO) use the same CCD detectors, which are STIe (formerly Tektronics) 
1024x1024, three-phase, multi-phased-pinned type with square pixels of 21 fim. The full- 
well capacity of the CCDs is about 150,000 electrons and the quantum efficiency is about 
0.3-0.5 in the optical range 500- 700 nm when using front side illumination. For EIT, the 
CCDs were thinned, back-side illuminated, and had anti-refiective coatings applied in order 
to keep the quantum efficiency high at 0.27-0.36 in the EUV 171.4-303.8A range. The 
CCDs have 40 non-imaging (over- or under-scan) pixels on each line which can be used for 
calibration and engineering purposes. 
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Figure 3.7: A schematic of tlie SOHO spacecraft and instruments (Domingo et al. 1995) 
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Figure 3.8: Scliematic of ttie EIT telescope ( Delaboudiniere et al.. 19951 



3.3.1 Extreme-Ultraviolet Imaging Telescope (EIT) 



EIT (Delaboudiniere et al.. 1995) is a Ritchey-Chretien telescope which images the corona 



and transition region on the solar disk, and up to 1.5 Rq above the solar limb. Light 



'Depends strongly on the distance of the object from the plane-of-sky 
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2002/02/27 01:00 



Figure 3.9: Observation showing active regions and magnetic loops as recorded by EIT in the 
Fe IX/X 171 hne. The temperature of this material is about 1 million K in the lower corona. 
Image courtesy of the |SOHO website[ 



enters the telescope through an initial filter of 700 A of cellulose sandwiched between 
two 1500 A thick aluminium films, which reject visible and IR radiation, and then passes 
through a quadrant selector (see Figure 3.8). Each quadrant of the primary and secondary 
mirrors has multilayer coatings of molybdenum and silicon deposited on them, each was 
optimised to reflect in a different wavelength bands centred on 304 A, 284 A, 195 A, and 
171 A. Probing the solar atmosphere at peak temperatures of 8.0 x 10^ K, 2.0 x 10^ K, 
1.6 X 10^ K, and 1.3 x 10^ K respectively (assuming quiet sun conditions). Finally, a set of 
filters just in front of the focal plane block long wavelength radiation before imaging, and 
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a shutter controls the exposure time. An image is formed on the CCD and as EIT has a 
field-of-view 45arcmin square, each pixel is 2.6arcsec in size. A typical image from EIT 
from the 171 A channel dominated by Fe ix and x lines is shown in Figure 3.9 



3.3.2 Large Angle Spectrometric Coronagraph (LASCO) 



The LASCO (Brueckner et al., 1995) instrument consist of three coronagraphs which im- 
age the corona from l.l-SOR© (CI: l.l-Si?©, C2: 1.5-6.0 C3: 3.7-30 i?©). The 
CI coronagraph has not operated since the loss and subsequent recovery of the spacecraft 
in 1998. C2 is an externally occulted Lyot type coronagraph. Light enters the C2 coron- 
agraph and encounters the external occulter, which was a new design a cone with sharp 
threads whose cone angle is slightly larger than the angle subtended by the Sun from LI 
(Figure 3.10). This design achieves a light rejection level of 1.5xlO~^, and at this level, 



light refracted by the next element in the path (the entrance aperture) becomes a concern. 
As a result, a polygonal serrated aperture was developed, with each side behaving as a knife 



edge, such that its direct diffraction avoids 01 (Figure 3.10). Light then passes through 



90 



3.3 Solar and Heliospheric Observatory (SOHO) 



an oversized diaphragm A2, preventing propagation of light scattered off the edges of the 
aperture, and enters the objective 02. 02 is a two element design which apodises light 
diffracted by 01, and ghost images created by 01 are blocked by the Lyot spot (a metallic 
layer deposited onto 03). The final element, 03, is a four element design which magnifies 
and focuses an image of the corona seen by 02 onto the CCD detector, and results in an 
angular resolution of 11.4arcmin per pixel. Two plane mirrors fold the optical path to 
reduce the over all instrument length. The shutter and filter wheel are mounted as close 
to the final pupil as possible, and behind 03, respectively. The polariser is placed closer 
to the image plane, in front of the CCD detector. The pointing error system detects any 
imbalance in the penumbra created by 01 using four photodiodes behind four symmetric 
holes located about 01. The differential output of paired diodes is amplified, coded, and 
telemetered to the ground where commands can be sent to move the pointing legs until 
satisfactory pointing is achieved. 

C3 is also an externally occulted Lyot type coronagraph. Light enters the instrument 
through AO, where the occulting disk shadows the coronagraph entrance aperture, Al, from 



direct sunlight (Figure 3.11 ). The occulting disk is an assembly of three disk on a common 
spindle configured to minimise diffracted sunlight falling on Al and the primary lens. The 
primary lens forms an image of the occulted corona and is followed by the internal occulter 
which blocks the image of the external occulter and halo. Behind this is the field lens, a 
multi-element anti-refiection coated lens which forms a collimated image of the corona on 
the relay lens. It also presents a sharp image of the instrument aperture Al to the Lyot 
stop. The elements of the relay lens are all placed behind this stop to minimise diffracted 
light falling on it. The relay lens focuses an image of the corona on to the CCD detector, 
and it also contains the Lyot spot. The resulting pixel resolution is 56arcmin per pixel. 
Between the relay lens and CCD lies a triple mechanism containing the filter, polariser, 
and shutter. Centred baffles are placed along the instrument so that the field, objective 
and relay lenses see only the rear of the baffles, or the walls shadowed by the preceding 
baffle. All the surfaces are anodised and coated in a dense black paint. 

Both C2 and C3 have a number of filters in the range 400 - 850 nm for C2, 400 - 1050 nm 
for C3 and polarisers at 0°, ±60°. The polarisers are used to obtain total brightness B, 
or polarised brightness pB, images via a combination of observations taken at polariser 



positions la = —60°, = 0°, /c = +60°, combined using (Billings, 1966): 



B='^{Ia + h + Ic) (3.15) 



^^=3L 



{la + h + hf - 3 {lah + IJc + hh) • (3.16) 
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Figure 3.11: Optical layout of the C3 coronagraph (Brueckner et al. 19951. 



A sample LASCO observation from C3 is shown in Figure 3.12 



3.4 Solar Terrestrial Relation Observatory (STEREO) 



The NASA STEREO (Kaiser et al. 2008) mission was lunched on a Delta 2 rocket on 



October 26, 2006. The mission consists of two near-identical spacecraft which after a 
series highly eccentric Earth orbits and Moon sling shots, were placed in orbit around 
the Sun. The Ahead spacecraft (STEREO-A) travels ahead of the Earth, slightly closer 
to the Sun and the Behind spacecraft (STEREO-B) lags behind Earth, slightly further 
from the Sun. The spacecraft separate at about ~ 22° per year from Earth, or ~ 45° as 
viewed from the Sun. Figure |3.14 shows the positions of the spacecraft at three different 
times. The STEREO mission was designed to study the causes and mechanisms behind 



92 



3.4 Solar Terrestrial Relation Observatory (STEREO) 




Figure 3.12: A typical C3 observation, the occulter, support pillar and non-imaging pixel 
have been set to zero. The white circle on the occulter indicates the Sun's size. A CME is 



visible off to the right of the occulter. Image courtesy of the SOHO website 



the initiation of CMEs, and then follow their propagation through the inner Heliosphere. 
STEREO will also be used to study the site of energetic acceleration and develop 3D 
models of solar wind properties such as magnetic topology, temperature, velocity, and 
density. As such, the STEREO spacecraft carry an almost identical payload consisting of 
optical, radio, as well as in situ particle and field instruments. These are divided into four 



suites: Sun-Earth Connection Coronal and Heliospheric Investigation (SECCHI; Howard 
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Figure 3.13: Schematic of tlie STEREO B spacecraft . Tiie positions of the some of the 



instruments and their detectors are indicated ( Kaiser et al. 2008 ) 



et Q/.|2008a ); In situ Measurements of PArticles and CME Transients (IMPACT; Luhmann 



et a?.|2008D ; PLAsma and SupraThermal Ion Composition (PLASTIC; [Calvin et a/.|2008[ ); 



and STEREO/WAVES (S/WAVES; [Bougeret et~al\^(m^ . Figure |3l3l shows one of the 
STEREO spacecraft and its instruments. 



3.4.1 Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) 

The remote sensing optical suite, SECCHI, consist of an extreme ultraviolet imager (EUVI) 
which images the solar disk in four wavelengths 171 A, 195 A, 284 A and 304 A, two coro- 
nagraphs (CORl and C0R2) which image the corona from 1.4-15i?0 in white-light, and 
the Heliospheric Imagers (HI) which image the inner Heliosphere (~4-88° elongation) in 
white-light. This combination of instruments can image from the solar surface to beyond 
1 AU from two perspectives. Combining the two perspectives allows for the possibility of 
3D reconstructions of solar features such as, loops, prominences and CMEs. 

Each of the scientific instruments on SECCHI uses a three-phase, back-illuminated, 
non-inverted mode (to ensure good full well capacity, 150k to 200k electrons) CCD model 
CCD42-40 manufactured by E2V in the United Kingdom. There are 2048x2052 image 
pixels, each measuring 13.5 ^m on a side and providing a total imaging area of 27.6 mm 
square. The total readout is 2148x2052 pixels, providing 100 columns of non-imaging 
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Figure 3.15: A cross section of the EUVI telescope (Wuelser et at. 2004) 



over- and under-scan regions for calibration and engineering purposes. The CCDs used 
for visible light detection (COR-1, COR-2, HI-1 and HI-2) have an anti-reflective coating 
on their backside (illuminated side). The quantum efficiency (QE) of these devices varies 
from roughly 80% at 500 nm, to 34% at 900 nm. The backside of the EUVI CCD has no 
coating in order to provide sensitivity at shorter wavelengths, with a quantum efficiency 
of 74% at 17.1 nm, and 70% at 30.3 nm. 



3.4.1.1 Extreme Ultraviolet Imager (EUVI) 



EUVI (Wuelser et al, 2004) is of similar design to that of SOHO/EIT except it provides 
higher resolution (40% better) and cadence capabilities. It is a Ritchey-Chretien telescope 
which images the solar disk and atmosphere up to 1.7 Rq. It provides pixel limited res- 
olution of l.Garcsec per pixel across the entire field-of-view. Light enters the telescope 
through the entrance filter (150 nm thick aluminium) which blocks undesired UV, optical 
and IR radiation, and passes through the aperture selector to one of the four quadrants. 
Each quadrant of the primary and secondary mirrors uses multilayer reflective coatings, 
optimised for four wavelength bands centred on 171 A, 195 A, 284 A and 304 A. The light 
continues through a filter wheel with redundant thin aluminium filters to remove the re- 
mainder of visible and IR radiation. The image is then formed on the CCD, and a rotating 



blade shutter controls exposure times. Figure 3.15 shows the layout of EUVI. A typical 



EUVI observation in the 195 A is shown in Figure 3.16 



3.4.1.2 CORl and COR2 

The inner SECCHI coronagraph, CORl, is a Lyot internally occulted refractive telescope 
which images the corona from 1.4 to 4i?0 at 656 nm (Ha) at 7.5arcsec per pixel (with 
onboard 2x2 binning). It is the first space-borne internally occulted refractive telescope 
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Figure 3.16: Image of the Sun, taken by the SECCHI Extreme Ultraviolet Imager (EUVI) 
on the STEREO Ahead observatory on March 21, 2011 at 02:40:30 UT. The 195 Angstrom 
bandpass is sensitive to the Fe xii ionisation state of iron, at a characteristic temperature of 
about 1.4 million degrees Kelvin. Image courtesy of the Stereo Science Centre. 



(in contrast to the internally reflective design of LASCO/Cl) and this design enables 
better spatial resolution closer to the limb than an external design. CORl is dominated 
by instrumentally scattered light which cannot be removed by Lyot principles, but as it 
is largely unpolarised it can be removed by talking polarised observations and calculating 
pB images. 



Figure 3.17 shows the opto- mechanical layout of the instrument. Light enters through 
the front aperture where the singlet objective lens focuses the solar image onto the occulter. 
The occulter is mounted on a stem at the centre of the field lens, and the tip of the occulter 
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Figure 3.17: A cut away view of the CORl coronagraph (Howard et al. 2008al 



is cone-shaped to direct sunlight into the surrounding hght trap. The hght trap is wedge- 
shaped, and as a result the light must reflect many times before it can escape and is thus 
effectively absorbed. Two doublet lenses act as a telephotolens and focus the coronal image 



onto the CCD while maintaining diffraction-limited resolution (Figure 3.17). A bandpass 



filter 22.5 nm wide centred on the Ha line at 656 nm is placed just behind the first doublet. 
The Lyot stop, spot, first doublet, and bandpass filter all form a single linear optical 



assembly (Figure 3.17). A linear polariser mounted on a hollow core motor rotational 
stage is located between the two doublets and a rotating blade shutter is located just in 



front of the focal plane detector assembly (Figure 3.17). 



The objective lens only focuses the solar image accurately at one wavelength (Ha), 
due to chromatic aberration, so the internal occulter is designed to block all other the 
solar photospheric light between 350-llOOnm. Diffracted light from the front aperture is 
focused onto a Lyot spot by the field lens, this removes the largest source of stray light. 
Placing baffles at various points between the front aperture and the Lyot stop removes 
additional stray light. To remove ghosting from the objective lens a Lyot spot is glued to 



the front surface of the doublet lens immediately behind the Lyot stop (Figure 3.17). A 



focal plane mask located between the shutter and detector removes light diffracted from 
the edge of the occulter. At orbital perigee at the design wavelength the resulting solar 
image is completely occulted to IARq and vignetted to 1.9 i?0. 
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Figure 3.18: A cut away view of the C0R2 coronagraph (Howard et al. 2008a|. 



The outer coronagraph C0R2 is an externally occulted Lyot coronagraph, of similar 
design to that of the SOHO/LASCO C2 and C3 coronagraphs, which images the corona in 



white light (650- 750 nm) from 2.5- 15 with 14.7 arcsec per pixel (2048x2048; Howard 



et al. 2008a). Solar radiation enters the instrument through the aperture AO (see Fig- 



ure 



3.18) where a three-disk occulter shades the objective lens from direct solar radiation 



and creates a deep shadow at the lens aperture Al. A mirror rejects heat back through the 
entrance aperture. The objective lens (01) focuses an image of the external occulter onto 
the internal occulter, and the field lens (02) focuses the Al aperture onto the A3 aperture 
or Lyot stop. The objective lens also creates an image of AO onto the field stop A2. The 
internal occulter and aperture stops block images of the brightly illuminated edges of the 
external occulter, AO, and Al apertures. Between the 02 and 03 elements lies a linear 
polariser mounted in a hollow-core motor. The third lens group forms an image of the 
corona onto the focal plane where the detector is located. A bandpass filter transmits 
radiation from 650- 750 nm (FWHM) optimising throughput, and a rotating blade shutter 



controls exposure time. A typical C0R2 observation is shown in Figure 3.19 



Both CORl and C0R2 take sequences of three polarised observations at 0° and ±60° for 
a complete polarisation sequence. This results in a cadence of 8 and 15 minutes respectively. 
These images are telemetered to the ground where they can be combined to give both total 
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Figure 3.19: Image of the solar corona, taken by the SECCHI outer coronagraph (C0R2) 
on the STEREO Behind observatory on February 15, 2011 at 03:54:33 UT. Image courtesy of 
the (Stereo Science Centrei 



brightness (B) and polarised brightness (pB) images. 
3.4.1.3 Heliospheric Imager (HI) 

The HI instrument is a combination of two refractive optical telescopes (HIl and HI2) 
with a multi-vane, multi-stage light rejection system which images the inner Heliosphere 
in white-light (HIl: 630-730 nm, HI2: 400-1000 nm ). HIl has a FOV of 20° degrees 
centred on an elongation of 13.28° degrees with 70arcsec pixel"^ plate scale, while HI2 
has a FOV of 70° degrees centred on an elongation of 53.36° degrees from the Sun with 2 
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arcmin pixel ^ plate scale ( Eyles et al. , 2009 ) . 



The HI concept is based on the laboratory measurements of Bufiington et al. (1996), 
who determined the scattering rejection as functions of the number of occulters and their 
angle below the occulting edge. The baffle system consists of the forward, internal, and 
perimeter baffles (see Figure 3.20[ a)). The forward baffle system uses Presnel diffraction 
with the positions and heights of the vanes arranged in an arc such that the (n + 1)*'* 
vane lies below the shadow line of the (n — 1)*'* and n*'* vanes (see Figure 3.21). The 
forward baffle is the primary source for light rejection in the instrument and provides the 
low levels of residual light needed of 10^^*^-10^^^, as computed using Fresnel's second order 
approximation of the Fresnel-Kirchhoff diffraction integral for a semi-infinite half-screen 



see (Figure 3.22 Right). These high levels of light rejection are vital to detect faint CMEs 
(see Figure [O2] Left) 

The internal baffle consists of five principal vanes and a number of secondary vanes 
(see Figure 3.20| [^a)). The system is designed to attenuate unwanted light entering the HI 
optical systems through multiple reflections in the baffles. Finally the perimeter baffles 
(side and rear baffles) are designed to shield the HI optical systems from reflections of 
photospheric light below the horizon defined by the baffles. 



The optical configurations of the two HI telescopes is shown in Figure 3.23 Neither 
of the telescopes is diffraction limited as such they are optimised to minimise the rms 
(root mean square) spot diameter over the operational temperature range. The design 
of the HIl optics was mainly constrained by the location of the first lens which must 
lie in the appropriate position in the shadow pattern of the forward baffle and the lens 
barrel and housing must be below the shadow line cast by the first baffle, so direct sunlight 
cannot reach these components. The spectral response is controlled by means of multilayer 
coatings on two internal lenes which have the minimum angles of incidence. All the lenses 
are coated with anti-reflective (AR) coatings to minimise ghost images. The HI2 optics 
were designed to give a wide field-of-view, wide bandpass and minimise ghosting. Due to 
the wide bandpass the effectiveness of the AR coatings is minimised. To combat this effect 
the optics were designed to spread out ghost images over large regions of the focal plane. 
So the ghost images from bright objects contribute to the overall background rather than 
generating multiple images, with their spacing dependant on the position of the object in 
the FOV. The lens barrels were optimised using ray-tracing software, with some cavities 
being added to act as light traps. All the internal surfaces and mountings were treated 
with either black copper oxide, or black chromium coatings. Due mainly to mechanical 
constraints neither of the HI telescopes have a shutter mechanism. 

The images from both HIl and HI2 are binned on-board to give 1024 x 1024 pixels at 
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Figure 3.20: (a) Heliosphereic Imager concept, (b) Cross sectional view through the instru- 



ment, showing the fields-of-view of the two telescopes (Eyles et al. 2009) 
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Figure 3.21: Concept of the cascade of knife-edge diffraction system, also shown are the 



positions of the entrance apertures for HIl, HI2 and the form of the rejection function (Eyles 



2009). 
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Figure 3.22: Left) The ficlds-of-view and geometry of the HI telescopes and anticipated 
intensity of a typical CME and coronal background. Right) Theoretical light rejection levels 
for HI in typical conditions as a function of distance below the shadow of first vane. Also 



shown are the approximate positions of HIl and HI2 optics ( Eyles et al. 2009 ) 



stated plate scale. In order to attain the required signal-to-noise ratio needed to image the 
extremely faint CMEs, the His need long exposure time images. Due to limited dynamic 
range, and to avoid cosmic ray swamping, a series of shorter exposures are taken and 
summed on-board: in the case of HIl 30 images of 40 s exposure each are summed (total 
exposure time 1200 s) and the resulting cadence is 40 minutes; for HI2 99 images with a 
50 s exposure time each are summed (total exposure time 4950 s) resulting in an image 
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Figure 3.23: Optical layout of the HIl (a) and HI2 (b) telescopes. The Schott glass type is 
indicated beside each element (Eyles et al. 2009). 



cadence of 2 hours. Neither HIl nor HI2 have shutters and hence all images must be 
corrected for the smearing caused by shutterless readout (i.e., each row sees the scene from 
the pixel below it). Some of the objects in the HI FOV are bright enough to saturate 
the CCD, which was designed to bleed along columns allowing these to be identified and 



corrected later (Eyles et al., 2009). A typical HIl observation is shown in Figure 3.24, and 



Figure 3.25 shows a CME as it progresses through the the HIl field-of-view and into HI2. 



3.5 Data Reduction 



The data were reduced and calibrated using SolarSoft library (Freeland &; Handy, 1998). 



Each of the instruments has its own data reduction software which corrects for various 
effects and produces science-grade observations. The data reduction can consist of the 
following: 
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Figure 3.24: Image of the solar corona, taken by the SECCHI inner Hehospheric Imager (HI- 
1) on the STEREO Behind observatory on February 15, 2011 at 14:09:34 UT. Image courtesy 
of the [Stereo Science Center| 

• On-board processing: numerical operations applied on board the spacecraft to 
keep the data within the valid range of the compression algorithm, which may include 
taking the square root of the image or dividing the image by 2, 3 or 4 multiple times, 
the operations are noted in the image header and corrected for when the data is 
reduced. 

• De-biasing is the removal of the electronic offset applied to each pixel. CCD elec- 
tronics are designed to produce a positive offset value to avoid issues with negative 
numbers. For most science grade CCDs the bias is derived from the overscan pixels 



105 



3. CME OBSERVATIONS AND INSTRUMENTATION 




106 



3.5 Data Reduction 



on the detector chip, recorded in the image header and subtracted off when the data 
is reduced. 



• Flat-Fielding is the correction of the image for non-uniform illumination (including 
vignetting) of the detector, and pixel-to-pixel variations (hot pixels) in the CCD 
response. Once a flat-field is derived all the images are divided by this to flatten the 
CCD response when that data is reduced. 



• Bleeding occurs when the a signal exceeds the full well capacity of a pixel. The 
excess charge "bleeds" into adjacent pixels in the same column, causing columns of 
saturated pixels and in severe cases there can be horizontal bleeding as well. These 
regions can be flagged and removed during data reductions. 



• Distortions the optics of telescopes produce distortions such as pin cushion or barrel 
distortion. The images can be geometrically corrected for these during date reduction. 



• Calibration involves converting the digital numbers (DN) to a physical unit such 
as mean solar brightness (MSB). 



For some instruments, backgrounds which aim to remove the static coronal signature 
can be subtracted during the reduction process (default for CORl) or after, to enhance 
the faint coronal signatures. In the case of the His which have no shutter, additional 
corrections need to be applied to account for pixel smearing during readout. Smearing is 
the effect of the extra exposure time a pixel will experience during readout and clear phase 
as a result of there being no shutter. In the case of HI this means a bright object (such 
as a star) will leave vertical trails above the position during read out. Also during the 
clear phase the pixels are exposed leaving vertical trails in the opposite direction. As the 
readout, exposure, and clear times are known this effect can be correct using 

I = T-^R (3.17) 

where I is the corrected image, T is a matrix with the exposure time along the diagonal, 
the readout time above the diagonal and the clear time below, R is the uncorrected image. 
Finally, the pointing information in the HI headers are updated with pointing information 
derived from the star fixes in the images, as dust impacts on the telecscopes can shift the 



optical alignment ( Brown et al. , 2009 ) . 
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Figure 3.26: Geometry of the azimuthal (zenithal) perspective projections (Left), the point 
of projection at P is one spherical radius from the centre of the sphere; the three important 
special cases (right) from Calabretta & Greisen (2002 1. 



3.6 Coordiante Systems 



The derivation of phyical coordinates from pixel coordinates is a multi step procedure. In 
the first step pixel coordinates are transformed to intermediate coordinates i.e., they are 
converted into the relevant units (meters, degrees, arcsec), but they are not adjusted for the 
reference point of the observations or projection/geometrical effects. These intermediate 
coordinates are then transformed by the various possible projections and rotated into 
celestial coordinates. The conversion from pixel to intermediate coordinates is given by: 

N 

Xi = Si'^ niij {pj - Tj) (3.18) 

i=i 

where pi refers to pixel coordinates, rj is the reference pixel, vriij is a linear transformation 
matrix and si is a scale function, with subscript i referring to pixel axes and subscript j 



referring to coordinate axis (Calabretta & Greisen, 2002). 



Traditionally almost all solar imaging telescopes use the gnomonic or 'TAN' projection. 
In this projection a feature at an angle 9 to the instrument's optical axis or bore-sight is 
projected onto the image plane at a radial distance Re given by: 

Re = Ftan{e) (3.19) 

where F is the focal length of the optical system, assuming symmetry about the optical 
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axis. The position of the feature in the image plane is given by: 

x = Re cos((/)) (3.20) 
y = Resm{(j)) (3.21) 

where 4> is the angular position of the feature about the optical axis. However, for the 
wide-angle optics of HI, this type of projection is not adequate, and it was found that a 
Azimuthal Perspective (AZP) could accurately represent the behaviour of the HI optics. 
In this projection R is given by: 

(^ + l)sin(a) 
/X + cos a 



where Fp is the paraxial focal length, and is the distortion parameter. Inspecting (3.22) 
it should be clear that the gnomonic projection is a special case of the AZP projection 
with /i = 0. Figure 3.26| shows the AZP projection for a unit sphere and three special 



cases. This figure shows why, for most solar telescopes which have a narrow field-of-view 
a gnomonic projection is suitable. Once the coordinates have been de-pojected they can 
be converted to the Sun-centred physical coordinates. 



3.6.1 Heliographic Coordinates 

The well-known solar coordinate system expresses the position of a feature on the solar 
surface in terms of latitude, 0, and longitude, and can be extended to 3D by adding 
the radial distance, r, from the centre of the Sun. There are two basic variations of the 
heliographic coordinate system called Stonyhurst and Carrington which both use the same 
rotation axis only differing in longitude. 



3.6.1.1 Stonyhurst 

The origin of the Stonyhurst coordinate system is at the intersection of the solar equator 
and central meridian as seen from Earth. B and ^> are given in degrees increasing north- 
wards and westwards, respectively, and r is either physical the absolute physical distance 



in m or relative to Rq as shown in Figure 3.27 



3.6.1.2 Carrington 

While Stonyhurst is fixed with respect to the Sun-Earth, line Carrington coordinates rotate 
at the mean solar rotation rate with the first rotation commencing on 9 November 1853 and 
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North pole 




East limb 



West limb 



South pole 



Figure 3.27: A diagram of the Sun showing the Stonyhurst coordinate system (Thompson 
[20f6| ). 

being sequentially numbered thereafter. Both of these systems have limitations, especially 
in representing features far off the solar disk. 

3.6.2 Heliocentric Coordinates 

Heliocentric coordinates express the true position of a feature in terms of physical units 
from the centre of the Sun. Heliocentric Aries Ecliptic (HAE), Heliocentric Earth Ecliptic 
(HEE) and Heliocentric Earth Equatorial (HEEQ) are all examples of heliocentric coor- 
dinate systems consisting of three mutually perpendicular axes X, Y and Z. HEEQ is 
closely related to Stonyhurst heliographic by 



e = imV-\ZHEEQ/^Xl^EQ + y^^^g), (3.23) 

$ = arg{XHEEQ,YHEEQ) 



Xheeq = r cos(G) cos($) (3.24) 
Yheeq = r cos(9) sin($) 
Zheeq = r sin(e) 

No single perspective solar observation can be truly accurately transformed into he- 
liocentric coordinates. Sun-Observer lines-of-sight are not truly parallel, so the physical 
position will depend somewhat on the distance along a line-of-sight, also no distinction is 
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Figure 3.28: A diagram of the Sun with hnes of constant Hehocentric cartesian (x, y) overlaid. 
The z axis would point out of the page ( Thompson^ ^2006 j ) . 



made between the possible projections. However even with these problems Heliocentric 
coordinates are useful, especially for representation of the results of 3D reconstruction and 
form the basis of the helioprojective coordinate system. 

Both of the coordinate systems discussed below are observer specific i.e. non-terrestrial 
observers will measure different coordinates to those measured from Earth. As such, in- 
formation about the position of the observer must be provided to properly define the 
coordinate system. 

3.6.2.1 Heliocentric-Cartesian 

Heliocentric-Cartesian coordinates are a true cartesian system with each of the axes being 



mutually perpendicular, and all lines of constant x {y or z) being parallel, see Figure 3.28 



The z-axis is defined to be parallel to the observer-Sun line and pointing towards the 
observer. The y-axis is defined to be perpendicular to z and in the plane containing the 
solar North pole axis and increasing northwards. The x-axis is perpendicular to both x 
and y and increases westward. The distance along each axis is expressed in a physical 
distance of meters, or relative to Rq 



3.6.2.2 Heliocentric-Radial 

Heliocentric-Radial coordinates share the same z axis as the previous coordinate system, 
but X and y are replaced with p and (j). The impact parameter (j) is the radial distance 
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from the z-axis, again expressed in a physical unit or relative to Rq. The position angle (f) 
is measured in degrees counterclockwise from solar north as viewed by the observer. 

3.6.3 Projected Coordinate Systems 

It has already been stated that observations from a single viewpoint can only approximate 
heliocentric coordinates. The observations are projected against the celestial sphere and 
this needs to be accounted for, in order to accurately describe the observations. Helio- 
projective coordinates mimic heliocentric coordinates but replace physical distances with 
angles. All angles are measured from Sun centre as seen by the observer. As the coordinates 
are spherical in nature they take into account exactly what projections are used. 

3.6.3.1 Helioprojective-Cartesian Coordinates 

This is the projected version of heliocentric cartesian coordinates where x and y have 
been replaced by 6x and 9y. The helioprojective equivalent of z is C = -C© — d surface 
of constant C, represent spheres centred on the observer, with the sphere of ^ = passing 
through the centre of the sun. The relationships between helioprojective-cartesian and 
heliocentric-carteian are 

X = dcos{6x) sin{6y), (3.25) 
y = dsm{9y), 

z = Dq — dcos{9x) cos{9y), 

d = V^^TFTCD^^^^, (3.26) 
9x = arg{DQ - z,x), 
9y = sin~^(y/d). 

where Dq is the observcr-Sun distance, and d is the observer feature distance. Close to 
the Sun, where the small angle approximation holds, the relationship reduces to 




(3.27) 



(3.28) 

These are the default coordinates for STEREO SECCHI observations. 
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3.6.3.2 Helioprojective-Radial Coordinates 

The hclioprojective version of helcentric-radial is where the impact parameter p is replaced 
with the angle Op. 
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Chapter 4 

Data Analysis and Techniques 



In this chapter the methods used to enhance and analyse the STEREO observations are 
described. CMEs can be extremely faint compared to the coronal background and in order 
to detect and track them image processing methods are applied to enhance the CMEs sig- 
nature. Once the CMEs are identified, 3D reconstruction techniques are applied to the two 
perspective views to derive their true position. Using the true position of the CME derived 
from a number of observations, the kinematics can be calculated and compared to the var- 
ious models. Due to the large number of ill- determined parameters in the drag models a 
rigorous statistical procedure was used to test if the drag models appropriately reproduce 



the observations. This chapter is based on methods described in Maloney, Gallagher 



McAteer. Solar Physics, 2009 and Byrne, Maloney, McAteer, Refojo & Gallagher, Nature 



Communications, \2010. 
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4.1 Image Processing 

In coronagraph observations, CMEs are observed as outwardly moving regions of stronger 
brightness intensity relative to the background corona. As CMEs can be extremely faint 
compared to the background, imaging processing techniques are applied to enhance and so 
enable their identification and tracking. Coronagraph observations are dominated by the 
K and F corona with an occasional star, comet or planet overcoming this signal. Thus the 
main goal in processing coronagraph images ia the removal of the static coronal signal, and 
the enhancement of the CME signal. This is accomplished through a number of techniques: 

• Background Subtraction: Long term (days to weeks) background subtraction 
removes instrumental stray light and static coronal signals such as the F and K 
corona. The backgrounds are derived by taking the average, mean, minimum, or 
some combination of these operations applied to the series of images. Short term 
(hours to days) variations can also be subtracted to remove further semi-transient 
features. 

• Running Difference: Running difference images are created by subtracting the 
previous image from the current image and thus help in enhancing moving features. 
Formally this can be defined as: 

I'{x,y)i = I{x,y)i - I{x,y)i-i (4.1) 

where Il{x, y) is the running difference intensity at pixel (x, y), I(x, y)i is the current 
pixel intensity, and I{x,y)i-i the previous pixel intensity. 

• Base Difference: Base difference images are formed by removing a base image from 
every subsequent image, usually the base image is a pre-event image. 

• Normalised Radial Gradient filter: It is well known that the intensity of radia- 
tion falls steeply falls moving away from the Sun, by a factor ~10^ between the limb 



and ~3i?0. The normalising-radial-gradient filter (NRGF) (Morgan et a/., 2006) 
attempts to account for this radial falloff using: 

l'(r,</<)= [/(r,<^)-/(r)(^)]Mr)(^) (4.2) 

where I'{r,(j)) is the NRGF intensity at height r and position angle (j), I{r,(j)) is 
the original intensity, and /(r)^^^, (T{r)(^^\^ are the mean and standard deviation of 
intensities calculated over all position angles at height r 
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The standard procedure used was to reduce the data as described in Section [33] to produce 
total brightness (B) images from which backgrounds are then subtracted. The images were 
then spacially median filtered and smoothed, and at this point the various different types of 



image were created. The effects of processing CORl observations are shown in Figure 4.1 



Figure 4.1 (a) shows the reduced data, with instrumental stray light apparent even on 



the occulter, as is the dominant K and F corona. Figure 4.1 (b) has had a background 



subtracted and a binary mask applied to remove all non-imaging pixels and some of the 



artefacts know as "nails". Figure 4.1 (c) shows the effect of removal of an additional short 



term background from (b) while (d) shows a base difference image, (d) a running difference 



image and (f) NRGF applied to (c). Figure 4.2 shows the same but for C0R2 observations. 



4.1.1 Image Processing for the Heliospheric Imager 

In the case of HI, the images are dominated by the static F corona, and upon removal of 
this stars, the HI can image down to a 13**^ magnitude star. The CME signal per pixel is 
often comparable to the background and far smaller than the contribution from stars. Thus 
the removal of the star field is critical for detecting CMEs far from the Sun. The exposure 



times for HIl, HI2 observations are 1,200s and 4,950s respectively (Section 3.4.1.3). The 
star field in the field-of-view of the instrument can be very well approximated as static over 
this time period, but there will be a shift (~1 pixel per image in HIl and HI2) due to the 
motion of the satellite. A standard running difference image would then be contaminated 
with large signals due to this moving star field (see Figure [43]^c) and Figure [44| ^c)). If the 
offset between two subsequent images can be found, then the images can be shifted and 
subtracted and only transient features will be left. Formally this can be defined as: 

I {x, y)^ = I (x, y)^-I{x + a,y + /3) ..^ (4.3) 

where a and /3 are the x and y offsets due to the motion of the satellite. However, there is 
a geometrical optical distortion in the HI images which makes this difficult. We developed 
a method to approximate this offset. A region (512x512) is extracted from both images 
and a local cross-correlation is performed in a running window (128x128) around these 
sub-images. The average offset from this procedure is used to shift the image. In a similar 
manner to the NRGF the HI images can be flattened to account for the radial falloff in 
intensity. The NRGF cannot be directly used as the HI images only sample a small range 
of (f). However, by summing down the columns of a HI image, and then median filtering 
this with respect to time a pseudo NRGF can be derived which flattens the HI images as 



shown in Figure 4.3 (d). Figure 4.4 (d) shows the results of this modified running difference 
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Figure 4.1: A CORl observation from 08 November 2007 05:52UT. Different processing 
steps are shown: (a) reduced data, (b) processed and standard background subtraction, (c) 
additional background subtraction, (d) base difference, (e) running difference, and (f) NRGF 
of (c). 
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Figure 4.2: A C0R2 observation from 09 November 2007 01:22UT. Different processing 
steps are shown: (a) reduced data, (b) processed and standard background subtraction, (c) 
additional background subtraction, (d) base difference, (e) running difference, and (f) NRGF 
of (c). 
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on HIl and HI2 observations. The effects of the distortion are especially clear in the HI2 
images. 

4.2 Three Dimensional Reconstruction 



The concept of stereoscopy was first suggested by Wheatstone (1839). The reconstruction 
of three dimensional (3D) information from images of an object observed from two different 
perspectives is well studied and the general reconstruction methods are well developed. 
The main technique used in this thesis is known as 'tie-pointing' and relies on the epipolar 
constraint. This states that any feature identified on an epipolar line in one STEREO 



image must lie along the same epipolar line in the other stereo image ( Inhester , 2006 ) . An 
epipolar line is the projection of the epipolar plane (the plane containing the two observers 
and the point P) in the stereo images (see Figure [45| ) . 

4.2.1 Tie-pointing 

One possible implementation of this method of stereoscopy is outlined below. We implicitly 
assume that two simultaneous observations of a feature (the point P; see top panel of 



Figure 4.6) from two different perspectives are available. If the feature is identified in one 



of the images (from A for example), the angles and 9y (bottom panel Figure 4.6) can 



be derived, and the only unknown is the depth of the feature or the distance along the 



line-of-sight d (see bottom panel of Figure 4.6). One can arbitrarily pick two values for d 



{do, di) thereby giving two sets of 3D coordinates [6^, 9y, do] and [9^, 9y, di] , which can 
then be transformed into Heliocentric-Cartesian as described in Section [3.6.3.11 

These coordinates (of the line-of-sight), as viewed from the second perspective, can 
then be derived by reversing the transform, and the epipolar line can be over-plotted on 
the second image. If the position of the same feature along this epipolar line can be found, 
this gives all the information needed to complete the reconstruction. The same procedure 
is carried out for this second point, selecting two values of d and transforming these into 
the reconstruction coordinate system. Now the problem is reduced to finding the point of 
intersection of two lines in 3D, which can be simplified to solving for the intersection point 
of two lines in 2D in two planes, for example the x — y and x — z plane (see bottom panel 



of Figure 4.6). 

Due to the nature of the HI instruments, there are often only observations available 
from one spacecraft as a result, an additional constraint is needed. In the HI FOV, the 
CMEs are at large distances from the Sun so we assume that the CME propagates pseudo- 
radially, continuing along the trajectory that it followed in the CORl/2 FOV. Based on this 
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Figure 4.3: A HIl observation from 10 November 2007 09:29UT. Different processing steps 
are shown: (a) reduced data, (b) standard background subtraction, (c) running difference, (d) 
'radial' filtered, (e) modified running difference ( |4.3| ), and (f) intensity histogram of (b) the 
solid vertical lines show the thresholds used to display (b) and the dashed line is the mean of 
image (b). 
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Figure 4.4: A HI2 observation from 12 November 2007 06:10UT. Different processing steps 
are shown: (a) reduced data, (b) standard background subtraction, (c) running difference and 
(d) modified running difference (4.3). 



assumption, we use the best fit of the trajectory CORl/2 data in the x — y (echptic) plane 
and constrain the CME to propagate along this fit. This fit line is treated exactly the same 
as the line-of-sight above. The equations are solved for the point of intersection between 
this fit line and the observed line-of-sight, yielding the x, y position. The z coordinate 
is then calculated by assuming two distances along the line-of-sight to find the equation 
of the line with respect to x or and substituting in the relevant coordinate calculated 
yields the corresponding z value. 
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Figure 4.5: Orientation of epipolar planes in space and the respective epipolar lines in the 



stereo images for two observers looking at the Sun ( Inhester 2006 ) 



4.2.2 Elliptical Tie-pointing 



The tie-pointing technique is best suited to study smah features close to the Sun, such 



as coronal loops ( Aschwanden et al. , 2008 ) . When applied to study CMEs a number of 



the assumptions do not hold. Firstly, as CMEs are large and curved the lines-of-sight 



will not intersect upon it but ahead of it, see Figure 4.7 and Figure 4.8 (a). Also, due 
to the motion of the CME and the Thomson scattering geometry, the same part of the 
CME cannot be rigorously identified from both perspectives. Consequently CMEs cannot 
be fully reconstructed by tie-pointing alone, though their position in 3D space can be 



localised by the intersection of sight lines tangent to the CME front (de Koning et al 



2009 Pizzo & Biesecker 2004) 



Characterising the CME front as an ellipse from both view points allows this localisation 
to be carried out for a series of epipolar planes or slices. Each slice corresponds to a 
quadrilateral in 3D. As the CME is known to be a curved object, an ellipse is inscribed in 
the quadrilateral such that it is it tangent to each side (line-of-sight), giving an optimal 
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Figure 4.6: (top) 3D representation of the geometry of the observations with the hnes-of-sight 
(LOS) from the two spacecraft (sohd) , spacecraft Sun centre hnes (dashed) and the point to be 
reconstructed P. (bottom) Projection in the x — y plane showing the angles and ^ where 
X corresponds to image pixel axis. The point P is found by solving the equations of the LOS 
for their intercept. The same can be done in the x — z plane for By. 
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reconstruction of tiie CME front in that slice (Byrne et al., 2010). By repeating this for 
a number of shces the entire CME can be reconstructed. Then, applying this method to 
series of observations, allows the 3D reconstruction to be studied as a function of time. 



For a more detailed discussion of this technique see the PhD thesis by Byrne (2010). 



4.3 Drag Modeling 



In Section 2.2.2 a number of possible drag models were shown and a generalised drag 



equation (2.52) was written: 



dvn 



dr 



oirR~ 



(4.4) 



where Vgw = /(i?, Wsty^™^*"*^^). There are a number of solar wind models which could be 



used in the drag model, for example the Parker solar wind as described in Section 1.1.3 
but the one chosen for this work was based on coronagraphic measurements of blobs in the 



solar wind (Sheeley et al, 1997): 



asymptotic /i _ (R-2.8)/8.1 
„,,, V J- c 



(4.5) 



It should be clear that this is a non-linear ordinary differential equation and has no an- 
alytical solutions. In oder to evaluate this equation, a numerical integration scheme was 
used, 4**^ order Runge-Kutta (RK4). This method is outlined below: 



f{t,y) where y{to 



2/0 



dy 
dt 



(4.6) 



then the RK4 solution to this problem is: 



1 



Vn+i =yn + -^{ki + 2A;2 + 2k3 + k^) 



tn+l = tn + h 



(4.7) 
(4.8) 
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where i/n+i is the RK4 approximation at tn+i and: 



h = hf{tn,yn), 
h = hf{tn + ^h, yn + ^ki), 
ks = hf{tn + ^h, yn + ]jk2), 
ki = hf{tn + h,yn + ks). 



(4.9) 
(4.10) 

(4.11) 
(4.12) 



The k are the changes in the y values estimated from the slope at the begining (A;i), 
midpoint {k2, /C2) and end (/C4), which are then averaged with greater weight given the 
estimates at the midpoint. The RK4 method is a 4*^* order method so the accumulated 
error is 0{h^) while the error per step is 0{h^). Figure |4.9| 



shows the results of this 



numerical integration of (4.4) for a number of values of a and (3 in the both the quadratic 
{5 = 2) and linear cases {5 = 1), as described in Section 2.2.2 



4.4 Bootstrap Technique 



The bootstrap technique (Efron & Tibshirani 1993 ) is part of a broader group of techniques 



knows as resampling methods. These methods allow a series of 'numerical' experiments to 
be carried on data already gathered. The basic idea is to take original data, -F, consisting 
of n independent and identically distributed (i.i.d.) points and resample it to obtain m 
sub samples and m estimates of 6, the parameter(s) of interest. 



F {X : Xi,X2, ...,Xn),i 



{X : Xi,Xi+i,Xj) 01 
F2 {X : Xi,Xi+i,Xj) 6*2 



(4.13) 



The simplest resampling method is known as the 'jack-knife'. In this method, given a set 



of observations from some distribution X = (xi, 2:2, 2:3, 



5 Xr^ 



each data point is in turn 



removed from the set and the measurements of interest the mean and standard deviation 
for example, are recalculated and stored. This allows n — 1 resamples to be computed. 
The next resampling method is knows as the 'm-out-of-n' or 'deleted-D jackknife' where 
m out of n data points are removed this allows (^) resamples to be carried out, however 
caution must be taken to ensure m is not made too large. The methods outlined above are 
generally applicable when the data is from a distribution, as the samples need to be i.i.d. 
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However there are other methods which can be apphed to non-i.i.d. data such as 'cross 
vaHdation'. In 'cross-validation' the data is spHt up into roughly k equal parts for each of 
the A;*'* parts, fit all other parts and compute 9*. Or, if testing a model, the model can be 
fitted to the other A; — 1 parts and then the prediction compared to the A;*'*. 

4.4.1 The Bootstrap 

The bootstrap can be thought of as an extension to the jackknife though its mathematical 
foundation and insight into the data are very different. Instead of removing m data points 
we randomly sample, with replacement, from the original sample making a data set of the 
same size. The bootstrap technique is outlined below 

1. Collect n data points to form X = xi, X2, ■ ■ ■ , Xn 

2. Randomly resample with replacement n from X to form X* the bootstrap sample. 

3. Compute parameters of interest 9* 

4. Repeat steps 1 and 2 N times 

N is chosen to be very large 1,000 to 10,000 depending on the requirements. 

4.4.2 Linear Regression 

In linear regression the problem is generally to find the correct values of some predictor 
to match the responses, for example the kinematic constants for a ball under parabolic 
motion. The data sets for linear regression consist of two parts: the predictor q and the 
response yi, and together they can be written 

Xi = (cj,yi), (4.14) 

where Ci is a 1 x p vector Cj = (cji, Ci2, ■ ■ ■ , Qp). For a given set Xj, the conditional 
expectation given the predictor is 

Hi = E(yi\ci) (i = l,2,...,n). (4.15) 

The term 'linear' applies to the key assumption that the expectation value /Xj is a linear 
combination of the components of the predictor Cj or, 

10 

fii = Ci/3 = ^ CijPj. (4.16) 
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(3 is the parameter vector or regression vector, /3 = /32, • • • , and the goal of 
regression is to determine these values. Going back to a ball under parabolic motion, the 
regression parameter could have the form (3 = (/io,fo,«o)) and in this case the predictor 
would be Cj = (1, \/2t^). This may seem to be nonlinear but again the linear term applies 
to the form of the expectation value (a linear combination of Cj) and the fact that is it a 
quadratic function of t doesn't matter. 



Regression problems are usually expressed as: 

Vi = caPi + Cj2/32 + • ■ ■ CipPp + €i = Ci(3 + ej for z = 1, 2, . . . , n (4.17) 

where are assumed to be errors from a random sample of a unknown distribution F with 
a mean or expectation value of 0, 

F (ei, 62, . . . , en) = e where ^^(e) = 0. (4.18) 

The goal is to estimate the parameter vector /3 from the pairs of (ci, yi), (ci, yi), . . . , (c„, y„). 
If we assume a trivial value of b for /3 the residual squared error (RSE) can be found 

n 

RSE(b) = J](y,-Cib)2. (4.19) 

i=l 

Finally the least-squares estimate of (3, or the value of b that minimises the RSE giving, 
/3' is 



RSE(/3') = min 

b 



.i=l 



(4.20) 



An extremely usefully consequence of the linear regression formulation is that we can write 
it as a scries of matrix operations. If we set a matrix C an n x p matrix with its i*'* row 
Cj (the design matrix), and let y be (y^, j/i, . . . , Vn)^-, then we can write 

C'^CfS = C^y, (4.21) 

which gives /3' as 

f3' = (C^C)-iC^y. (4.22) 
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4.4.3 Bootstrap Linear Regression 



The model for linear regression has two components: the parameter vector /3, and F the 
probability distribution of the errors, 

P = {(3,F) (4.23) 

The general bootstrap algorithm outlined above requires that we estimate P, we have 
least-squares estimates of /3, (3' but no measure of F. If /3 was known we could directly 
calculate the errors (or residuals) ei = yi — Cj/3 but we don't know /3. However we can get 
an approximate error distribution using /3', 

e^ = yi- Cil3'. (4.24) 

The simplest approximation to F is the empirical distribution of e^, so F' is this distribution 
with a probability of 1/n and each residual. 

With our estimate of the error distribution F' we can now calculate P' = {/3' , F') or 
P' — > X*, which must be the same as P ^ x. To generate the x* we select a random 
sample of bootstrap error terms 

F'->(el,6^,...,e;) = e*, (4.25) 

where each is e* is a randomly drawn from e'. The bootstrap responses are then generated 
according to 

y* = aP' + e* (4.26) 

notice that is a constant for all data points. Finally, the bootstrap least-squares estimate 
/3'* is the b that minimises of the RSE of the bootstrap data, 



5;;(y*-c,/3')* = min 



i=l 



1=1 



(4.27) 



Errors and confidence intervals may be constructed from the bootstrap distributions of the 
of 13*. 

It should be clear that there are two possible ways to use the bootstrap, with linear 
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regression the Xj 



(cj, Ui) pair in which case the bootstrap sample is: 



— {{^h 5 yii)i {^121 2/12)) ■ ■ ■ ) i^ini yin)} 



(4.28) 



or bootstrap the residuals in which case the bootstrap sample is: 



{(ci, ci/3' + e',J, (ci, ci/3' + , . . . , (ci, ci/3' + e^J} 



(4.29) 



where 11,12, ■■■ ,in are random samples of the integers from 1 to n. The choice of which to 
use depends on the problem at hand, but in the limit as n, the number of samples grows 
large, both tend to the same result. 

4.4.4 General Bootstrap fitting 

Bootstrapping the residuals can also be applied in a broader manner in the following way: 

1. Make an initial fit of a model y = /(x, 0) to the data yi (this need not be via 
least-squares), and obtain fit parameters 6 

2. Construct an empirical residual distribution e^ = y^ — /(x, 0') 

3. Construct abootstrap residual sample e* by randomly sampling from the equally 
weighted empirical distribution e'. 

4. Make a bootstrap data set y^ = /(x, 0') + e* 

5. Fit a model to the bootstrap data y* = /(x, 0*), and obtain bootstrap estimate of 



6. Repeat steps 3-5 N times 

Where the fitting of the model to the data can be carried out by any number of means. 
With computational power of modern computers N is often chose to be very large, e.g. 
10,000-100,000. 
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observer 1 




observer 1 




observer 1 



Figxixe 1: Backprojection to reconstruct point-like, curve-like and surface-like objects to clemonstrate the 
different conditions of solvability. 



2 



Figure 4.7: A schematic of tie-pointing a curved 3D surface within the epipolar geometry. 
As the Unes-of-sight are tangent to different surfaces they may not insect upon it hmiting the 
accuracy of the reconstruction (Inhester 2006). 
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Figure 4.9: The results of numerically integrating the drag equation using input values for a 
number of possible scenarios. 
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Chapter 5 

CME Trajectories in Three 
Dimensions 



Before STEREO, single view-point observations obscured the true trajectories of CMEs. 
The STEREO mission provides two perspective views of CMEs in the inner Heliosphere. 
This, for the first time, enables the reconstruction of the 3D trajectories of CMEs in a range 
of heliocentric elongations (^2-88 degrees). In this work, a number of CMEs, which were 
simultaneously observed by CORl and C0R2 from both STEREO satellites, were selected. 
These observations were then used to derive the CMEs trajectories in 3D out to ^15 Rq. 
The reconstructions using CORl and 2 observations support a radial propagation model. 
Assuming pseudo-radial propagation at large distances from the Sun, the CMEs 3D posi- 
tions were extrapolated into the HI field- of-view. This enables the 3D trajectories of CMEs 
to be reconstructed from ^1.^- 2^0 Rq. It was demonstrated that CMEs undergo accel- 
eration in the inner Heliosphere; CMEs slower that the solar wind are accelerated, while 
CMEs faster than the solar wind are decelerated, both tending to the solar wind velocity. 
This effect will have a significant impact on the forecasting of CME arrival times at 1 AU 



and hence space weather. This chapter is based on work published in \Maloney, Gallagher 
& McAteer] Solar Physics, \200^ 
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5.1 Introduction 



Prior to the launch of STEREO, CMEs could only be routinely observed in white light up 



to 32 Rq in the plane-of-sky using the SOHO LASCO (Section 3.3.2), and rarely to Earth 



using the Solar Mass Ejection Imager (SMEI) instrument (Howard et al., 2007 Jackson 



et al. , 2004). CME related shocks are observed in radio data, but rely on density models to 



relate the observations to heights in the corona, and provide no information on the direction 
of propagation. Some work has been done using radio observations from two observatories 



to reconstruct 3D positions from radio triangulation (Reiner et al.: 1998 2007). Another 



possibility to track CMEs over large distances is during fortunate quadrature observations 



such as SOHO-Sun- Ulysses (Bemporad et al. , 2003). CMEs have also been tracked out into 



the Heliosphere when they interact with various spacecraft such as Wind, ACE, Ulysses, 



Cassini, Voyager 1 and 2 (Gopalswamy et al, 2005b, Lario et al, 2005, Richardson et al. 



2005, and references therein). 

The identification of events over interplanetary time/distance scales when these are 
not available requires the comparison of coronagraph observations with in situ data. The 
relationship between images and the in situ data is complex: coronagraphs image solar 
radiation which has been Thomson scattered by electrons in the CME and in-situ mea- 
surements provide the actual densities and magnetic fields of a track through the CME. 
All of these problems are compounded when multiple events occur in a short time frame. 

In an attempt to overcome some of these problems, statistical studies have been carried 
out to try to extract physics from the observations, however these types of study can miss 



vital event-to-event variations ( Gopalswamy et al. , 2005a ) . In a recent paper Owens Sz 



Cargill (2004), compared a number of empirical models to predict the arrival time of 



CMEs at 1 AU. They cited the main sources of error as projection effects, assuming that 
solar wind conditions are the same for all CMEs, and the determination of the point of 
intersection of the CME with the spacecraft with respect to in situ measurements. 



Coronagraph observations are subject to projection effects (Section 3.1.2) which are 



one of the leading sources of uncertainty in determining kinematics ( Howard et al. , 2008b ) 



Other methods of deriving height from observations such as the "point-P" , "fixed-(?!>" , and 
"harmonic mean" have been used, but a recently study of these methods on simulated data 
shows the problems associated with them. Figure \5A\ shows the height, error in the height, 
and velocities derived from synthetic observations from the various methods, from which 
it is clear none of them accurately produce the CME kinematics. 

Various techniques have been developed to overcome some of these limitations such as 
tomographic techniques which use the fact that the observations are taken in a rotating 
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Figure 5.1: Position (top), error (center) and speed (bottom) of a simulated CME and as 
derived from ttie syntlietic SECCHI images from a 3D MHD simulation with the different 



methods (Lugaz et al. 2009) 



reference frame to reconstruct the 3D nature of the corona (at the loss of time resolution; 



Prazin (2000)). Another possibility is the polarimetric technique of Moran & Davila (2004). 



Other methods that use a priori knowledge of the system can also be use to infer the 3D 
position, such as pre-assuming the CME geometry. This type of analysis, known as forward 
modelling, has been applied both to single view point and stereoscopic observations. In 
the later case, the model parameters are modified in either a semi-automatic or manual 



fashion to simultaneously fit the observations from both perspectives ( Boursier et al. 2009 



Thernisien et al. 2009). 



Stereoscopy has been used in astrophysics for a long time, with one of the earliest ap- 
plications being the determination of the distances to near-by stars using parallax. Using 
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stereoscopy techniques, the reconstruction of the 3D coordinates of features identified in 



images from two vantage points is a well-defined linear problem (Inhester, 2006). The 
biggest challenge in using this method is the identification and matching of features in the 
stereo image pairs. Stereoscopy has been successfully applied to on-disk features such as 



coronal loops ( Aschwanden et al. , 2008 ) and extended to prominences and CMEs ( Liewer 



eTaTI [20091 |Mierla et aT\ [20081 [SHVastava et al\ [20091 |Temmer et al.\ [2008t [Thompson 



2011). Also, assuming CMEs travel at constant velocity at large distances from the Sun, 



the trajectory can be derived by inverting the equation for the fixed-(/> case ( Sheeley et al. 



2008 ) . For a recent review of the methods which have been applied to coronagraph obser- 



vations see 



Mierla et al. (2010). 



In order to better understand the acceleration and propagation of CMEs the true 3D 
trajectory needs to be known to facilitate comparison between observations and theory. 
From a space weather perspective, knowing the 3D trajectory will lead directly to better 
predictions and also indirectly through a better understanding of CME acceleration and 
propagation. 

In this chapter the trajectories four CMEs are reconstructed in 3D. The observations 



are presented in Section 5.2 The true 3D trajectories of the four events are presented 
Section [5.3[ Section 5.4 contains a discussion of the results, the conclusions drawn, and 
future work. 



5.2 Observations 



The dates of the observations are 2007 November 18, 2008 March 25, 2008 April 9 and 
finally 2008 April 10. The CME catalogue hosted at the Coordinated Data Analysis Work- 
shop (CDAMfl [Copalswamy eral^ ( [2009[ )) provided heights derived from SOHO/LASCO 
observations which were deprojected and used to corroborate the accuracy of the 3D re- 



construction. Figures [5.2[ [5.3[ [5.4[ [5.5[ show the reduced and processed observation from 
CORl, C0R2, HIl and HI2. 



5.2.1 Uncertainties in the Reconstructed Heights 

In order to compare our results to theory we need to have an estimate of the uncertainties 
in our reconstruction. We use a statistical approach to find the spread of the data which 
we can then use as an estimate of the uncertainty. The 3D height data is de-trended with 
a third order polynomial and the moments of the distribution of the residuals are then 

' http:/ /cdaw. gsfc.nasa.gov/CMEJist 
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Figure 5.2: Sample observations from 2007 October 08-13 event from: (a) CORl Behind 

13:45 UT, (b) CORl Ahead 13:45 UT, (c) C0R2 Behind 22:23 UT, (d) C0R2 Ahead 22:22 UT, 
(e) Hll Behind 2007-10-10 15:29 UT, (f) HI2 Behind 2007-10-13 04:10 UT. 
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Figure 5.3: Sample observations from 2008 Marcii 25-28 event from: (a) CORl Behind 19:15, 

(b) CORl Ahead 19:15 UT, (c) C0R2 Behind 22:22 UT, (d) C0R2 Ahead 22:22 UT, (e) HIl 
Ahead 2008-03-26 06:09 UT, (f) HI2 Ahead 2008-03-27 14:09 UT. 
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Figure 5.4: Sample observations from 2008 April 09-10 event from: (a) CORl Behind 
11:25 UT, (b) CORl Ahead 11:25 UT, (c) C0R2 Behind 13:22, (d) C0R2 Ahead 13:22 UT, (e) 
HIl Behind 2008-04-10 05:29 UT. 
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Figure 5.5: Sample observations from 2008 April 09-13 event from: (a) CORl Behind 

15:05 UT, (b) CORl Ahead 15:05 UT, (c) C0R2 Behind 23:22 UT, (d) C0R2 Ahead 23:22 UT, 
(e) HIl Behind 2008-04-11 09:29 UT, (f) HI2 Behind 2008-04-12 18:09 UT. 
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calculated. The standard deviation (cr) of the distribution is taken as a good measure of 
the uncertainty. As a check on the statistics, the number of data points greater than one a 
from the mean were calculated and compared to the 31.8% which we expect if the data is 
normally distributed. If the number of points in the individual fields-of-view are too small 



this approach will not work. In Figure 5.6 we show data from the 2007 October 09 event 
for C0R2 and HI2 as examples. The top plots show the height-time data as well as the 
best fit line using a 3"^*^ order polynomial. The bottom plots show the de-trended data, 
mean, and one standard deviation {a) either side of the mean. The value of the la and 
the number of points greater than la from the mean are indicated on the plots. 



5.3 Results 

The observed apex of each CME was identified using point and click methods on the data. 
These points were then used to calculate the 3D coordinates and thereby plot the 3D 
trajectories. From these we were able to calculate the CMEs launch angles with respect to 
the ecliptic plane, and in the ecliptic plane with respect to the Sun-Earth line. To highlight 
the advantages these new data give, we estimated the velocity of the CME as viewed 
from A and B using a plane-of-sky assumption, as well as from the 3D reconstruction. 
These estimates were calculated using the first and last data point in the C0R2 FOV and 
v = Ah/At where h and t are the heights and times of the first and last observations. We 
also estimated the velocities on larger time/distance scales by using the first C0R2 data 
point and the last data point (from HIl or HI2) in the same manner. 

Using the derived launch angles observations from other spacecraft can be de-projected 
and compared to the heights derived from the 3D reconstruction. Figure |5.7| shows one 
such comparison of the deprojected heights derived from LASCO C2 and C3 taken from 
the CDAW catalogue and STEREO C0R2 and HIl. Within errors (approximately the 
symbol size) the reconstructed and de-projected heights agree in the range of the range 5- 
18 Rq. This demonstrates the CME reconstruction method is suitable and gives accurate 
CME height and launch angles. 



5.3.1 Event 1: 2007 October 08-13 

This event was first observed in CORl B on 2007 October 8 at 10:15 UT off the western 
limb, but was extremely faint in images from both spacecraft so the CORl data were not 
used for the 3D reconstruction. The first data point used for the 3D reconstruction was 
observed at 15:52 UT on the same day by C0R2 A and B instruments, the last data point 
was observed at 04:01 UT on 2007 October 13 by the HI 2 B instrument. From the 3D 
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Figure 5.6: Height-time plots for the 2007 October 09 event in the C0R2 and HI2 FOV, the 
hne is the best fit to the data with a 3"^ order polynomial (top row). The de-trended data, 
standard deviation a and the percentage of points out side of a one a error (bottom row) . The 



dashed lines are plotted at ± Icr from the mean (Maloney et al. 2009). 



data, the launch angle in the x-y plane was found to be 56° (where 0° corresponds to an 
Earth-directed CME, and 90° to a CME off the western limb), thus this was a front side 



event (see Figure 5.8). The CME also had an out of ecliptic component of 11° (where 0° 
corresponds to in-ecliptic while 90° corresponds to Solar North) . The velocity seen from A 
was found to be 190kms~^, from B 157kms~^, and from the 3D reconstruction 216kms~^. 
The mean velocity over the entire event was found to be 430 km s^^. The CMEs true height 
was tracked from 5.7-239.4 Rq (Figure [5^ top right) with an estimated error of 0.15 R0, 



0.25 Rq (Figure 5.6, bottom left) and 2R0 (Figure 5.6, bottom right) in C0R2, HIl and 
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Figure 5.7: STEREO 3D CME Height from C0R2 and HIl and LASCO CDAW deprojected 
height using the propagation angle derived (from the 3D reconstruction) for the CME. The 
X-axis is the time of day in hours. 



HI2 respectively. The CMEs morphology evolves from a nearly circular loop in CORl B 
and C0R2 A and B, to an elliptical structure with much larger extent in the vertical than 



the horizontal direction in HIl (see Figure 5.2). In HI2, the front evolves from nearly linear 



to a concave front where the wings/flank of the CME are ahead of what was the CME 
apex (see Figure [53]) . There is some indication of a multiple loop CME in C0R2 A, HIl 



B, and HI 2 B (see Figure 5.2). 



5.3.2 Event 2: 2008 March 25-27 

The second event was first observed on 2008 March 25 at 18:55 UT in both CORl A and 
B off the eastern limb. The last data point used was observed at 02:09 UT on 2008 March 
27. The CME was well observed in A and B CORl but was somewhat faint and diffuse in 
C0R2. The CME launch angle in the x-y plane was found to be —74° which corresponds 
to a front-side event, while the angle with respect to the ecliptic plane was —22° (see 



Figure 5.9). The velocity from A was found to be 863 kms ^, from B 1009 kms ^, and 
from the 3D reconstruction 1020 km s^^. The velocity over the entire event was found to 
be 494 km s"^ The CMEs height was tracked from 1.9- 139.7 R0 with an estimated error 
of 0.3 R0 in HIl, and 0.6 R© in HI2. There were not enough data points in the CORl/2 
FOV to estimate the error. The morphologies from both perspectives are very similar a 
nearly circular profile consisting of one faint loop which propagates all the way out into 



the COR 2 FOV (see Figure 5.3). In the HIl FOV there appears to be a second loop-like 
feature behind the main front which is expanding at a similar rate as the main front. The 
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Figure 5.8: CME trajectory for the 2007 October 08-13 event (Figure [5^: Cut in the x- 
y plane (top left), cut in the y-z plane (top right), cut in the x-z plane (bottom left), 3D 
trajectory (bottom right) from Maloney et al. (2009). 



CME is tracked until about halfway across the HI2 FOV before it becomes too difficult 
to identify. The second loop-like structure could be associated with a second flux rope, or 
perhaps prominence material. 



5.3.3 Event 3: 2008 April 09-10 

This event was first observed on 2008 April 08 in both CORl A and B at 10:05 off the 
western limb toward the South. The last data point was observed at 15:29 UT on the 2008 
April 10 in HIl B. The CME is well observed in A and B images from both CORl and 
C0R2. Due to the large out-of-ecliptic angle the CME does not pass through the FOV of 
HI2 B. The CME launch angle in the x-y plane was found to be —71° which corresponds to 
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Figure 5.9: CME trajectory for the 2008 March 25-27 event (Figure 5.3): cut in tlie x-y plane 
(top left), cut in the y-z plane (top right), cut in the x-z plane (bottom left), 3D trajectory 
(bottom right) fromlMaloney el a/.l (|2009|). 



a front side event, while the out-of- ecliptic plane angle was —60°. The velocity from A was 
found to be 504kms~^, from B 513 kms~^, and from the 3D reconstruction 542 kms~^. The 
velocity over the entire event was found to be 354 km s"^ The CMEs height was tracked 
from 1.9- 50.4 Rq with an estimated error of 0.5 in HIl. There were not enough data 
points in the CORl/2 to estimate the errors. The CME appears to be narrower in CORl 
A than in CORl B. When the CME enters the C0R2 FOV it appears the same size in 



A and B (see Figure 5.4). The CME has a more elliptical shape in HIl and continues to 
propagate at a steep angle. There is no obvious trend of the trajectory back towards the 



ecliptic plane, which one may have expected (see Figure 5.10). 
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Figure 5.10: CME trajectory for the 2008 April 09-10 event (Figure 5.4): cut in the x-y plane 
(top left), cut in the y-z plane (top right), cut in the x-z plane (bottom left), 3D trajectory 
(bottom right) fromlMaloney et al.\i2009\. 



5.3.4 Event 4: 2008 April 09-13 

This event was first observed off the eastern limb on the 2007 October 9 at 15:05 UT by 
both CORl instruments. The last data point was observed by HI2 A at 00:09 UT on 2008 
April 13. The CME was extremely faint in both A and B CORl, but a small number of 
data points were able to be recorded. The CME was reasonably clear in C0R2 A, but very 
faint in B. As a result the CME front was not visible and for this event a dark linear feature 
was tracked instead. The CME launch angle in the x-y plane was found to be —51° which 



corresponds to a front-side event, while the out-of-ecliptic angle was 16° (see Figure 5.11). 



The velocity from A was found to be 109 km s ^, from B 217 km s ^, and from the 3D 
reconstruction 189 km s^^. The velocity over the entire event was found to be 338kms~^. 
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Figure 5.11: CME trajectory for the 2008 April 09-13 event (Figure 5.5 ): cut in the x-y plane 
(top left), cut in the y-z plane (top right), cut in the x-z plane (bottom left), 3D trajectory 
(bottom right). 



The CMEs height was tracked from 2.2- 128.0 R0 with an estimated error of 0.8 Rq in 
HIl and 1.1 R0 in HI2 (there were not enough data points in the CORl/2 to estimate the 
errors). The CME appears very narrow in C0R2 B images, and could be a loop structure 
with its axis nearly perpendicular to the plane of sky. The CME has an elliptical shape in 



HIl with some material in the centre (see Figure 5.5). This shape is held into HI2 where 
the CME is lost as it crosses in front of the Milk Way, where the background subtraction 
technique does not work as well. 
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5.4 Discussion and Conclusions 



We have shown it is possible to reconstruct the 3D trajectories of CMEs in the inner He- 
hosphere using STEREO/SECCHI observations. These 3D reconstructions give the true 
position of the CMEs studied in the range between ~2-240 R© as a function of time. The 
CORl/2 data indicates that a radial propagation model is appropriate once the CME is 
above a few solar radii. This is in agreement with what we expect, as the forces acting on 
CMEs fall off steeply with increasing radius (especially forces which would act perpendic- 
ular to the propagation direction in the ecliptic). We have also shown that the velocities 
that would be derived from either of the spacecraft alone can be very different to that of 
the 3D velocity. Also we have shown that CMEs do undergo acceleration in the Heliosphere 
as the velocities in the C0R2 FOV can be very different to those calculated using the first 
and last data points. The reconstructions can also be of significant value for space weather 



forecasting. Using an empirical formula such as used in Gopalswamy et al. (2005a) we can 
try and estimate the arrival time of these CMEs at 1 AU. The arrival time of a CME at 
one 1 AU is given by: 

r = (5.1) 

a 

where 

T 

with 1/(1 AU) the velocity at lAU, U the velocity in low corona, r the time between the 
two measurements, and d in this case 1 AU. Now as we have no in-situ data at 1 AU we 
use the fact that in all the events the CMEs have tended towards a typical solar wind 
value. This should be a good proxy for the velocity at 1 AU. Using the 3D velocity and the 
velocity that differs the most, we can make estimates of the CME arrival time r. For the 
first event r was calculated to be 5.39 days using the velocity derived from the B spacecraft, 
and 5.11 day from the 3D velocity. The difference between the two calculated arrival times 
for the events are Ar =0.28 days, Ar =0.39 days, Ar =0.17 days, and At =0.10 days for 
the four events respectively. Unfortunately, as none of the CMEs were detected in-situ 
near 1 AU we can not tell if these differences correspond to improvements or not. Another 
major advantage is the ability to predict if the nose or flank of the CME will interact with 



Earth -this was identified as one of the major sources of error in Owens & Cargill (2004). 
The geo-effectiveness will be affected by which part of the CME interacts with the Earth, 
and also the arrival times will be affected as the CME flank should take longer to arrive 
at Earth. 

There are a number of outstanding issues with the reconstructions. Some of these are 
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the results of the reconstruction method, while others are related to the data used or its 
interpretation. The stereoscopy method requires that we can identify the same feature in 
both images and this is extremely difficult in some cases as the CME is very faint or has 
a very different appearance from each spacecraft. For example, CORl observations for 
the first event were so faint the CME could not be identified in the A images (se Figure 



5.2), and in both the third and fourth events the identification of the same feature was 



extremely difficult as the CME has a different appearance in each image (see Figures 5.4 



and 5.5). This is evident in the reconstructions as they are less contiguous, and there 
are large changes in the CME trajectory which do not make physical sense as there is no 
known source for such large forces out in the Heliosphere. 

The different sensitivities and light rejection levels of the different instruments mean 
the CME may appear bigger or smaller just due to instrumental effects. It is believed that 
the discontinuities of the first event are due to this type of effect. Also as the front of the 
CME imaged is a line-of-sight integration through the CME (flux rope or magnetic bubble), 
the front seen from each spacecraft will be different. This is especially important in the HI 
FOV as CMEs can expand up to large sizes and we only have one point of view, but will 
also have an effect in CORl/2. This means the positions, and thus heights derived, will 
always be an upper limit on the actual values. If a model shape for CMEs is used - such 
as an ellipse or circle - we can estimate the size of this effect. A CME which passes one of 
the spacecraft taking in situ measurements would allow a comparison of the reconstructed 
CME path and the actual CME path and provide a welcome test of the methods. 

A number of interesting observations arise from these data, such as the possible multi- 
loop nature of two of the CMEs, the complex structure seen at the rear of the CMEs 
(especially in HIl), and finally the distortion of the CME front in Event 1. Firstly events 
1 and 2 show possible multi-loop structures (see Figures 5.2 and 5.3): these could be 
attributed to a second flux rope or possibly due to prominence material. These features 
are only visible because of the unique features of the HIl and HI2 instruments. Secondly, 
most of the events show complex structure at the rear of the CME. This is very apparent 



in Figures |5.2[ |5.4| and |5.5[ This type of structure could be interpreted as the rear part of 
a flux rope if the flux rope axis was into, or out of, the image plane. Finally, the distortion 
of the CME front in event 1 as it progresses through the Heliosphere could be caused by a 
number of processes. This may be a Thomson scattering effect i.e., we are seeing different 
parts of the CME from different planes, and it only appears as though the CME front is 
distorted with the flanks ahead of the nose. The other possibility is that as the CME has 
expanded to extremely large size at this point and the nose and flanks are experiencing 
very different ambient solar wind conditions. As such, the CME flank may be in the high 
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speed solar wind region (hence accelerated by aerodynamic drag), or the nose may be in 
a relative density-enhanced region due to the heliospheric current sheet and experiencing 
more drag than the CME flanks. This effect has previously suggested to occur and is 
known as 'pancaking'. 

These reconstructions are the first steps towards detailed study of CME kinematics in 
the heiloshpere. Using these trajectories, we can derive the full 3D CME kinematics which 
will not be subject to projection effects. When we combine these with CME mass estimates 
we will have all the information needed to study the forces acting on CMEs from ~2 Rq 
to beyond 1 AU. We plan to compare the derived forces with those predicted by the flux 



rope, 'snow plough', and aerodynamic drag models (Cargill, 2004 Vrsnak, 2006 Vrsnak 



& Gopalswamy, 2002). Having the 3D information will finally allow for a full comparison 



between the various theories and the data. In doing so we will be able to gain insight into 
some the interesting observations discussed here. 
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Chapter 6 

CME Kinematics in Three 
Dimensions 



The forces governing CME evolution during both their early acceleration, and later prop- 
agation are still unclear, although plane-of-sky coronagraph measurements have provided 
some insight into their kinematics near the Sun (<32Rq). The dual perspectives of the 
STEREO spacecraft are used to derive the 3D kinematics of CMEs over a range of he- 
liocentric distances (^2 - 250 Rq). Evidence for solar wind (SW) drag forces acting in 
interplanetary space, with a fast CME decelerated, and a slow CME accelerated, towards 
typical SW velocities was found. It was also found that the fast CME showed a linear 
dependence on the velocity difference between the CME and the SW, while the slow CME 
showed a quadratic dependence. The differing forms of drag for the two CMEs indicate 
the forces and thus mechanisms responsible for their acceleration may be different. Fur- 
ther, using an advanced 3D reconstruction technique and statistically rigorous fitting, we 
were able to show that a CME underwent aerodynamic drag at distances greater than 7Rq. 
This chapter is based on work published in Maloney & Gallagher , Solar Physics, 2010 and 
Byrne, Maloney, McAteer, Refojo & Gallagher , Nature Communications, 201(\ 
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6.1 Introduction 



The kinematic evolution of CMEs can be broken into three phases: initiation, accelera- 



tion, and propagation (Zhang et al., 2001). During the acceleration phase magnetic forces 



such as the Lorentz force are thought to dominate and drive the eruption (Section 2.2.1). 
During the propagation phase, the initial acceleration has ceased, and the CME motion is 
dominated by the interaction between the SW and the CME (Section 2.2.2). The "snow 
plough" , aerodynamic drag, and flux-rope models all aim to explain the motion of CMEs 
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2010). 



Before the launch of the STEREO mission, synoptic white-light CME observations were 
hmited to 32 Rq using LASCO, while SMEI sometimes tracked CMEs to Earth (~215 Rq). 



In radio observations, fast CMEs which drove shocks could be tracked to Earth (Reiner 



et al. 2007). Interplanetary Scintillation (IPS) observations provided density and velocity 



measurements for both CMEs and the SW from 50 Rq to beyond 1 AU, and using tomo- 
graphic techniques, can give 3D information ( |Manoharan[ |2006l |2010[ ). CMEs are also 
observed in in situ measurements with WIND and ACE at LI (~1 AU), and occasionally 
CMEs can be tracked up to very large distances of up to 5 AU using additional spacecraft 
( |Tappin[ [20061 ). 

Numerical modelling has been used to study CME propagation with numerous ap- 
proaches such as ID Hydro simulations, 2.5D MHD simulations, and full 3D MHD sim- 



ulations ( 
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Cargill et al. (1996), Cargill & Schmidt (2002), and Cargill (2004) used MHD simulations 



to show that aerodynamic drag was an appropriate description of a flux rope in a mag- 
netohydrodynamic flow, and found that the drag coefficient was normally around unity, 
but can range between approximately 0.0-10.0. He also found that there was a feedback 
interaction between the flow and flux rope which could significantly distort the flux rope 



depending on the strength of the magnetic field (see Figure 6.1). Smith et al. (2009) us- 
ing Hakamada-Akasofu-Fry version 2 (HAFv2; Fry et ar]|2003 ), a modified kinematic code, 
simulated a series of events, predicting up to 89% of the observed events. Falkenberg et al. 



(2010) using the ENLIL (Odstrcil &: Pizzo, 1999a| [b) code with a cone model, which is a 
3D time-dependent MHD code, modelled an observed CME and compared the simulation 
to the measured in situ data, and found it showed good agreement given the correct input 
parameters. While ENLIL is a MHD code, the CMEs are simulated by time dependent 
velocity, density and temperature perturbations at the inner boundary, with no magnetic 
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Figure 6.1: Magnetic field lines projected in the x-z plane (a), (c), and (e). Contours of 
constant Bd = By — Byo/Byo where Byo is the initial external field strength (b), (d), and (f). 
Solid and dashed lines represent Bd less than and greater than zero. Time is normalised to 



the Alfven wave crossing time of the flux rope — 2ao /vao ( Cargill et al. 1996 ) . 
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field component similar to HAVv2. The fact both ENLIL and HAVv2 do not model the 
CME magnetic field but still give accurate predictions may indicate that the magnetic field 
is not crucial in determining CME dynamics, while Cargill's results indicate it is for CME 
morphology. 

Statistical studies comparing in situ with white light observations indicate a trend of 



CME velocity converging towards the SW velocity as they propagate to 1 AU ( Gopalswamy 



2006). Other studies based on white light observations have indicated that aerodynamic 



drag of some form may explain this trend (Shanmugaraju et al. 2009: Vrsnak, 2001 ). Radio 
observations suggest that a linear form of aerodynamic drag may be most appropriate for 



fast CMEs ( Reiner et al. , 2003 ) . Tappin ( 2006 ) showed that acceleration can continue 



far out (5AU) into the Heliosphere. However, these studies are subject to the difficulties 
associated with the observations they are based on. For example, white light observations 
were limited to single, narrow, fixed, view-points meaning only observations of the inner 



Heliosphere could be made, and even these were subject to projection effects (Howard 



et al., 2008b). Also, linking features in imaging and in situ observations is complex and 



can be ambiguous, a problem exacerbated during periods of high activity. In the case of 
numerical simulations their complexity can make it hard to extract which effects are the 
most important, possibly obscuring the important underlying physics. 

Using STEREO observations, a number of papers have been published which extract 
3D information and study CMEs at extended heliocentric distances, over-coming some of 



the difficulties outlined above. Davis et al. (2009) identified a CME in HIl and HI2, and 



using a constant velocity assumption ( Sheeley et al. , 2008 1 they derived the speed and 



trajectory of the CME. The predicated arrival time, based on the speed derived, agreed 



with the in situ observations. Wood et al. (2009) used the "point-P" and "fixed-t/)" methods 



to derive the height, speed, and direction from elongation measurements out to distances 
of ~12Oi?0. A recent paper by|Liu et al] (|2010[) tracked a CME to ~15Oi?0 in 3D using 



J-maps from both spacecraft to triangulate the CMEs position in 3D. On the other hand. 



Maloney et al. (2009) tracked the trajectory of CME apexes in 3D using triangulation. 



some as far as 240 Rq. 

A general equation describing the motion of a CME may be written: 

Dv 



P 



Dt 



j X B -VP + pg + Fd 



(6.1) 



where the first term on the right is the Lorentz force, the second term is the force due to 
gas pressure, the third term is gravitational force, and the final term is the drag force. An 
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equation describing the motion of a CME in the drag dominated regime may be written: 



M- 



-l/2CDPswAcme{v 



(6.2) 



where Cd is the drag coefficient, psw is the solar wind density, Acme is the CME area 
and Vsw is the solar wind velocity, and M is the CME mass. We use a parametric drag 



model similar to that of Vrsnak Sz Gopalswamy (2002) with the added parameter 6, which 



determines if the drag is quadratic (aerodynamic) or linear (Stokes) (see Section 2.1.2). 
This parametric form collapses the complex dependences of the CME area (Acme), mass 
(M) and the solar wind density (psw) into a power-law which depends on heliospheric 



distance R. Thus (6.2) can be written: 



dv 
dt 



-aR-^iv-Vsn^Y 



(6.3) 



where a, 13, and 6 are constants. 

The CME trajectories were derived using the tie-point method (Section 4.2.1) and from 
these the kinematics were calculated. These kinematics are then used to investigate the 
effects of drag on the CME. We present the reconstructed CME kinematics for four events: 
one acceleration (CME 1), one decelerating (CME 2), and one with constant velocity (CME 
3). For the final event (CME 4), the entire CME front is reconstructed, and we consider 
the midpoint kinematics of the reconstructed front this event showed a more complex 



acceleration profile. The observations are described in Section 6.2. 1[ and their analysis is 



described in Section 6.2.2 The results for each event are presented in Section 6.3 Finally, 



a discussion of the results followed by the conclusion is presented in Section 6.4 



6.2 Observations and Data Analysis 
6.2.1 Observations 

The trajectories of four CMEs were reconstructed using observations from STEREO SEC- 
CHI observed during: 2007 October 8-13 (CME 1), 2008 March 25-27 (CME 2), 2008 
April 9-12 (CME 3), and 2008 December 12-15 (CME 4). The observations were re- 



duced and processed as described in Section 3.5 and Section [4Tj Figure [6^2] shows reduced 



observations from the 2008 March 28 event where the CME is simultaneously observed in 
both CORl and C0R2 in both the Ahead and Behind spacecraft, but only in from the 



Ahead spacecraft in HI. Figure 6.3 shows observation of the 2008 December 12 event. 



CMEs 1-3 were observed in either the inner coronagraph (CORl) or outer coronagraph 
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Figure 6.3: (a) Indicates the STEREO spacecraft locations, separated by an angle of 86.7 at 
the time of the event, (b) Shows the prominence eruption observed in EUVTB off the northwest 
limb from approximately 03:00 UT, which is considered to be the inner material of the CME. A 
multiscale edge detection and corresponding ellipse characterisation are over-plotted in CORl. 



(c) Shows t he CME propagation towards Earth from STEREO-A and STEREO-B ( [Byrne 
efaIl|20T0l). 
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(C0R2) simultaneously by both STEREO- A and STEREO-B. From these images the CME 
apex was localised via tie-pointing (Section |4.2.1 ), and by tracking it through a series of 
images the trajectory was built up. In these events the CME was only observed in HI 
by one spacecraft, so an additional constraint is required to localise the CME apex: we 
assumed that the CME will continue along the same path, with respect to solar longitude, 
as it did in the CORl/2 field-of-view. CME 4 was observed by both A and B by ah SECCHI 
instruments. For this event the entire CME front was reconstructed using the elliptical 
tie-pointing (Section 4.2.2). From the 3D data the CMEs kinematics were derived. 



6.2.2 Kinematic Modelling 

The kinematics were only fit during the time interval we believe that drag is at play and 
the observations are accurate. There was evidence for early acceleration which was not 
attributed to drag and as a result these data were not fit. Also, some observations of CMEs 
tracked far into HIl or HI2 field-of-view were excluded from fitting as identification of the 
CME became ambiguous. 

The kinematics of CMEs 1 and 2 were fitted, via a least-squares method, with a para- 



metric model for the drag (6.3). In order to test which form of drag is most suitable (linear 
or quadratic), we fitted (6.3) with 5 set to 1, and then separately with 5 equal to 2. A 
number of the model parameters can be fixed from the observations, such as the CME 
height and velocity. We assume that the CME tends to the SW speed, which was taken to 
be where the velocity plateaus. The model parameters obtained from the fitting were then 



compared with previous results from Vrsnak Sz Gopalswamy ( 2002 ) . From this comparison 



we infer which model best reproduces the kinematics, and hence is the most appropriate. 
Both the fast and slow CMEs (CME 1 and CME 2) were analysed using this method. The 
intermediate CME (CME 3) was fit with a constant acceleration model as this seemed 
most apt since there was no apparent acceleration. 

In the case of CME 4 the kinematics were analysed using a bootstrapping method 



(Section |4.4[ ). An initial fit of (6.3) was carried out and used to calculate the residuals of 
the fit. These were then used to create the bootstrap samples by randomly sampling them 
and adding them back on to original fit, creating a bootstrap sample. All of the model 
parameters were allowed to vary within physically appropriate limits for a bootstrap run of 
10,000 iterations. The observationally known parameters such as CME initial height and 
speed were compared to the resulting bootstrap values as a verification of the method. 
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6.3 Results 



CME 1 and CME 2 were fitted in three ways with a, /3, and 5 all allowed to vary (black 
line) , with fixed 5 of two (magenta line) , and one (orange line) . In both cases the free fitting 



returned values that were not comparable to previous studies (Vrsnak, 2001) showing the 



need to apply so form of physical constraint. The constrained fit parameters for the events 



are given in Table 6.1, CME 3 was fit with a constant acceleration model (black line). 



CME 4 was fit using the bootstrap method with /iq, vq, Vsw, ct, /3, and 5 allowed to vary. 



6.3.1 CME 1 (2007 October 8-13) 



Figure 6.4 a)-(c) shows the kinematics for the accelerating CME. This CME was first ob- 
served at 15:05 UT on 2007 October 8 off the west limb, and was found to be propagating 
at an angle of 56° from the Sun-Earth line. Figure |6.4[ a) shows the height of the CME. 
Figure [6^ b) shows the velocity profile which clearly shows the CME is undergoing accel- 
eration, initial velocity of ~150 km s~^, and final velocity of ~450 km s~^. There may be 
two acceleration regimes; an early increased acceleration phase (before 18:00 UT on the 
October 8), followed by a drag acceleration. The early acceleration can be attributed to a 
magnetic driving force, and so was not fitted with the drag model. Later, when the CME 
reached the centre of the HI2 field-of-view, determining the front position became difficult, 
so this region was not fitted. Figure [6^c) shows the acceleration profile of the event. The 
5 = 2 (orange) fit gives the lowest chi-squared value. 



6.3.2 CME 2 (2008 March 25-27) 

The kinematics from the decelerating CME are shown in Figure [6^ d)-(f). This CME 
was first observed at 18:55 UT on 2008 March 25 off the east limb and was found to 
be propagating at an angle of -82° from the Sun-Earth line. Figure |6.5[ d) and (e) show 
the height and velocity profiles, and the velocity profile clearly demonstrates the CME is 
undergoing deceleration. The CME had an initial velocity (HIl) ~800kms^^ and final 
velocity of ~375kms^^. Due to the high speed of this CME, it was only observed in a 
small number of frames in CORl and C0R2. As a result, the kinematics were difficult to 
quantify in these instruments. However, there appears to have been an early acceleration 
feature. The deceleration in the HIl and HI2 field-of-view continued until the CME reaches 
a near-constant velocity, and travels at this velocity throughout the rest of the field-of-view. 
Figure [6l5| ^d) shows the acceleration profile of the event. The 6 = 1 (magenta) fit gives the 
lowest chi-squared value. 



161 



6. CME KINEMATICS IN THREE DIMENSIONS 




-40 r , 1 

09-Oct 10-Oct 11-Oct 12-Oct 13-Oct 
Start Tirr-e (08-0ct-07 12:00:00) 



Figure 6.4: Kinematics of the 2007 October 8 event, (a) height, (b) velocity, and (c) accel- 
eration. Vertical dashed line indicates separation between early and late phase acceleration. 
Horizontal dot-dash line indicates the assumed SW velocity. The black line corresponds to a 
fit varying a, {3 and 5, the magenta fit has 5 — 2 while the orange fit has 5 ~ \ 
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Start Time (25-Ma--08 18:00:00) 

Figure 6.5: Kinematics of the 2008 March 25 event, (a) height, (b) velocity, and (c) accel- 
eration. Vertical dashed line indicates separation between early and late phase acceleration. 
Horizontal dot-dash line indicates the assumed SW velocity. The black line corresponds to a 
fit varying a, j3 and 5, the magenta fit has 5 = 2 while the orange fit has 5 = 1 
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Figure 6.6: Kinematics of the 2008 April 09 event (a) height (b) velocity (c) acceleration. 
This event shows an early acceleration (to left if dashed line) which levels of to a scatter about 
typical solar wind speeds. This event was fit with a constant acceleration the resulting fit 
parameters are ho — 22Rsun, uo = 334kms~^, and a=— O.lSms^^. The assumed solar wind 
value is indicated by the horizontal dot-dash line. 
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Table 6.1: Fit parameters for the accelerating events. CME 1, Vgw = 450 km s Vc 
233kms~i , ft-o = 12 7?©. CME 2, Vsw = 325kms-\ v„ne = 702kms-i , /iq = AARq. 



6.3.3 CME 3 (2008 April 9-12) 



In Figure 6.6 we show the kinematics of the constant velocity CME. This CME was first 
observed at 15:05 UT on 2008 April 09 off the east limb and was found to be propagating 
at an angle of -73° from the Sun-Earth line. Figure |6.6[a) shows the height of the CME. 



Figure 6.6 ^b) shows the velocity profile which has a scatter about ~300kms~^. Again, 
there may be some evidence in the CORl/2 observations for an early acceleration phase, 
but due the event's poorly observable features at this early stage, it is hard to quantify 
this. The departure from the fit after April 12 20:00 UT is thought to be due to error in the 
reconstruction, n as the CME apex becomes too faint to identify. As this event shows no 
obvious acceleration it was not fitted with the drag model but with a constant acceleration 
model h{t) = /iq + ^o^ + ^f^at^ (thin black line). Figure 6.6 'c) shows the acceleration 



profile, the fit values (/iq = 22 R^-^n, fo = 334kms ^, and a = — 0.18ms ■^) are consistent 
with negligible acceleration throughout the field-of-view. 

6.3.4 CME 4 (2008 December 12-15) 

The prominence associated with this CME was observed at 50-55° north from 03:00 UT 
in EUVI-B images, obtained in the 304 A passband, in the northeast from the perspective 
of STEREO- A, and off the northwest limb from STEREO-B. The prominence is considered 
to be the inner material of the CME, which was first observed in CORl-B at 05:35 UT 



(Figure 6.3). This CME was Earth-directed thus ~0° from the Sun-Earth line. The CME 
was initially propagating at 40° north from the ecliptic, but due non-radial propagation this 
tended towards zero as the CME attained a height of ~5Oi?0 as shown in Figure [6?7| b). 
This non-radial motion was found to be well characterised by a power law. The CME was 
also found to expand, its angular width increasing from approximately 33° to nearly 70° 
(Figure [6.7[c)). Its expansion was also well characterised by a power law. Figure [6?7|a) 
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Figure 6.7: Kinematics of the 2008 December 12 event, (a) velocity, (b) declination, and 
(c) angular width. The velocity in (a) is derived from the CME midpoint and the fit os an 
aerodynamic drag model (Byrne et al. 2010). 
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shows the velocity profile for the event, the CME was initially rapidly accelerated from ap- 
proximately 100 -300 km s~^ between 2-5 i?0, before gradually rising to a scatter around 
a value of 550kms~^. The acceleration peaks at approximately 100 ms^^ at a height of 
~3i?0 and then decreases. The early acceleration is attributed to a Lorentz force, while 
the subsequent increase in velocity, at heights above 7 Rq , is predicted by theory to result 



from the effects of drag. From inspection of the Wang-Sheeley-Arge model ( Arge & Pizzo 



2000 Wang & Sheeley, 1990) it appears that the CME is in a high speed solar wind stream 



of ~550kms ^. 

The initial CME height, CME velocity, asymptotic solar wind speed, a, /3 and S are 



obtained from a bootstrap distributions for those parameters (Figure 6.8). These provide 



the best fit to the observations, as well as confidence intervals of the parameters. Best-fit 



values for a and /3 were found to be (4.55+|^?) x 10"^ and 



-2.02lgg5, which agree with 



values found in previous modelling work Vrsnak (2001 ). The best-fit value for the exponent 



of the velocity difference between CME and the solar wind, 5, was found to be 2.27^q jq, 
which is clear evidence that aerodynamic drag (6=2) functions during the propagation of 
the CME in interplanetary space. 

The drag model provides an asymptotic CME velocity of 555]J]42^ kms~^ when extrapo- 
lated to 1 AU, which predicts the CME to arrive ~ 1 day before the Advanced Composition 
Explorer (ACE) or WIND spacecraft detected it at the LI point. We investigate this dis- 
crepancy by using our 3D reconstruction to simulate the continued propagation of the 
CME from the Alfven radius (~21.5i?0) to Earth, using the ENLIL with Cone Model 
at NASA's Community Coordinated Modeling Centeij^ The height, velocity, and width 
from the 3D reconstruction were used as initial conditions for the simulation, it was found 
that the CME was actually slowed to ~342kms~^ at lAU. This is a result of its inter- 
action with an upstream, slow-speed, solar wind flow at distances beyond 5Oi?0. This 
CME velocity is consistent with in situ measurements of solar wind speed (~330kms~^) 
from the ACE and WIND spacecraft at LI. Tracking the peak density of the CME front 
from the ENLIL simulation gives an arrival time at LI of 08:09 UT on 16 December 2008. 
Accounting for the offset in CME front heights between our 3D reconstruction and ENLIL 
simulation at distances of 21.5- 46 ii0 gives an arrival time in the range of 08:09- 13:20 UT 
on 16 December 2008. This prediction interval agrees well with the earliest derived arrival 



times of the CME front (Davis et al. . 2009 Liu et al. . 2010). This was identified as plasma 



pileup ahead of the magnetic cloud flux rope in the in-situ data of both ACE and WIND 



(Figure 1.16). The CME subsequently impacted the Earth. 



^ http: / /ccmc. gsfc.nasa.gov/ 
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Figure 6.8: Bootstrap distributions from a 10,000 iteration run for the parameters of (6.3) 
namely: asymptotic solar wind speed, CME initial speed, CME initial height, a, /?, and 6. 
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6.4 Discussion &; Conclusions 



We have shown that it is possible to derive the 3D kinematics of features (the CME apex 
in this case) in the inner Heliosphere (~2-25Oi?0) using STEREO observations. The 3D 
kinematics are free from the projection effects of traditional 2D kinematics but may contain 
artifacts from the 3D reconstruction method (e.g, Maloney et al.pOOQ ) and other sources. 
Both CME 1 and CME 2 showed two regimes in their velocity profiles, a low down (< 15 Rq) 
early rapid acceleration (in comparison to later values), followed by a gradual acceleration 
far from the Sun (>3Oi?0). The early acceleration is thought to be due to a magnetic 
driving force, as the solar wind velocity low in the corona (f<jio(< 10 R0) <268kms~^, 
Sheeley et a/.||1997 ) is lower than the velocity already attained by the CMEs in both cases. 
Here we assume that the later acceleration is due to the interaction between the SW and 
the CME, as in each case the CME attains a final velocity close to typical values for the 
solar wind. 

Considering CME 2 in Figure [6^ d)-(f), it can clearly be seen that the velocity levels 
off to a constant value typical of the solar wind. We interpret this as the CME reaching 
the local solar wind speed, and as a result, the force acting on the CME going to zero. 
For CME 1 Figure [6^a)-(c), the velocity initially increases, however, there is a plateau 
towards the end after April 11 6:00 UT which occurs at SW like speeds. The height 
measurements towards the end are very scattered and show a rapid increase. This is 
most likely due to losing the front to the background noise and triangulating a different 
feature. CME 3 propagates at a roughly constant velocity, which is consistent with the 
drag interpretation. The CME appears to have already attained the local SW speed and 
therefore is not accelerated. The fitting results show that a linear dependence produces a 
better fit for the fast event (CME 2), while a quadratic dependence better fits the slow event 
(CME 1). The differing range of the interaction CME 1 ~12Oi?0 and CME 2 ~ 80 7^0 
may be explained by the suggestion that wide, low mass, CMEs are more affected by drag 



than narrow massive CMEs (Vrsnak et al, 2010). 



Reiner et al. ( 2003 ) suggest that for fast events, a linear model of drag better reproduces 



the kinematics, which agrees with our findings. Vrsnak (2001) also suggested that a linear 



dependence might be appropriate. From a theoretical perspective, a quadratic dependence 
corresponds to aerodynamic drag, while a linear dependence suggests Stokes' or creeping 
drag. It is not currently clear which model is more physically correct. The fit parameters 



obtained do not agree with those found by Vrsnak (2001) and while our values are not 



unphysical, it is not clear why they differ so much from the previous studies. 
The mechanism behind the apparent differing forms of drag, linear (5 



1) and 
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quadratic {5 = 2), for the slow and fast event are unclear. The application of any hy- 
drodynamic theory such as drag, to a CME, may be missing vital physics. Could the 
magnetic properties play a role modifying the form of the drag (reconnection, suppression 
of turbulence, wave energy transport)? For example, Cargill et al. (1996) showed that, 
depending on the orientation of the flux rope and background magnetic field (aligned or 
non-aligned), the drag coefficient can vary between 0.0-10.0. They also found that the 
magnetic field of the flux rope is important for its survival as it propagates. Further, which 
form of drag is correct for a CME in the SW, the low Reynolds number viscous dominated 
Stokes' drag, or the high Reynolds number turbulence dominated aerodynamic drag? In 
order to address these questions a larger sample study is needed in order to verify that 
these effects are a recurring and observable phenomena and also to build up the statistics. 

In one case, that of CME 4, a bootstrap analysis allowed us to conclude that the 
CME acceleration was consistent with aerodynamic drag. However using this model, or a 
ballistic, model to extrapolate an arrival time at Earth to compare to the in situ data, gave 
an arrival time about a day too early. An ENLIL simulation of the CME indicated that the 
CME was decelerated by a slow speed wind stream ahead of the CME, at distances greater 
than 50 Rq . The importance of the dynamic interaction between CMEs and the solar wind 
is highlighted by the initial acceleration and subsequent deceleration of the CME. 

We have shown that it is possible to derive the true 3D kinematics for a number of 
CMEs in the inner Heliosphere. Based on this, we have been able to conclusively show 
that CMEs undergo acceleration in the inner Heliosphere, and more specifically, that due 
to its range and strength, this acceleration is believed to be the result of some form of drag. 
In one case we were able to show this acceleration was due to aerodynamic drag. Drag 
acceleration has important implications for space weather predictions, and for the analysis 
techniques which assume CMEs travel at constant velocity through the Heliosphere. The 
HI observations of CMEs in the Heliosphere provide a unique and limited opportunity to 
study the propagation of CMEs, and to understand the coupling between the solar wind 
and CMEs. 
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Chapter 7 

Direct Imaging of a CME Driven 
Shock 



Fast CMEs generate standing, or bow shocks, as they propagate through the corona and 
solar wind. Although CME-driven shocks have previously been detected indirectly via their 
emission at radio frequencies and measured directly in situ, direct imaging has remained 
elusive due to their low contrast at optical wavelengths. Here we report the first images of 
a CME-driven shock as it propagates through interplanetary space from 8Rq to 120 Rq 
(0.5 AU), using observations from the STEREO HI. The CME was measured to have a 
velocity of ^1000 km s~^ and a Mach number of 4-l^l-2, while the shock front stand-off 
distance (/S.) was found to increase linearly to ^20 Rq at 0.5 AU. The normalised standoff 
distance (A/Dq) showed reasonable agreement with semi- empirical relations, where Dq is 
the CME radius. However, when normalised using the radius of curvature, A/Rq did not 
agree well with theory, implying that Rq was under- estimated by a factor of ^3-8. This is 
most likely due to the difficulty in estimating the larger radius of curvature along the CME 
axis from the observations, which provide only a cross-sectional view of the CME. This 
chapter is based on work published in Maloney & Gallagher\ The Astrophysical Journal 



Letters, \20Tl 
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Figure 7.1: Diagram of the various quantities used to describe the shock and CME. 



7.1 Introduction 



Bow shocks occur when a blunt object moves relative to a medium at supersonic speeds 



( Rathakrishnan , 2010). These shocks are formed across many scales and in different con- 



ditions: from astrophysical shocks such as planetary bow shocks (Slavin Sz Holzer, 1981) 



or the shock at the edge of the Heliosphere (van Buren et ai, 1995), to shocks generated 



by the re-entry of the Apollo mission capsules (Glass 1977). CMEs which travel faster 



than the local fast magnetosonic velocity (with respect to the solar wind velocity) produce 
such standing shocks in the frame of the CME ( Stewart et al. . 1974a|b| ) . Interplanetary 
(IP) CME-driven shocks have previously been detected in radio observations as Type II 
bursts and using in-situ measurements. Direct imaging of shocks, on the other hand, has 



remained elusive, primarily due their low contrast ( Gopalswamy et al. , 2008 ; Vourlidas Sz 



Ontiveros 2009). 



In principle the jump conditions in combination with the MHD or fluid equations 



allow any shock to be modelled and understood (Section 2.3). However, the simulations 
are very complex, and it is often necessary and some times more insightful to work with 
analytical solutions or semi-emprical relations. There may also be unknown properties and 
assumptions have to be made about the underlying system such as it geometry. This is 
more so in the case of MHD shocks as the addition to the shock and flow properties, the 
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magnetic fields must be known. 



The sliape, size, and standoff distance of a shock are controlled by several factors: the 
shape and size of the obstacle; the velocity difference between the obstacle and the medium 
with respect to the sonic speed (i.e., the Mach number); and the properties of the medium, 
such as the ratio of specific heats (7) and the magnetic field. Imaging observations of a 
shock give no information about the magnetic field strengths or directions, or the flows 
speeds or directions. As a result the jump conditions are not useful in this case. However, 
imaging observations allow the speed of the CME to be measured, and, using a model, its 
relative speed to the background medium and Mach number can be derived. The standoff 
distance, CME shock and shock shape can also be readily measured from the imaging 
observations. 



Relationships between the shock standoff distance and the compression ratio have been 
derived by a number of different authors. The well known semi-empirical relationship of 



Seiff (1962) has the form: 



0.78 



Pu 



(7.1) 



A 

Do Pd 

which was derived for a spherical object, where A is the shock standoff distance, Dq is 
distance from the centre to the nose of the obstacle (in this case the radius), and pu, Pd 
are the densities upstream and downstream of the shock respectively. Using gas-dynamic 



theory, Spreiter et al. (1966) demonstrated that Pu/Pd could be written in terms of the 



upstream sonic Mach number, Ms, and the ratio of specific heats 7: 



A 



1.1 



(7-1)M2 + 2 



(7.2) 



(7 + l)M| ■ 

The increase in the coefficient in the standoff relations from 0.78 to l.l is due to the 



fact that the object under consideration (Earth's magnetosphere) in (7.2) is more blunt 



than a sphere; specifically, it is an elongated ellipse. Neither rt7.1^ or (7.2) behave as 



expected at low Mach numbers, where the shock should move to a large standoff distance. 



A modification which corrects for this enables (7.2) to be written in the form: 



A 



1.1 



(7-l)M,^ + 2 
(7 + l)(M|-l)^ 



(7.3) 



where the additional term in the denominator ensures the shock moves to a large distance 



as the Mach number approaches unity (Farris & Russell, 1994). They also suggested that 



using the obstacle radius of curvature rather than radius would be more suitable as it 
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accounts for the shape of the obstacle, resulting in: 



A 

Ik) 



0.81 



(7-1)M2 + 2 



(7 + l)(M2-l) 



(7.4) 



where Rq is the obstacle radius of curvature. 



In general, a conic section can be represented by y{x)^ = 2R{D — x) + b{D — x)'^ where 
b is the bluntness (6 < —1: blunt elliptic; b = —1: spherical; — 1 < 6 < 0: elongated 
ecliptic; 6 = 0: parabolic; b > 0: hyperbolic). The shape of the shock fronts are known to 
be represented by a modified conic section. One such parameterisation of the shock front. 



from Verigin et al. (2003a), is: 



2/2 (x) = 2Rs{Ds 



x) + 



X 



[1 + 



M|-l 



l + ds{Ds-x)/Rs' 



(7.5) 



where bs is the bluntness of the shock, and ds is related to the asymptotic downstream 



slope or Mach cone (see Figure 7.1). 



The relationships between the standoff distance and the Mach number have been inves- 
tigated from a number of perspectives, including numerical modelling, analytical relations. 



laboratory experiments and in-situ measurements of planetary bow shocks (Spreiter & Sta- 



hara, 1980 1995). These have shown that, in general, the semi-empirical relations provide 



an adequate description of shocks, with the low Mach regime being an exception (Verigin 



et al., 2003b). Depending on the physical context the sonic Mach number (Mg) can be 
replaced with the magnetosonic Mach number (Alms), when dealing with plasmas such as 
the solar wind and CMEs. It has been shown that using gasdynamic relations works well 
when dealing with magnetised plasmas when the MHD mach numbers are high. It also 
provides a good approximation when the Alfven (Ma) or fast magnetosonic (Mms) Mach 
numbers are low and these Mach numbers are substituted for the gasdynamic {Mg) Mach 



numbers (Fairfield et al., 2001). The reason this makes a good approximation can be seen 



by looking at the Rankine-Hugoniot equations or the jump conditions (Section 2.3): in the 
case of both gas-dynamic and MHD shocks the mass continuity relation is the same: 



[pUn] = 0. 



(7.6) 



This is the most important quantity in determining the standoff distance. The standoff 
decreases until the mass flux flowing around the body matches the mass flux crossing the 
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shock. 

Standoff distances of CME-driven shocks have been investigated from an in-situ per- 



spective by many authors (e.g., Lepping et al, 2008 Odstrcil et al.. 2005 Russell & Mul- 



ligan , [2002 ) . Russell & Muhigan (2002) found the shock standoff distance (A; thickness of 



magnetosheath) was of the order of 21 R© at 1 AU. Lepping et al. (2008) derived an average 
A of about 8 R© at 1 AU. However, when considering the CME radius (flux rope radius) as 
Do, the typical A expected from ([7^ is about SRq at 1 AU. iRussell & Mulliganl (|2002l) 



proposed that (7.4) may be more suited as it accounts for the fact the CME front may not 
be circular and that the radius of curvature at the nose is a dominant factor in determining 
the standoff distance. However, they found that (7.4) did not fit the observations either, 
and speculated this may be due to observational effect of only measuring one of the radii 
of curvature of the CME. The underlying structure of a CME is believed to be a flux rope: 
which has two characteristic curvatures; a smaller one due the curvature perpendicular to 
its axis (the radius when viewed as a cross-section), and the larger curvature along the 
axis. 

We investigate if the shock relations hold for a CME-driven IP shock. Specifically, we 
use direct observations of a CME-driven shock observed in C0R2 and HIl instruments 
of SECCHI on STEREO. In Section [T^H we present SECCHI observations of the CME 
and resulting shock, and describe the analysis technique. The results of our analysis are 



presented in Section 7.3, We discuss our results and state our conclusions in Section [7^ 



7.2 Observations and Data Analysis 



The CME analysed first appeared in the CORl (Section [sXO ) coronagraph images from 



STEREO-B at 15:55 UT on 2008 April 5. It was most likely associated with a B-class 
flare from NOAA active region 10987, which was just behind the west limb as viewed from 



Earth. Figure 7.1 shows the CME as it propagates out from the Sun into the different 
instruments' fields-of-view from 8R0 to 120 R©. The CME was visible in both A and B 



spacecraft in the inner and outer coronagraphs, but was only visible in HIl (Section 3.4.1.3) 
from STEREO B. The CME propagation direction was found to be ~106° west of the Sun- 
Earth line, the spacecraft were at a separation angle of 48° degrees (from each other). The 
shock is visible as a curved brightness enhancement in both the C0R2 observations in 



Figure 7.1 and also the HIl B observations. The assumption was made that the curved 
front is a shock, but there were no radio or in situ data available to corroborate. However, 
due to the CME's velocity (~1000kms~^) and the smoothness and position of the feature 



ahead of the CME, it can be argued that this is a legitimate assumption (Bemporad k, 
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Figure 7.3: (a) 3D reconstruction of the CME and shock front viewed perpendicular to the 
propagation direction with the initial ellipse fits overplotted. (b) Data transformed into a 
coordinate system centred on the initial ellipse fit. Overplotted is the subsequent fit to the 



shock front using ( 


7.5 


& Gallagher 


2011 





Units have been normalised with respect to Dq on the right ( Maloney 



Mancusol 120101 lOntiveros & Vourlidasl 12009 ). 



For each observation in which the CME or shock was visible the front was identified. 
A number of points along this front were then manually chosen. For the observations 
where the CME or shock was observed from both spacecraft, the front was localised in 



three dimensions using the tie-point method (Section 4.2). As the CME or shock was 
only observed by one spacecraft in the HI field-of-view, we used the additional assumption 
of pseudo-radial propagation, based on the direction derived from CORl and C0R2 to 



localise the front (Section 4.2). The resulting data consisted of a series of points in 3D for 
the CME and shock for each observation time. Figure 7.3 a) shows the 3D reconstruction of 
both the shock front and CME front, viewed perpendicularly to the direction of propagation 
(assumed to be a cross-section). The techniques for deriving the 3D coordinates of features 
in the CORl/2, and especially the HI field-of-view are not without error. In the case of 
the event studied here, the CME was close (<10°) to the plane-of-sky of STEREO A. As 
a result, errors in position should be small. 

In order to compare with relationships in Section |7.1[ the data were transformed into 
a coordinate system centred on the CME. To accomplish this, each CME front was fitted 
with an ellipse. The centre coordinates of these fits were then used to collapse all the data 
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on to a common coordinate system centred on the CME. The shock front was fitted with 



(7.5), which gave the shock properties such as the shock standoff distance A, the Mach 



number M, and the radius of curvature at the nose of the obstacle Rq- Figure 7.3 a) shows 



data and the initial fit, Figure [T^b) shows the shifted data and the shock fit using (7.5). 
The fast magnetosonic Mach number was calculated using Mms = (vcme — Vsw)/'Vmsi where 
Vcme is the CME velocity, Vsw is the solar wind velocity and Vms is the fast magnetosonic 
speed. Since Vsw and Vms were not known at the position of the CME, a model corona 
was used to evaluate them. This was based on the Parker solar wind solution with a 
simple dipolar magnetic field of the form B{r) = Bo^RQ/r)'^, where Bq was 2.2 G at the 
solar surface ( Mann et al. , 2003[ > . For each of the paired CME and shock observations the 
standoff distances A {=Ds — Dq) were obtained by three different means: (i) using the 3D 
coordinates of the furthest point (max(/i), where h = y^x^ + + z^) on the shock and the 
CME as hshk and heme respectively; (ii) the previous method applied to the data in the 
common coordinate system giving Dq and Ds', and (iii) the front fitting procedure also 
applied to produce standoff distances. However the results of method (i) cannot be used 
with the relations from Section 7.1 as they are not in a CME/obstacle centred coordinate 
system, though the results from method (ii) and (iii) can be compared to (7.2), (7.3) and 
(pl). 



7.3 Results 

A summary of the shock properties derived from the observations as a function of time is 
shown in Figure [7^ a)-(f). With the exception of the CME {heme) and shock heights (/i^/ifc), 
all the properties have been derived from the data, collapsed on to a common coordinate 
system with respect to the CME. The gap between the first three data points and the 
others is a result of both the CME and shock leaving the C0R2 field-of-view and entering 
the HIl field-of-view. The contrast between shock and background in the first three and 
last three observations is extremely low, making identification of the shock difficult. As 
a result, these points are not reliable, and should be neglected. Figure |7.5[ a) shows the 
derived heights of the CME and shock as they were tracked from 8 R to 120 R (0.5 AU) . 



Using a linear fit to hghk — heme (=A) versus heme (Figure 7.4), the extrapolated standoff 
distance at Earth was found to be ~54R0. Figure [T^ b) shows the distance to nose of 
the CME {Do) and shock {Dg) front directly measured (filled symbols); also shown are 
the values derived from fits to the shock and CME front (hollow symbols). The increasing 
offset between the two is due to the differing centres of their coordinate systems, as one 
is elliptic and the other is parabolic. Figure [T^c) shows the standoff distance A derived 
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Figure 7.4: Shock stand-off distance (A) against CME Height and a hnear fit to the data. 



using Do and (filled symbols) and from the fits to the fronts (hollow symbols). Both 
are in general agreement and show an increase with time. The standoff distance normalised 
using Do is shown in Figure [7^ d). The normalised standoff distance is roughly constant 
with a mean value of 0.37 it 0.09 Do~^- The standoff distance normalised to the radius of 
curvature at the nose of the CME (i?o) is shown in Figure [T^ e). The curvature could 
only be derived from the front fitting, and as such only hollow data points are shown. 
Figure |7.5[ f) then shows the magnetosonic Mach number (M^s) derived using: (i) the 
CME speed in conjunction with the coronal model (filled symbols); and (ii) the shock 



front fitted using (7.5) (hollow symbols). The mean Mach number from the coronal model 
was 3.8ib0.6, while a value of 4.4ibl.6 was found using the front fitting method. The mean 
Mach number from both methods was 4. lit 1.2. 



Figure [T^a) shows the relationship between the normalised standoff distance (A/Do) 
and Mach number (M^s) for a number of models. The Mach numbers were calculated using 
the coronal model (filled symbols), and front fitting (hollow symbols). The normalised 
standoff distances were calculated using measured values of Do and Ds (filled symbols). 
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Figure 7.5: Shock properties derived directly from the observations (filled symbols) and from 
fits to the shock and CME (hollow symbols) as a function of time, (a) The maximum height 
of the CME front (triangles) and shock front (circles), (b) The distance to front of CME (Dq) 
and shock (Ds) in CME centred coordinate system, (c) The shock standoff distance A. (d) 
The normalised standoff distance {A/Dq)- (e) Standoff distance (A) normalised by the radius 
of curvature of the CME (Rq). (f) The Mach number (M) derived the CME velocity and 
model for the corona (filled circles) and from the fits to the shock front (hollow circles) from 



Maloney & Gallagher (2011). 
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and fits to the CME and shock fronts (hollow symbols). Both show good general agreement 
between our observations and the models (<20%). The model of Seiff (1962) shows the 
poorest agreement, although this is not unexpected as it was derived for a circular obstacle 
and the CME is quite blunt compared to a circle. Figure |7.6[ b) shows the relationship 
between the standoff distance normalised by the radius of curvature of the CME (A/Rq) 
and Mach number for a number of models. In this case, Rq can only be derived from the 
front fitting. These values are then plotted as a function of the Mach numbers derived 
using both methods described above (hence, each value of A/ Dp appears twice). Our 
results do not agree with the expected relationship (7.4) and indicate that the radius of 
curvature Rq is underestimated by a factor of Ri3~8. One possible reason for this is that 
we have not considered the effect of the magnetic field of the CME, and solar wind effects 
on the shock. However, one would expect if this had a significant effect it would also affect 
the other relationship. It should be noted that the fast magnetosonic velocity and sonic 
velocity calculated from our model differ by less than 7% after excluding the first three 
data points, as mentioned earlier. This also suggests that the magnetic field should not 
play a major role. A more likely reason is due to an observational effect similar to that 



suggested by Russell & Mulligan (2002), where only one radius of curvature of the CME 
is observed. The observations provide a cross-sectional view of the CME along one of its 
axes. As a result, we have no information on the curvature along other CME axes. 



7.4 Discussion and Conclusions 

For the first time, we have imaged a CME-driven shock in white light at large distances from 
the Sun. The shock was tracked from 8 Rq to 120 R© (0.5 AU) before it became too faint to 
be unambiguously identifed. The CME was measured to have a velocity of ~1000kms~^ 
and a Mach number of 4.1ibl.2, while the shock front stand-off distance (A) was found to 
increase linearly to ~20 Rq at 0.5 AU. The normalised standoff distance [A/Dq) was found 
to be roughly constant with a mean of 0.37ib0.09 Dq^^- The normalised standoff distance 
derived using Dq and Ds-, and its relation to the Mach number [Mms)-, were compared to 
previous relationships and showed reasonable agreement. The normalised standoff distance 
[A/ Do) and Mach number were also derived by fitting the CME and shock front, which 
agreed well with theory and our other method of estimation. The fitting also allowed us to 
find the CME radius of curvature (Rq), enabling us to test the relationship between A/ Rq 
and the Mach number. In this case the derived ratios did not agree with the theoretical 
predictions and showed a significant deviation. 

The faint nature of the shock front made its identification challenging, and thus the 
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Figure 7.6: (a) Shock standoff distance normalised to Dq as a function of Mach number, (b) 
Shock standoff normalised by Rq as a function of Mach number. Also shown are the results 
of a number of semi-empirical models. Filled symbols indicate values derived using a coronal 



model, while hollow symbols indicate values derived using front fitting (Maloney & Gallagher 
20U\. 



front location and characterisation showed some scatter (Figure 7.5). For example, the 



Mach numbers in Figure [7^ f), show a large amount of variability especially from the front 
fitting. The standoff distances in Figure [T^ c) show the same trend, and the two different 
methods give similar results. It should be noted that the first three, and last three, data 
points show large deviations from the rest of the data for a number of derived properties. 
These correspond to very low contrast observations, and hence should be ignored. The 
Mach number derived from our coronal model and the CME position and speed agrees 
with that derived from the shock front fitting method. This is a good indication that our 
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methods accurately describe the shock even in the presence of large uncertainties. 

Both sets of data for the normalised shock standoff distance A/Do versus Mach num- 
ber (Mms) derived directly, and from front fitting, show good general agreement (Figure 
|7.6| ^a)). The standoff distance normalised by the CME radius of curvature [A/Rq) versus 
Mach number (Mms) from either the fits, or derived directly, do not agree with any of 
the relationships (Figure [7.6[ b)). Assuming that a CME can be modeled as a flux rope, 
it should have two radii of curvature. Our observations are a measure of a combination of 
these, which depends on the orientation of the fiux rope. This observational affect implies 
that we may only be measuring the smaller of the two, and this leads to the underesti- 
mation of Ro- Finally, the general agreement between the Mach number derived from 
our model and the Mach number derived from the fits suggest that the fitting is not the 
source of the problem. Using a mean Mach number of 4 to derive a value of 0.26 from the 
semi-empical ratio A/Rq, and the standoff distance at Earth calculated above (54 R©), we 
can estimate the radius of curvature of this CME at Earth to be ~2O7R0 (~0.95AU). 
This value could be tested using multipoint in situ measurements. 

Imaging observations of CME-driven shocks opens up a new avenue for studying their 
fundamental properties. This type of observation will be highly complimentary to radio 
and in situ measurements. A complete picture of the shock could then be constructed, and 
the derived properties from the different observations could be compared and contrasted. 
Furthermore, the analysis presented here will be applicable to future observations of shocks. 
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Chapter 8 

Discussion &: Future Work 



The goal of my thesis research was to increase our understanding of CME evolution in the 
inner Heliosphere. To achieve this goal, a 3D reconstruction technique was developed and 
applied to CMEs far from the Sun and shown to he effective. From these reconstructions, 
the true CME kinematics could he extracted and studied. These kinematics were compared 
to theoretical models and conclusions on their applicahility drawn. A CME-driven shock 
was imaged and its properties compared to semi-empirical models and shown to agree well. 
This chapter presents the main results and conclusions of this thesis and outlines possihle 
future work. 
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8.1 Principal Results 

The purpose of this thesis research was to exploit the unique observations of the STEREO 
mission to further our understanding of CME evolution, focusing especially on their prop- 
agation at large distances from the Sun. This was achieved by extending and applying 

3D reconstruction techniques to CMEs imaged by the SECCHI suite of instruments from 
both STEREO A and B. From these 3D reconstructions the true CME ki!d, free from the 
limiting uncertainties of projection effects. Investigating these kinematics revealed CMEs 
do undergo acceleration far from the Sun. The acceleration is consistent with some form 
of drag, and in one case could be shown to be due to aerodynamic drag. The advanced 
image processing methods used, allowed the detection and tracking of a CME-driven shock 
at large distances from the Sun. Analysing this, it was found that the shock stand-off 
distance agreed well with semi-empirical models. 

The principal results from these studies may be summarised as follows: 

8.1.1 3D CME trajectories 

• The STEREO observations allowed CMEs to be tracked from the solar surface to 
beyond 1 AU. The modified running difference technique allows the extremely faint 
CME signature, in HI ficlds-of-view, to be identified and tracked. This, for the first 
time, enabled CMEs to be tracked in white light from the Sun to the Earth. 

• The dual vantage point observations of STEREO were used to reconstruct the 3D 
trajectories of CMEs in the CORl and C0R2 field of view using a triangulation 
technique, namely tie-pointing. These reconstructions supported a pseudo-radial 
propagation model for CMEs at large distances from the Sun. This enabled the 3D 
reconstruction to be extended into the HI field of view, as often CMEs are only ob- 
served by either HI A or HI B due to their propagation direction. The assumption 
of pseudo-radial propagation provided the additional constraint necessary to trian- 
gulate the CME in the HI ficld-of-view. Thus, CME trajectories were reconstructed 
in the range between 2 - 240 Rq . 

• From examining the CME trajectories, initial estimates of the CME kinematics were 
made. These demonstrated that CMEs were being accelerated, with fast CMEs 
decelerated, slow CMEs accelerated, and both tending towards the solar wind speed. 
Previously, this acceleration could only be inferred, however using STEREO and the 
3D reconstruction methods it was directly shown. 
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• The morphology of CMEs far from the Sun was revealed and showed their complex 
structure - especially at the rear of some CMEs - and the possible multi-loop struc- 
ture of others. As CMEs propagate, their morphology evolves and, in one case, there 
was clear evidence for CME 'pancaking'. 

8.1.2 3D CME Kinematics and Drag Modelling 

• Using the 3D trajectories, the true kinematics of CME apexes were derived in the 
inner heliosphere, < 250 i?0, which were free from projection effects of traditional 
observations. 

• The detailed analysis of the CME kinematics revealed that CMEs often contain multi- 
ple phases of acceleration: low down (<15 Rq) rapid acceleration, and a more gradual 
acceleration far from the Sun (>15i?0). The rapid acceleration was attributed to 
magnetic forces, as the speed attained by the CMEs was already in excess of the solar 
wind at the same positions. The later acceleration was posited to be due to drag as 
the CMEs tended towards solar wind like speeds. 

• In order to test if this acceleration was due to drag the CME kinematics were fitted 
with drag models. In the case of a fast CME the kinematics were best reproduced 
by a linear drag model, while in the case of a slow CME a quadratic drag model was 
deemed best. The differing forms of drag for that the two CMEs indicated the forces 
responsible for their acceleration might have been different. 

• An entire CME front was reconstructed for the 2008 December 12 event using the 
elliptical tie-pointing technique. The kinematics of the midpoint of the CME front 
were extracted and studied. This CME was deflected into a non-radial trajectory 
which was well described by a power law. The angular width of the CME was found 
to increase super-radially and was also well described by a power law. The CME 
kinematics showed an early acceleration phase below 5 Rq , with the CMEs speed in- 
creasing from approximately 100 to 300 km s~^, the acceleration peaked at 100 m 

at ~ 3i?0. This was followed by a more gradual acceleration up to a scatter around 
550kms~^. The kinematics beyond 7 Rq were fitted with a parameterised drag equa- 
tion. Using a bootstrapping technique it was shown that this CMEs acceleration was 
due to aerodynamic drag. Using the final CME speed gave an estimated arrival time 
at LI a day too early. This was investigated by running the ENLIL with cone model 
simulation, constrained by the CME properties derived from the 3D reconstruction. 
The simulation showed the CME was decelerated by a slow speed stream ahead of 
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it, slowing the CME to approximately 330 km s - thus it predicted an arrival time 
within hours of the actual arrival time. 

• The dynamic interaction between CMEs and the solar wind was shown to contain 
multiple phases as the CME propagated through different solar wind regions. This is 
an important unknown factor in determining CME propagation. Accurate forecasting 
of CME arrival times without knowledge of the solar wind conditions the CME would 
encounter was shown to be difficult. 



8.1.3 CME-driven Shock 

• A CME-driven shock was identified and tracked from 8 Rq to 120 Rq (0.5 AU). Shocks 
have been imaged close to the Sun previously, and routinely studied indirectly by the 
radio emission they generate, however due to their low emission at optical wave- 
lengths this was the first time a shock was directly imaged and analysed at such 
large distances from the Sun. 

• The Mach number for the shock was found to be 4.1 it 1.2 and the shock stand-off 
distance was found to increase linearly to 2Oi?0 at 0.5 AU. 

• The normalised stand-off distance was found to be roughly constant with a mean 
value of 0.37 ± 0.09 D^^. The shock stand-off distance normalised by the CME 
radius agrees with the semi-emprical model (<%20). However the shock stand-off 
distance normalised by the CME radius of curvature differed significantly from the 
semi-emprical relations. This indicated that the CME radius of curvature was under 
estimated by a factor of ~ 3 - 8. This was found to be most likely due an observational 
effect, as only a cross-section of the CME is observed, thus we cannot estimate the 
most likely larger radius of curvature along the flux rope axis. 

• The radius of curvature of the CME along the flux rope axis was estimated to be 
~ 205 Rq (0.95 AU) at 1 AU, while the stand-off distance was estimated to be approx- 
imately 54 Rq . These values can be compared to multi point in situ measurements 
to gauge their validity. 



8.2 Future Work 

The methods developed and implemented in this thesis are the first to produce true CME 
kinematics at large distances from the Sun. Only a small number of events have been 
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analysed, and in order to rigorously answer the open questions on CME propagation a 
large number of events need to be analysed and statistical measures computed. In order 
to achieve this an event catalog needs to be compiled and studied. This is a very time 
consuming task which is subject to large user bias. Using advanced imaging processing 
it may be possible to automate or semi-automate this process. This would firstly remove 
the user bias, and also allow a much larger number of events to be analysed. Numerous 
methods have been used to extract 3D information from the STEREO observation. An 
inter-comparison of their accuracy, applicability (or lack thereof) in studying CMEs at large 
distance from the Sun would be of great value. More realistic modelling of the ambient 
solar wind would allow a better comparison or the drag models to the derived kinematics. 



8.2.1 Event Catalogue 

Analysing a small number of events has led to some interesting results. In order to rigor- 
ously make statistically significantly conclusions about what forces govern CMEs during 
their later propagation a larger number of events must be analysed in an identical man- 
ner. To facilitate this large scale analysis, software was written which, using only on start 
and end times for each instrument, could download, reduce and process the observations 
producing standardised quicklook movies and data (base difference, running difference, 
etc). Further software was written (xlmager) which allows the A and B observations to be 
analysed at the same time as shown in Figure |8.1[ This software also allows the scaling of 
the data, ellipse characterisation, and recording and editing of data points along the CME 
front. Two implementations of the tie-pointing algorithm have been written - one which 
uses the ellipse characterisations, and the other which interpolates between the data. The 
two methods are necessary as the front of the CME is not always well characterised by an 



ellipse, as shown in Figure 8.1 These methods take the output from the xlmager code and 
automatically extract the 3D trajectories and kinematics for the event, and store them for 
later use. 

Currently, a catalogue of 87 events between November 2007 and June 2011 has been 
compiled. Of these events, 15 have had a preliminary analysis performed. The remaining 
events have yet to be analysed. Once this is completed, analysis on the statistics of the 
kinematics can begin. Interesting questions like how many CMEs are accelerated, and 
by how much, could then be answered. Also, each event could be analysed using the 
bootstrap method to see if all slow and fast events show a different dependance on the 
difference between the CME and solar wind speed linear versus quadratic or if there is a 
distribution of values. A number of improvements could be made to the software, including 
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better handing of outliers in the kinematics, more robust 3D reconstruction, and better 
displaying and scaling methods to ease identification. 



8.2.2 Image Processing and Automation 



The faint nature of most CMEs is one of the key obstacles to studying their kinematics. 
The modified running difference technique was developed to enhance the CME signal. 



Various other image processing methods such as curvelets (Candes & Donoho, 1999 Starck 



et al. , 2002 ) have also been investigated. The idea of using curvelets is to represent a 



features as a superposition of functions of various lengths and widths obeying the width ^ 
length'^ scaling law. The first curvelets achieved by first decomposing the image into sub- 



bands or scales, and on each scale a local ridgelet Candes & Donoho (1999) transform is 



performed giving the curvelet decomposition. More recent curvelet implementations give 
the decomposition coefficients in continous frequency as: 



c{j,l,k) 



.1) 



at a scale 2 orientation I and position x^'' = Rg^^C^ -^ki, 2 ^^'^k2). The curvelet is obtain 
by rotation and translation of a mother curvelet ipj: 



'4'j,i,k{x) = ^j{Rei{x 



where Rqi is a rotation by 9i radians and 6i = 27r2 '■'/^l/. The waveform of ipj is defined in 
terms of its Fourier transform in polar coordinates: 



2b'/2]0^ 
27r ' 



(8.3) 



This type of multi-scale analysis has already been performed on coronagraph images ( Byrne 



et al., 2009: Gallagher et al., 2011) as shown in Figure 8.2 Figure 8.3 shows a curvelet 



decomposition of an HIl image, in which each image is a reconstruction using one scale. 
The curvelet transform gives much more information than just scale, such as direction and 
angle. By combining filtering across scales and properties it should be possible to enhance 
the signal of CME compared to the background corona and stellar field. 

Inspecting the HI images shows that there are clearly two size and intensity distributions 
present: small circular bright features (stars), and large some what amorphous blobs of 



very low intensity (CMEs) see Figure 8.4 (top row). Separating these two distributions 
seemed to lend itself to a two step process: firstly thresholding to perform the intensity 
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separation, and secondly, morphological operators to provide the size/shape separation. An 
initial threshold in intensity space gives a binary mask where all points below the threshold 
have been set to zero, and all points above set to unity. A morphological operator which 
draws conclusions on how a predefined shape (structuring element or kernel) fits or misses 
the shapes in the image is then applied to detect the desired property. Application of 
the morphological erode operation with a small circular kernel removes all small isolated 
features. The erode operation for binary data can be written: 

C = AeB= n (8.4) 
beB 

where A Q B is the erosion of image A by structure element B, and (^)-fe represents a 
translation of A by b. In the case of grayscale values the minimum difference of the set is 
used. Following this, a morphological dilate operator with a slightly larger circular kernel 



serves to form contiguous blocks and leaves only the CME and largest stars (Figure 8.4 
bottom left). The dilate operation for binary data can be written: 

C = A®B= U {A)i, (8.5) 

beB 

where A (B B is the dilation of image A by structure element B, and {A)fj represents a 
translation of A by b. In the case of grayscale values the maximum difference of the set 
is used. A contour operation allows the pixels corresponding to the front to be extracted. 



these can be over-plotted and the original image (Figure 8.4; bottom right). This technique 
proved to work relatively well, but the same threshold can not be used on an entire image 
sequence, and attempts to find a metric that gives a threshold automatically have thus 
far failed. However a variable threshold interactively changed by the user, in addition to 
selection of the appropriate contour, could provided an aid to the analysis. This type of 
semi-automation could be included in future versions of the xlmager software. 



8.2.3 3D Reconstruction 

Numerous methods have been developed to derive the CME positions from STEREO ob- 
servations. A detailed inter-comparison between all the available methods has not been 
done as yet. An especially interesting comparison would be between the two triangulation 
methods, tie-pointing (TP) and elliptical tie-pointing (FTP). It is known that in the TP 
method the reconstructed point will not lie on the CME, but either ahead of, or behind 
it, while for the FTP method the point should lie on the CME front. The difference be- 
tween the two methods has not yet been either quantitatively or qualitatively evaluated. 
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Studying the difference between the two methods as function of CME propagation angle 
and distance could be used to estimate these differences. Ideally, robust software to apply 
all the different methods to the extracted CME front would be developed, allowing an 
inter-comparison for any event studied. Understanding the effects of the various methods 
on the CME position, and hence kinematics, would enable better constraints on the errors 
used in the fitting process, thereby increasing confidence in the fit results. 

In the future a combination of forward modelling combined, with the 3D reconstruc- 
tions, will lead to a more complete picture of the CME observations in terms of the CMEs 



structure and its kinematics. It was shown in Section 3.1.1 that using single viewpoint ob- 
servations in combination with forward modelling was ambiguous. That is, the orientation 
of the model axis could not be uniquely determined, and, in most cases, two perpendicular 
orientations could both reproduce the data adequately. Using two viewpoints of STEREO 
should remove this ambiguity. Some work has already been done in this area. Figure 



shows synthetic coronagraph images of a graduated cylindrical shell (GSC) model as ob- 
served by A and B (Mierla et al. , 2009). Simultaneously constraining the fit from A and B 
observations should provide a unique best fit. There are many benefits of using a forward 
model like the GSC, for example better interpretation of the observed structures and more 
detailed kinematics can be extracted. It is well known that the velocity of a CME front 
is a combination of the CME velocity and its expansion velocity: using a model like the 
GSC the centre of mass velocity and expansion velocity can be derived. When comparing 
the CME kinematics to models the centre of mass velocity is the preferable measure. 



8.2.4 Kinematic Modelling 

The first obstacle in kinematic modelling is identifying and tracking the CME, the next 
problem is the accurate determination of its kinematics. Given height-time data, the first 
and second derivatives with respect to time need to be calculated, giving the velocity 
and acceleration respectively. As the data is not continuous, numerical approximations 
to the derivative operator must be used. Traditional methods for determining numerical 
derivatives such as forward, reverse, and centred difference have serious limitations. One 
serious limitation is the calculation of errors on the derivatives which depend inversely 
on the step size - linearly for forward and reversed differences and quadratically for the 
centred difference ~ double the candance and the error is quadrupled. This is extremely 
counter intuitive: as the number of data points is increased the errors also increase, even 
though more information is being added. More advanced methods like 3-Point Lagrangian 
interpolation suffer from the same problem, an error that depends inversely on step size. 
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The bootstrap technique provides an accurate and robust way to test if the kinematics fit a 
profile or model, but requires that these be known beforehand. Some methods exist which 
may allow better kinematics to be obtained using only the data itself, such as the Kalman 
filter and inversion methods. These methods should be investigated in the future to see if 
they perform better than the current methods. 

Once accurate kinematics have been derived, the next step is to compare them to 
theoretical models such as the drag model. In the drag model a number of the unknown 
parameters were absorbed into a power law parameterisation. Any knowledge of these 
parameters would allow a more accurate test of the drag model. Observation based models 
such as ENLIL and HAFv2 give the ambient solar wind conditions in the inner Heliosphere. 



Figure 8.6 is a snap shot in time from an ENLIL simulation showing solar wind density 
and velocity at various cuts through the simulation domain. Once the CME trajectory is 
known the solar wind properties such as density and speed can be extracted along the CME 
path. These can then be input into the drag models, which would be a great improvement 
to the current analytic solar wind models used. The ambient solar wind conditions have 
been shown to play a vital role in determining CME kinematics. Better constrains on 
these conditions would allow a better comparison between the observed kinematics and 
drag model predictions. 

Also ENLIL and HAFv2 allow CMEs to be simulated. The CME are simulated by time 
dependent changes to boundary condition (density, temperature, velocity). Comparing the 
ID drag models to these complex 3D simulations would allow a better understanding of 
which forces play the most important role and where they do so. This could give insight 
into which type of drag force is acting and the model best suited to study it. 
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Figure 6: Filtered raw images using an isotropic wavelet (left) and curvelets (right) by 
removing coefficients most probably due to noise. 
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Figure 7: Contrast enhanced raw images by amplifying the isotropic wavelets coefficients 
(left) and the curvelet coefficients (right) at the finer scales. 
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Figure 8.2: Corongraph Images Filter using an Isotopic Wavelet (left) and a Curvelet right 
(IGallagher et oLl 120111). 
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Figure 8.3: Curvelet decomposition of a HIl image showing different decomposition scales. 
The units are image pixels. 
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Figure 8.4: Original HIl image scale to 95% of the maximum intensity (top left) and the same 
image scaled to 0.5% of the maximum intensity (top right). Resulting image after application of 
a morphological open and close operators with different size kernels (bottom left) the detected 
CME front over plotted (bottom right). 
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Figure 8.5: Synthetic Graduated Cylindrical Shell STEREO Observations from the perspec- 



tive of A and B (Mierla et al. 2009) 
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Figure 8.6: ENLIL simulation results showing the velocity and density from NASA's Com- 
munity Coordinated Modeling Center (CCMC). 
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Propagation of Coronal Mass Ejections in 
the Inner Heliosphere 

[Shane A. Maloney, B.Sc. (Hons.)| 

Solar Coronal mass ejections (CMEs) are large-scale ejections of plasma and magnetic 
field from the corona, which propagate through interplanetary space. CMEs are the most 
significant drivers of adverse space weather on Earth, but the physics governing their 
propagation through the Heliosphere is not well understood. This is mainly due to the 
limited fields-of-view and plane-of-sky projected nature of previous observations. The 
Solar Terrestrial Relations Observatory (STEREO) mission launched in October 2006, was 
designed to overcome these limitations. 

In this thesis, a method for the full three dimensional (3D) reconstruction of the tra- 
jectories of CMEs using STEREO was developed. Observations of CMEs close to the 
Sun (<15i?0) were used to derive the CMEs trajectories in 3D. These reconstructions 
supported a pseudo-radial propagation model. Assuming pseudo-radial propagation, the 
CME trajectories were extrapolated to large distances from the Sun (15-240 -R0). It was 
found that CMEs slower than the solar wind were accelerated, while CMEs faster than the 
solar wind were decelerated, with both tending to the solar wind velocity. 

Using the 3D trajectories, the true kinematics were derived, which were free from pro- 
jection effects. Evidence for solar wind (SW) drag forces acting in interplanetary space 
were found, with a fast CME decelerated and a slow CME accelerated toward typical SW 
velocities. It was also found that the fast CME showed a linear dependence on the velocity 
difference between the CME and the SW, while the slow CME showed a quadratic depen- 
dence. The differing forms of drag for the two CMEs indicated the forces responsible for 
their acceleration may have been different. Also, using a new elliptical tie-pointing tech- 
nique the entire front of a CME was reconstructed in 3D. This enabled the quantification 
of its defiected trajectory, increasing angular width, and propagation from 2 to 46i?0 (0.2 
AU). Beyond 7 Rq, its motion was shown to be determined by aerodynamic drag. Using 
the reconstruction as an input for a 3D magnetohydrodynamic simulation, an accurate 
arrival time at the LI Lagrangian point near Earth was determined. 

CMEs are known to generate bow shocks as they propagate through the corona and SW. 
Although CME-driven shocks have previously been detected indirectly via their emission at 
radio frequencies, direct imaging has remained elusive due to their low contrast at optical 
wavelengths. Using STEREO observations, the first images of a CME-driven shock as it 
propagates through interplanetary space from 8 Rq to 120 Rq (0.5 AU) were captured. The 
CME was measured to have a velocity of ~ 1000 km s~^ and a Mach number of 4.1 it 1.2, 
while the shock front standoff distance (A) was found to increase linearly to ~ 20 Rq at 0.5 
AU. The normalised standoff distance (A/Dq) showed reasonable agreement with semi- 
empirical relations, where Dq is the CME radius. However, when normalised using the 
radius of curvature {A/Rq), the standoff distance did not agree well with theory, implying 
that Ro was underestimated by a factor of ~ 3-8. This is most likely due to the difficulty 
in estimating the larger radius of curvature along the CME axis from the observations, 
which provide only a cross-sectional view of the CME. The radius of curvature of the CME 
at 1 AU was estimated to be ~ 0.95 AU 



